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Infrared  signature  models  covering  the  3.5  to  5-micrometer  spectral 
region  have  been  developed  for  events  likely  to  produce  false  alarms 
in  airborne  IR  warning  receivers  defending  against  SAM  and  AAM  threats. 
The  models  can  be  used  to  predict  the  spectral  radiation  of  clouds, 
ocean  sunglitter,  uniform  specular  sunglints,  hot  gas  effluents  repre- 
senting industrial  emissions  and  threat  plumes,  and  the  apparent  spectral 
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intensities  of  hot  graybody  sources.  The  spectral  resolution  is 
limited  by  the  radiative  transfer  codes  utilized  (LOWTRAN  4 and  ATLES) 
except  in  the  blue  spike  region  (2380-2400  reciprocal  centimeters),  for 
which  a new  pseudo  line-by-line  procedure  was  developed  yielding  fully- 
resolved  spectra.  All  of  the  models  have  been  successfully  validated 
against  measurements  and/or  results  of  more  elaborate  theoretical 
calculation.  Included  in  this  report  are  new  AFGL  and  EPA  spectral 
data  on  clouds,  sunglints  and  industrial  emissions. 


PREFACE 


Useful  discussions,  data  and  computer  codes  were  supplied  during  this 
effort  by  Eric  Shettle  and  Brian  Sandford  of  the  Air  Force  Geophysics  Laboratory, 
Warren  Wiscombe  of  the  National  Center  for  Atmospheric  Research,  Horace  Dry  of  R&D 
Associates,  William  Herget  of  the  Environmental  Protection  Agency,  Frederick 
Simmons  and  Stephen  Young  of  the  Aerospace  Corporation,  and  William  Keslar  and 
G.  Behrens  of  U.S.  Steel  Corporation.  Fred  Smith  and  Paul  Eitner  of  Science 
Applications,  Ann  Arbor,  MI,  performed  the  work  reported  in  Sections  6 through  8 
under  Lub contract.. for.. Honey well;  ^hey  also  assisted  in  developing  a code  for  the 
blue  spike  model  described  in  Section  5.  Richard  Sanderson  of  the  Air  Force 
Avionics  Laboratory  was  project  monitor  for  the  study. 
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SECTION  I 
INTRODUCTION 


Detection  of  valid  threats  and  attainment  of  very  low  sensor  false  alarm 
rate  remain  essential  objectives  for  infrared  threat  warning  receivers. 
Conventional  single-color  systems  using  spatial  and  temporal  discriminants  have 
failed  to  demonstrate  adequate  protection  commensurate  with  insensitivity  to 
false  threats.  It  is  now  generally  recognized  that  superior  warning  receiver 
performance  may  be  attainable  through  spectral  discrimination.  However,  there  has 
not  existed  an  adequate  data  base  or  analytical  standard  of  evaluation  for 
measuring  the  performance  improvement  to  be  realized  through  multispectral 
discrimination.  These  deficiencies  are  addressed  in  two  study  contracts  awarded 
by  the  Air  Force  Avionics  Laboratory:  the  IR  Signature  Studies  documented  in  this 
report  and  the  IR  Receiver  Analysis  in  progress  at  Nichols  Research. 

The  long-term  objective  of  these  two  studies  and  the  planned  IR  Signa- 
ture Studies  II  procurement  is  to  develop  for  the  Air  Force  a performance 
assessment  model  applicable  to  warning  receivers  employing  a wide  range  of 
spectral  discrimination  techniques.  This  capability  will  allow  the  effectiveness 
of  any  existing  or  proposed  technique  to  be  estimated  from  computer  analysis 
rather  than  costly  hardware  development  and  testing. 

The  purpose  of  the  IR  Signature  Studies  has  been  to  identify  events 
which  constitute  potential  false  threats  in  the  3.5  to  5-micrometer  spectral 
region,  and  to  select,  develop  and  validate  models  to  represent  their  spectral 
signatures  at  resolutions  sufficient  to  evaluate  both  spectral  line  correlation 
and  broad  band  discrimination  techniques.  The  false  threats  of  interest  are  those 
which  could  be  encountered  by  airborne  sensors  operating  in  the  5 kft  to  30  kft 
altitude  range  and  defending  against  AAM  and  SAM  threats.  It  was  a requirement  of 
the  study  that  the  models  be  capable  of  predicting  the  signatures  at  essentially 
infinite  spectral  resolution  in  the  blue  spike  region  of  the  real  threats,  since 
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potential  discrimination  techniques  could  utilize  the  regular  line  structure  of 
the  blue  spike  threat  emission.  Elsewhere  in  the  3.5  to  5— micrometer  spectral 
region,  lower  resolution  band-model  techniques  were  deemed  adequate,  with  the 
possible  exception  of  regions  in  which  regularly-spaced  HC1  or  CO  emission  lines 
offer  additional  discrimination  potential.  The  contract  called  for  similar 
modeling  of  AAM  and  SAM  spectral  signatures  based  on  simple  generic  models  that 
assume  homogeneous  isothermal  plumes.  The  study  emphasizes  false  threat  spectral 
signatures  as  compared  to  threat  signatures. 

The  model  validation  task  involved  the  collection  of  measured  spectral 
signatures,  and  a requirement  of  the  study  was  to  identify  to  the  Air  Force 
Geophysics  Laboratory  specific  scenes/locations  which  should  be  measured  as  part 
of  the  AFGL  Background  Measurement  Program  (Sandford  et  al^).  Generally  the  type 
of  data  asked  for,  and  supplied  by  AFGL,  was  not  available  in  the  literature  or 
from  other  sources.  The  AFGL  spectra,  representing  mostly  clouds,  sunglints  from 
lakes  and  ponds,  and  industrial  site  emissions,  were  extremely  useful  in  the  model 
validations.  Many  of  these  are  reproduced  in  this  report,  since  AFGL  does  not 
plan  to  publish  them. 

Section  2 summarizes  the  study  results  and  conclusions.  The  remaining 
sections  and  Appendixes  define  the  signature  models  selected  and  developed  in  the 
study,  and  the  validations  of  the  models. 
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SECTION  II 


SUMMARY  AND  CONCLUSIONS 

This  report  describes  a set  of  models  that  accurately  characterize  the 
spectrally  resolved  infrared  signatures  of  AAM  and  SAM  threat  models,  and  of 
background  events  likely  to  cause  false  alarms  in  airborne  IR  threat  warning 
receivers.  The  general  prescriptions,  equations,  identifications  of  utilized 
(available)  codes,  and  literature  references  constitute  a fairly  complete  defini- 
tion of  the  models  and  the  procedures  required  in  their  implementation. 

The  potential  false  threats  identified  in  the  study  are: 

Clouds  (solar  scattering  by). 

Ocean  sunglitter, 

Protected  small  lakes  and  pools,  swamps,  discharge  basins, 
artifacts  (specular  solar  reflections  by), 

Industrial  sites: 

Flares , 

Stack  plumes, 

Refinery  emissions, 

Hot  solid  objects, 

Fires : 

Incendiary  fires, 

Refuse  fires  (battlefield,  industrial,  agricultural), 

Brush  and  forest  fires, 

Munitions: 

Muzzle  flash, 

Shell  bursts, 

Flares , 

Munition-induced  explosions, 

Battlefield  smokes, 

Vehicular  emissions. 
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The  model  that  predicts  cloud  spectral  signatures  is  the  result  of  an 
extensive  mathematical  analysis  that  produced  new  approximate  solutions  of  the 
radiative  transfer  equation.  The  model  is  fairly  general,  and  its  predictions  are 
in  excellent  agreement  with  AFGL  measurements. 

The  ocean  sunglitter  model  represents  a generalization  of  the  well- 
known  Cox-Munk  model  to  include  spectral  variations  calculated  from  the  Fresnel 
spectral  reflection  coefficients  for  a clean  water  surface.  A simplified  version 
of  the  model  gives  the  spectral  signatures  of  specular  sunglints  from  smooth  water 
surfaces.  Both  versions  were  successfully  validated  against  AFGL  measurements. 

A new  model  and  computer  code  (BLUESPIKE)  were  developed  to  calculate 
the  fully-resolved  spectrum  of  the  blue  spike  region  of  threats  and  false  threats. 
Line-by-line  and  band-model  techniques  were  combined  in  the  model  to  achieve  a 
computational  efficiency  much  higher  than  is  typical  in  conventional  line-by-line 
calculations.  Very  little  accuracy  was  sacrificed,  as  demonstrated  by  comparisons 
to  line-by-line  calculations  and  high  resolution  measurements. 

A substantial  part  of  the  study  effort  consisted  of  acquiring, 
modifying,  and  testing  radiative  transmission/emission  computer  codes  needed  in 
the  cloud  and  water  reflection  models,  and  needed  in  generating  signatures  of 
gaseous  threats/false  threats  outside  the  blue  spike  region.  We  developed  a 
modified  version  of  LOWIRAN  4 that  is  applicable  to  the  folded-path  geometries 
involved  in  solar  reflections  from  clouds  and  water  surfaces.  We  extended  the 
capabilities  of  the  ATLES  band-model  code,  validated  the  code  against 
measurements,  and  increased  its  efficiency. 

With  respect  to  gaseous  emission  sources,  we  have  concluded  that  the 
procedures  used  in  this  study  are  adequate  for  modeling  events  likely  to  cause 
false  alarms  in  warning  receivers.  These  validated  procedures  are: 

a)  Use  of  ATLES  for  moderate  resolution  hot  gas  emission 

b)  Use  of  BLUESPIKE  for  high  resolution  gaseous  emission 
in  the  blue  spike  region,  and 
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c)  Use  of  LOWTRAN  4 for  low  resolution  graybody  emission 
sources. 

While  generally  adequate  for  most  warning  receiver  evaluations,  the  ATLES  and 
LOWTRAN  codes  are  of  course  limited  in  resolution  and  thus  not  useful  for  modeling 
discrimination  schemes  requiring  high  spectral  resolution. 

In  the  course  of  the  study  a number  of  possible  discrimination  tech- 
niques were  identified.  In  general,  these  techniques  rely  on  inference  of 
physical  source  characteristics  from  the  observed  spectra.  The  potentially  useful 
inferences  and  methods  demonstrated  include: 

1)  Temperatures  of  graybody  sources  (from  3.5  to  4.2 
continue  emission) 

2)  Identification  of  CC^  gaseous  emission  (from  the  blue 
spike) 

3)  Thickness  determination  for  C0„  plume  (from  red  wing 

emission)  2 

4)  Identification  of  various  gaseous  species  in  the 
plume  including  SO  HC1  and  CO  (from  spectral 
emission  lines). 

In  the  order  given,  each  inference  requires  higher  resolution  instrumentation, 
varying  from  40  cm  for  a temperature  determination  to  1 cm-1  for  determination 
of  CO  or  HC1  in  the  source. 

A detailed  analysis  of  discrimination  techniques  has  not  been  performed 
here;  however,  it  does  appear  that  the  above  methods  will  be  useful  for  dis- 
crimination of  aircraft  or  missiles  from  false  threats. 

We  could  not  find  models  (or  adequate  data  for  the  development  of 

models)  for  predicting  the  spectral  signatures  of  fires,  munitions-related  events 
cr  battlefield  smokes.  The  data  collected  on  these  sources  is  given  in  Sections  7 
and  8. 


SECTION  III 
CLOUD  RADIANCE  MODEL 

The  cloud  radiance  model  developed  in  this  study  predicts  the  infrared 
spectral  radiance  field  above  thick  sunlit  clouds.  It  can  probably  be  proved  that 
the  maximum  reflected  radiance  above  a cloud  exceeds  or  equals  the  maximum 
diffusely  transmitted  radiance,  and  it  is  a fact  that  the  reflected  radiance 
increases  with  the  thickness  of  the  cloud.  In  this  sense  the  model  predicts  the 
brightest  IR  spectral  signatures  to  be  expected  from  clouds.  The  model  could  be 
generalized  to  deal  with  arbitrarily  thin  clouds  and  the  transmitted  radiance 
field  by  working  out  additional  solutions  according  to  the  methods  defined  in 
Section  3.2. 

The  model  has  three  parts.  The  first  is  an  approximate  solution  of  the 
integro-dif ferential  radiative  transfer  equation  describing  the  reflected 
radiance  field  in  terms  of  the  incident  flux  and  the  aerosol  scattering  parameters 

of  the  cloud.  The  solution  was  obtained  by  applying  the  double-delta  approxi- 

. 2 ... 
raation  of  Turner  to  a nonconservative  scattering  medium;  Turner  treated  the 

conservative  case  only.  The  solution  can  be  used  to  predict  atmospheric  back- 

scatter  (for  a homogeneous  atmosphere)  as  well  as  cloud  scattering.  This  part  of 

the  model  is  described  in  Sections  3.1  through  3.3. 

The  second  part  of  the  model  is  an  extensive  data  set  describing  the 
basic  scattering  and  extinction  properties  of  different  size  distributions  of 
water  droplets  and  ice  crystals.  The  data  consists  of: 

Kq(v),  the  extinction  cross  section  per  particle 

Uq(v),  the  single  scattering  albedo  (scattering  coef- 
ficient t sum  of  absorption  and  scattering  coef- 
ficients) 


g(v), 


the  scattering  phase  function  asymmetry  param- 


The  Mie  series  computations  required  to  obtain  these  parameters  was  performed  by 
Eric  Shettle  of  AFGL  using  the  log-normal  size  distribution: 


dN(r) 

dr 


(log  r - log  r ) 


where  o = 0.3,  and  r is  the  mean  radius  of  the  particles  (the  computations  assume 

m 

spherical  water  droplets  and  spherical  ice  "crystals"}  The  wat^r  and  ice  indices 
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of  refraction  were  taken  from  Hale  and  Querry  . Dr.  Shettle  provided  10  sets  of 
Kq(v) , Wq(V)  and  g(v)  corresponding  to  5 different  sizes  of  water  droplets  and  ice 
crystals.  In  each  set,  which  covers  the  spectral  region  2.6  to  5-micrometer,  the 
wavenumber  spacing  is  sufficient  to  define  completely  the  spectral  variations  in 
the  three  quantities.  The  data  are  given  in  Appendix  A.  The  scattering  phase 
function  was  not  computed  since  this  would  have  involved  a considerable 
expenditure  of  computer  time  for  the  number  of  wavelengths  involved.  This  means 
that  a phase  function  with  g as  an  adjustable  parameter  must  be  assumed  (see 
Section  3.3). 

The  third  part  of  the  model  is  a modified  version  of  the  LOWTRAN  4 
computer  code.^  This  part  of  the  model  accounts  for  atmospheric  absorption  of  the 
incident  and  reflected  radiation,  for  atmospheric  absorption  of  thermal  radiation 
emitted  by  the  cloud,  and  for  atmospheric  emission. 

Figure  1 shows  the  transmission/emission  geometry  used  in  the  mod- 
ified LOWTRAN  code.  Point  A is  the  location  where  the  cloud  spectral  radiance  is 
to  be  determined.  The  code  computes  the  spectral  transraittances  of  paths  SCA  and 
c.k,  and  the  spectral  radiance  of  path  CA.  It  also  computes  the  blackbody  spectral 
emission  of  the  cloud,  multiplies  it  by  the  transmittance  of  path  CA,  and  adds  the 
result  to  the  spectral  radiance  of  path  CA  to  obtain  the  thermal  contribution  to 
the  predicted  spectral  signature. 
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Figure  1.  Geometry  Used  in  Modified  LOWTRAN  4 Code 

The  solar  spectrum  reflected  by  the  cloud  in  a given  direction  is 
performed  according  to  the  calculations  described  in  Section  3.2.  This  spectrum 
is  multiplied  by  the  spectral  transmittance  of  path  SCA,  and  added  to  the  thermal 
spectrum  computed  by  the  modified  LOWTRAN  code.  Note  that  the  model  does  not 
include  effects  of  scattering  by  the  atmosphere  outside  the  cloud  (except  for 
extinction),  reflection  of  atmospheric  thermal  radiation  by  the  cloud,  scattering 
of  cloud  thermal  emission  within  the  cloud,  or  molecular  absorption  inside  the 
cloud.  These  processes  are  assumed  to  have  relatively  small  effects  in  the  cloud 
spectral  radiance  over  the  3 to  5-micrometer  region. 

3.1  DEFINITIONS  AND  FORMAL  SOLUTION 

Figure  2 defines  the  geometry  used  to  characterize  the  radiance  field 
and  the  scattering  layer  (atmosphere  and/or  cloud).  The  scattering  layer  is 
assumed  to  be  plane-parallel  and  homogeneous.  By  homogeneous  we  mean  that  the 
aerosol  scattering  properties,  namely  the  single  scattering  albedo  and  phase 
function,  are  independent  of  height  z or  optical  thickness  t.  Optical  thickness 
is  one  of  the  independent  variables  used  to  define  the  radiance  field;  hence,  the 
variation  t(z)  and  the  vertical  distribution  of  the  density  of  scatterers  is 
arbitrary . 
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Solar  radiation  is  incident  at  t * 0,  and  the  total  optical  thickness  of 
the  atmosphere  and/or  cloud  is  Tg.  The  coordinates  of  the  sun  will  be  denoted  by 
Pq  and  4>g»  where  pg  is  the  cosine  of  the  solar  zenith  angle  and  $q  is  the  solar 
azimuth  direction.  The  field  coordinates  of  direct  solar  or  diffuse  (scattered) 
radiation  are  (t,  p,  <P)  with  the  convention  that  positive  p ■ cos  0 corresponds  to 
"upwelling"  radiation  and  negative  p to  "downwel ling"  radiation.  Thus,  the  sun's 
rays  or  a diffuse  stream  parallel  to  the  sun's  rays  has  the  direction  —Pg,  <f>g>  and 
a stream  towards  the  sun  has  the  direction  Pg,  $g  + tt.  We  will  refer  to  these  two 
directions  in  calculating  the  diffuse  radiance  I (t,p  <j>)  > but  the  final  result 
will  represent  only  the  upwelling  component;  i.e.,  p will  be  restricted  to 
positive  values. 

As  mentioned  earlier,  emission  by  the  atmosphere/cloud  is  handled 
separately,  with  the  assumption  that  there  is  no  scottering  of  the  emitted 
radiation. 


The  radiative  transfer  equation  for  the  diffuse  radiance  I (i,p,<f>)  in  a 
plane-parallel,  one-dimensional,  homogeneous,  nonemitting  atmosphere  is: 


dl  / a', 

P ( t » P , 4>) 


I ( T , U , ♦) 


p(p,<1>,p\  <t>' ) 


I(t,p'  , <t>'  )dp  ' d4> ' 


Es(t)  p(p,<f>,-Pg,<t>0) 


(1) 


where  ujg  is  the  albedo  for  single  scattering  (ratio  of  scattering  coefficient  to 
sum  of  absorption  and  scattering  coefficients), 


B (z)  _ 8(z) 

0 ct(z)  + 8 (z)  ~ k ( z ) 


constant, 


(2) 


and  t is  the  optical  thickness  of  the  atmosphere  above  level  z; 
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k(z ' )dz ' 


The  constant  funi^or/pfp^ ) is  the  single-scattering  phase  function  for  scattering 
angle  pg  ■ cos  8g)  which  is  defined  by  the  directions  of  the  incident  ( y * > <J>' ) and 
scattered  ( p , 4>)  rays  & t is  normalized  according  to: 


p(cos  0) 


dw(9) 


The  second  term  of  Eq.  (1)  represents  the  scattering  of  direct  solar  flux 
(irradiance)  Eg(r)  = Eq  exp  (— t/Uq)  which  has  penetrated  to  level  r;  Eq  denotes 
the  extraterrestrial  solar  irradiance  on  a plane  normal  to  the  sun's  rays. 


By  multiplying  Eq.  (1)  by  the  integrating  factor  expf-t/u),  dividing  by 
P and  integrating  from  Tq  to  t,  we  get,  for  the  upwelling  radiance 


I(t,M)  - I(t0,vi,4>)  exp  (tq~t)/p 


2 it  1 


/o  -1  *4 


j-(t,-t)/pi(t,  |Wi^l)dTldyl<1^ 


4ir(p+p 


Vo  Eo  p(M.-p0*V  "t/Pq  "To/yo  -(t0't)/p 


By  integrating  Eq.  (l)  from  0 to  t we  could  also  obtain  an  expression  for  the 
dovnwelling  radiance,  but  as  stated,  we  will  obtain  a final  solution  for  the 
upwelling  case  only.  We  emphasize  that  I(t,p,  ) in  Eq.  (1)  and  in  the  integral 
term  in  Eq.  (5)  represents  the  radiance  in  any  direction. 


Our  objective  is  to  reduce  Eq.  (5)  tc  an  approximate  closed-form 
solution  using  a technique  described  by  Turner  . This  will  be  accomplished  by 
introducing  approximations  in  Eq.  (1)  which  reduce  it  to  a pair  of  first-order 
differential  equations,  whose  solution  yield  a first  approximation  to  I(t,u,4>). 
This  result  will  then  be  substituted  into  Equation  3.5  to  obtain  a second 
approximation  to  I(t,u,<J>).  The  first  approximation  is  basically  a two-stream 
model  that  involves  the  diffuse  fluxes  (irradiances)  in  .he  positive  and  negative 
z directions,  which  are  defined  by 


E+(t) 


y I(T,y,4>)dpd<t>, 


(6) 


and  E (t) 


U I ( T ,y  ,4>)dud<|>, 


respectively.  The  corresponding  upward  and  downward  average  diffuse  radiances  are 


<I(t)>± 


I(t  ,y , <j>)dyd<{> 


(7) 


3.2  APPLICATION  OF  THE  DOUBLE-DELTA  APPROXIMATION 

We  will  outline  the  steps  in  deriving  the  radiance  field  in  the  double- 
delta approximation  and  will  write  out  only  a few  of  the  intermediate  results.  We 
will  present  the  final  result  in  a simplified  form  applicable  to  optically  thick 
clouds , 


The  method  derives  its  name  from  the  form  used  to  represent  the  phase 
function  in  the  first  approximation.  We  assume  initially  that  the  phase  function 
consists  of  forward  and  backward  delta  functions;  i.e., 
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p(p,<j>,p'  ,4)' ) - 4iTT)6(»j-y'  )6($-<J>’ ) 


+ 4TT(l-n)6(y+y ' )6(4>-^,+1t) , (8) 

where  n is  the  fraction  of  radiation  scattered  into  a forward  hemisphere. 
Consistent  with  this  representation  of  the  phase  function,  we  initially  approx- 
imate the  diffuse  radiance  field  by  delta-functions  directed  towards  and  away  from 
the  sun: 


I(T,y,<t>) 


E+(t)6(p-m0)6(4>  -i|>0+it) 


+ E_(T)<$(p+y0)<S(<|>-4>0)  , 


(9) 


where  E+(x)  is  defined  by  Eq.  (6).  Substituting  Eq.  (8)  into  Eq.  (1)  we  obtain 

v ■ ( l-o>0n) i ( t , y ,<$> ) - w0(i-n)i(T,-n4+n) 


-OLE  (t) 
0 s 


n6(M+w0)6(<t>-$0)  + ( 1-J1)S(p-vi0)6(<J>-4>q+iO 


(10) 


Substitution  of  (9)  into  (10)  yields  the  coupled  linear  differential  equations 


dE_ 

dr 


(l-u^n) 


u»0(i-n) 


E+(t) 


E+(  T ) 


U)  (1-0) 

~ E (t)  - «-(l-0)E  (t) 

U0  - 0 s 

(1-0)  0) 

e (t)  + o>„nE  (t). 

- 0 s 


(11) 


If  the  albedo  of  the  surface  at  x = Tq  is  zero,  or  if  the 
atmosphere/cloud  has  sufficiently  large  optical  thickness,  one  of  the  boundary 
conditions  on  Eq . (11)  is 


(12) 


The  other  boundary  condition,  which  always  applies,  is 

E_(0)  “ 0;  (13) 

i.e.,  the  downward  diffuse  irradiance  is  zero  at  the  upper  boundary  (t  = 0).  The 
solutions  of  Eq.  (11)  for  these  boundary  conditions  will  be  denoted  E^(t)  and  E' 
(t).  After  obtaining  these  solutions,  we  can  relax  the  condition  that  the  surface 
albedo  equal  zero.  That  is,  we  can  obtain  a second  solution,  denoted  E^(t)  and 
E"(t),  which  assumes  no  solar  input,  and,  via  appropriate  boundary  conditions, 
represents  the  reradiation  of  (reflection  by)  the  lower  boundary.  The  two 
solutions  can  be  superimposed  to  obtain  the  total  irradiance,  and  then  the  total 
radiance . 


With  the  boundary  conditions  (12)  and  (13),  the  solution  of  (11)  is 


P0E0w0(l-n)  sinh 


e;(t) 


E'(t) 


X(t0“t) 


(l-<^n)  sinh  ATq  + PqX  cosh  Xt^ 


P0E0 


(I-WqTI)  sinh  [X(VT)]  + pQX  cosh  [x(t0“t)J 
(l-0)on)  sinh  XTq  + PqX  cosh  Xt^ 


(14) 

-t/w. 


-e 


where 


X - (i/v0) 


|u-<^n)2  - o>02  d-n)2J 


1/2 


(15) 


Note  that  X * 0 for  Wq  ■ 1;  for  this  case  Eq.  (14)  must  be  evaluated  in  the  limit 
as  X >0.  If  (12)  and  (13)  are  the  appropriate  boundary  conditions,  then  the 
upwelling  radiance  in  the  second  approximation  is  obtained  by  substituting  (14) 
into  (9)  and  then  (9)  into  (5).  This  solution  is  given  in  Section  3.3 
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When  the  underlying  surface  is  diffuse  and  has  non-zero  albedo  p,  the 
solution  for  the  upwelling  radiance  can  be  obtained  as  follows:  we  assume  no 
solar  input  and  an  isotropic  source  at  the  bottom  of  the  atmosphere,  representing 
the  diffusely  reflected  radiation.  It  is  reasonable  under  these  conditions  to 
assume  that  the  radiance  field  is  nearly  isotropic,  i.e.,  that 


I(T,M,<j>)  cs  <I(t)>±  =*  2E±(t)/21T, 


where  the  second  relationship  follows  from  the  definitions  3.6  and  3.7.  Note 
that,  in  general, 

<yl(r)>  = E (t)/27T  (17) 


according  to  these  definitions.  Thus,  by  averaging  Eq.  (10)  first  over  the  upward 
hemisphere  and  then  over  the  downward  hemisphere,  and  setting  Eg  ” 0,  we  obtain 
the  pair  of  equations 

dE 

- 2(l-(U0n)E+(T)  - 2to0(l-n)E_(r) 

dE_ 

~=  2u)0(i-n)E+(r)-  2(i-w0n)E_(T). 


The  boundary  conditions  are 


E.<V  ■ p 


E (0)  ■ 0 


fc'V  * “oEs<Vl  ■ pk(’o>  ♦ EoVV1 

L -I  / a 


where  E^.(t)  is  the  previously  obtained  solution  for  p ■ 0 (second  of  Eq.  (14)). 

Let  E|(x)  denote  the  solutions  of  Eq.  (18)  with  the  boundary  conditions 
(19).  Then  the  average  value  of  the  nearly-isotropic  radiance  field,  according  to 
(16),  is 
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(20) 


1 

2 


[<I(t»+  ♦ <I(t)>_ 


+ 


FV(t) 


/2tt 


Eq.  (9),  with  E|(t)  in  place  of  E+(t),  gives  the  solar-generated  anisotropic 
field,  in  the  first  approximation.  The  sum  of  this  field  and  the  one  arising  from 
surface  reflection  is 


i(t,p,4>)  - -J-  Je; 

+ E^(r)6(p+P0)6(<fi-4>0) 
+ E^'(t)  + E”(t)J  /2tt. 


( t )6  (y-yQ)  6(<^ -^0+1T) 


(21) 


This  result  would  be  substituted  into  Eq.  (5)  to  obtain  the  second  approximation 
to  the  upwelling  radiance.  Because  of  the  forms  of  (21)  and  the  solutions  E'  and 
E",  the  integrations  can  be  performed  to  yield  a closed-form  result,  without 
specification  of  the  scattering  phase  function.  Thus,  any  appropriate  phase 
function  may  be  used  in  place  of  the  double-delta  function  given  by  Eq.  (8). 


3.3  SOLUTION  APPLICABLE  TO  OPTICALLY  THICK  CLOUDS 


Very  little  diffused  sunlight  emerges  from  the  undersides  of  thick  water 
or  ice  clouds  in  the  3 to  5-microjeter  spectral  region,  where  the  single- 
scattering albedo  is  always  appreciably  less  than  unity.  Hence,  it  is  reasonable, 
for  the  purposes  of  this  study,  to  assume  that  the  surface  albedo  is  zero,  and  to 
utilize  the  simpler  solution  that  results  from  this  assumption. 


For  P ■ 0,  the  second  approximation  to  the  upwelling  radiance  I(t,v,  ) 
is  obtained  by  substituting  (14)  into  (9)  and  then  (9)  into  (5).  Note  that  the 
first  term  of  (5)  is  zero  when  p 3 0.  Evaluation  of  the  integrals  in  (5)  yields 
the  desired  result: 


-ri^rn! 


“0EOrt’ 


sinh  Xtq  + PqX  cosh  Xtq 


where 


P(ys) 
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- cosh  X(tq-t))  + sinh  X(Tq-t) 


. 2,2 

1-y  X 


I ‘-“o”  A . 
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inh  X(Tq-t) 
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fraction  of  radiance  scattered  by  the  aerosols 
into  a forward  hemisphere 
single-scattering  albedo 

total  optical  thickness  of  atmosphere/c loud 
extraterrestrial  solar  irradiance  (or  spectral 
irradiance)  on  plane  normal  to  the  sun's  rays 

<‘V  [(131^?o2<1-'')2]11/2 

p yu0  -^/(l-y  )(l-y^)cos(i)i-(t>n)j 
pjjPU0  + ■y/(l-y2)(l-y^)cos(<J)-(})l))J 
single-scattering  phase  function 
cosine  of  solar  zenith  angle 
azimuth  of  sun 

cosine  of  zenith  angle  of  upwelling  radiance 

(0  < v < 1) 

azimuth  of  upwelling  radiance 

optical  thickness  (level)  at  which  upwelling 
radiance  is  observed. 


Note  that  the  units  of  I(t,p,$)  are  those  of  Eq  times  sr  . The  arguments  of  the 
phase  functions  P.  and  P.  are  simply  the  cosines  of  the  angles  (scattering  angles) 
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between  (y,$)  and  (y^,  n + 4>g)  and  between  (y,<f>)  and  (-y^,^),  respectively.  For 
0)q  = 1,  Eq.(22)  must  be  evaluated  in  limit  X 0. 

The  scattering  phase  function  that  we  elected  to  use  with  Eq.  (22)  is 
the  Henyey-Greenstein  function 

P<V  = PH-G(ys)  “ ^ l~g2)  ( l+g2-2gyg )_3/^2  , (23) 

where  g is  an  anisotropy  parameter  defined  as  the  average  value  of  yg,  the  cosine 
of  the  scattering  angle: 
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PH-G<"») 


dy  d<f> 
s 
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The  anisotropy  parameter  n in  Eq.  (22)  is  related  to  g by 


(24) 


n - 


1±£ 

2g 


d+g2)1/2 


(25) 


All  of  the  cloud  spectral  radiance  calculations  presented  in  the 
following  two  sections  were  obtained  using  Eqs.  (22),  (23)  and  (25).  The 
calculated  spectra  always  represent  radiance  fields  above  the  clouds:  hence 
t » 0.  Note  that  the  geometry  of  the  scattering  calculation  (before  we  apply 
corrections  for  atmospheric  molecular  absorption/emission)  is  defined  completely 
by  the  three  parameters  y^,  y and 

As  mentioned  earlier,  the  aqrosol  scattering  parameters  used  for 
calculations  consist  of  the  AFGL-supplied  spectral  data  Wq(V),  g(v)  and  Kq(v)  (the 
extinction  coefficient  per  particle).  The  data  sets  corresponding  to  different 
size  distributions  of  water  droplets  and  ice  spheres  are  given  in  Appendix  A.  In 
performing  the  cloud  spectral  reflection  calculations  we  first  select  the  cloud 
type  (water  or  ice),  the  cloud  optical  thickness  at  2.6  micrometers,  and  the  mean 


radius  of  the  cloud  particles.  Then,  using  the  appropriate  data  set,  we  convert 
g(V)  to  0(V)  via  Eq.  (25),  and  convert  Kq(v)  to  Tq(\>)  using 

T0(V)  “ ^Kq(v)/K0(2.6  pm)J  tq(2.6  W1")  • (26) 

For  the  calculations  we  chose  Tq(2.6  pm)  = 10,  which  insures  that  the  cloud  will 
' optically  thick  over  the  entire  2.6  to  5-micrometer  region  for  any  selected 
cloud  type  or  mean  radius.  Thus,  it  was  not  necessary  to  perform  the  scaling 
defined  by  Eq.  (26),  but  this  operation  was  included  in  the  computer  code  so  that 
it  would  have  a greater  range  of  application  (i.e.,  to  arbitrary  cloud  optical 
thickness  and  P * 0). 

3.4  EXAMPLE  CALCULATIONS 

Figures  3 through  show  some  of  the  spatial  radiance  variations 

predicted  by  the  model.  The  calculations  do  not  include  atmospheric  absorption  or 
emission,  or  cloud  emission.  Figure  3 shows  the  variation  of  l(0,p,4>)  with  0 * 
cos  *p  at  wavelength  3.9  micrometers,  for  0q  * cos  ^p^  a 60  degrees  and  <{>—  <5^  * 0 
or  180  degrees.  The  cloud  consists  of  water  droplets  with  r^  ■ 4 micrometers. 
Note  that  the  radiance  increases  with  but  for  > 4 the  cloud  is  optically 
thick  and  there  is  essentially  no  further  increase.  Note  also  that  the  radiance 
is  a maximum  in  the  backscatter  direction  for  optically  thick  clouds,  but  for  thin 
clouds  the  largest  radiance  occurs  at  0 ■ 90°,  “ 0.  Figure  4 is  the  same 

as  Figure  3 except  * 16  micrometers.  Figure  5 shows  that  the  reflected 
radiance  can  have  a strong  wavelength  dependence  through  0)q  and  g.  Figures  6 
and  7 show  the  azimuth  variation  of  the  radiance  for  three  different  solar 
zenith  angles;  Figure  6 represents  a water  cloud  with  r^  » 16  micrometers  at  the 
wavelength  3.9  micrometers;  Figure  7 represents  an  ice  cloud  with  rffl  ■ 
32  micrometers  at  the  wavelength  4.0  micrometers. 
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3.5 


VALIDATION  OP  THE  CLOUD  RADIANCE  MODEL 


Seventeen  cloud  radiance  spectra  were  selected  from  data  tapes  provided 
by  AFGL,  for  the  purpose  of  validating  the  cloud  radiance  model.  These  data,  from 
the  Background  Measurement  Program,1  were  obtained  by  AFGL  in  support  of  the  IR 
Signature  Studies  program,  and  to  our  knowledge  have  not  been  reported  elsewhere. 
AFGL  provided  all  of  the  data  defining  the  solar  illumination  and  measurement 
geometries,  as  well  as  radiosonde  data  closely  matching  the  time  and  location  of 
measurements.  The  clouds  were  thick  in  all  but  one  of  the  measurements;  for  all  of 
these  cases  the  only  model  parameters  that  were  unspecified  are  the  mean  radius  of 
the  cloud  water  droplets  or  ice  crystals  and  the  cloud  temperature.  Cloud 
temperatures  were  assumed  equal  to  the  radiosonde  temperature  at  the  cloud  top 
altitude . 


Table  1 gives  the  geometrical  data  and  cloud  descriptions  provided  by 
AFGL  with  the  radiance  spectra.  Note  that  in  Mission  727  the  aircraft  observed 
"transmitted"  radiance,  i.e.,  it  viewed  the  brightest  section  of  the  edge  of  a 
cirrus  cloud  from  an  altitude  below  the  cloud.  While  the  cloud  radiance  model  is 
not  intended  to  describe  the  transmitted  radiance  field,  it  is  possible  to  define 
a reflection  geometry  that  represents  the  same  net  scattering  angle  as  the  AFGL 
measurement  geometry  (other  arguments  can  be  given  to  support  the  definition  of  an 
equivalent  reflection  geometry  for  this  particular  case;  see  Section  3.5.2). 
Table  1 also  gives  the  assumed  mean  radius  of  the  cloud  particles, 

Each  of  the  three  different  missions  represented^  in  Table  l cor- 
responds to  a particular  time  and  location,  and  is  associated  with  a radiosonde 
data  set  (pressure,  temperature  and  dew  point  versus  altitude).  Some  of  supplied 
radiosondes  had  no  dew  point  data  for  tne  higher  levels.  In  these  cases  we 
augmented  the  radiosonde  dew  points  using  average  mid-latitude  values  given  in  the 
Air  Force  Handbook  of  Geophysics.^  Appendix  B gives  the  complete  profiles  used  in 
the  calculations,  i.e.,  the  atmospheric  data  used  as  inputs  to  the  modified 
LOWTRAN  4 code.  (Appendix  B also  contains  augmented  profiles  for  the  AFGL 
missions  that  provided  data  relevant  to  the  other  models  developed  in  this  study). 


3.5.1 


Mission  726/7 


Figure  8 shows  four  AFGL  spectra  of  a low  stratus  cloud  dock  observed 
over  the  ocean  nerr  San  Louis  Obispo,  California,  on  28  September  1977,  and  the 
corresponding  calculated  spectra  (the  lower-resolution  curves).  The  logarithmic 
ordinate  scale  applies  to  record  178  of  Mission  726  /7 ; the  other  three  spectra, 
records  176,  175  and  177,  have  been  displaced  downward  by  factors  of  2,  4 and  8, 
respectively.  The  brightest  spectrum  is  record  178,  corresponding  to  a relative 
solar  azimuth  direction  of  5,51  degrees,  while  the  second  brightest  is  record  176, 
corresponding  to  a solar  azimuth  of  175.5  degrees;  i.e.,  "forward"  scattering 
results  in  higher  radiances  then  "backscattering"  for  these  particular  solar 
elevation  and  line-of-sight  depression  angles.  The  N-0  band  at  2563  cm  1 and  the 

—l  L 

HDO  bands  near  2720  cm  are  prominent  features  in  these  and  many  other  of  the 
AFGL  cloud  spectra.  The  4.3-micrometer  CO^  band  is  not  a strong  feature  in  these 
spectra  because  thermal  emission  dominates  reflected  solar  radiation  in  this 
spectral  region,  and  because  there  is  a nearly  isothermal  layer  immediately  below 
the  aircraft  (see  Appendix  B). 

The  calculated  cloud  radiance  spectra  are  shown  separately  (and  on  a 
common  ordinate  scale)  in  Figure  9,  Results  are  shown  for  mean  water  droplet 
radii  of  8 and  16  micrometers.  The  calculated  spectra  in  Figure  8 correspond  to 
a mean  radius  of  16  micrometers.  It  is  evident  that  the  larger  radius  gives 
better  agreement  wi^h  the  AFGL  data. 

■■  ' _1 

The  AFGL  spectral  radiances  below  2100  cm  in  Figure  8 are  probably 

too  large  as  a result  of  known  calibration  inaccuracies.  The  calculated 
’ ~\  _j 

radiances  between  2400  cm  and  2500  cm  are  too  large,  probably  because  of  a 
deficiency  in  the  LOWTRAN  model,  which  seems  always  to  predict  too  high  a 
transmittance  in  this  portion  of  the  wing  of  the  4.3-micrometer  COj  band.  We  have 
found  that  the  LOWTRAN  code  also  tends  consistently  to  overestimate  the 

transmittance  in  the  most  transparent  region  of  the  3. 2-micrometer  H„0  band,  near 

—1  1 — 1 

3150  cm  . Differences  between  the  measured  and  calculated  spectra  near  2560  cm 

clearly  indicate  that  the  atmosphere  contains  significantly  more  N^O  than  is 

assumed  in  the  LOWTRAN  code. 
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Considering  that  Che  differences  noted  above  are  all  explainable,  at 
least  qualitatively,  and  are  attributable  mostly  to  the  atmospheric  transmission 
part  of  the  calculation,  it  appears  that  scattering  model  developed  in  Section  3.3 
gives  essentially  correct  results  for  the  conditions  represented  by  mission  726/7. 

3.5.2  Mission  727/9 

As  mentioned,  the  AFGL  aircraft  observed  "transmitted"  radiance,  while 
the  model  describes  only  "reflected"  radiance.  The  underside  of  the  cirrus  cloud 
was  observed  at  the  same  elevation  as  the  sun  (14.25  deg)  but  at  an  azimuth 
approximately  5 degrees  different  from  the  solar  azimuth.  Thus,  the  net 
scattering  angle  undergone  by  the  observed  radiation  is  roughly  5 degrees.  It  was 
reported  that  the  AFGL  interferometer  observed  the  brightest  part  of  the  underside 
of  the  cloud  (the  so-called  "silver  lining"). 

Figure  10  shows  the  actual  measurement  geometry  (1)  and  the  scatter- 
ing geometry  assumed  in  calculating  the  cloud  spectral  radiance  (2).  That  is,  in 
the  scattering  calculation  (which  determines  the  ratio  of  reflected  to  incident 
spectral  radiance  at  the  top  of  the  cloud)  we  assumed  the  solar  elevation  and 
1 ine-of-sight  depression  angles  were  both  2.5  degrees,  which  corresponds  to  a net 
scattering  angle  of  5 degrees.  However,  in  the  atmospheric  transmittance  and 
radiance  calculations  the  actual  measurement  geometry  was  used.  Because  the  net 
scattering  angles  are  approximately  equal  for  the  two  geometries,  one  might  expect 
the  calculated  spectral  radiance  to  equal  the  measured  spectral  radiance,  provided 
the  average  optical  thickness  traversed  by  received  photons  is  the  same  for  both 
geometries . 

The  lower  part  of  Figure  10  indicates  the  general  trends  to  be 
expected  in  the  transmitted  and  reflected  radiances  as  a function  of  distance  from 
the  cloud  edge  (increasing  optical  thickness).  Normally,  the  transmitted  diffuse 
radiation  will  first  increase  due  to  the  increased  number  of  scatterers,  but  will 
eventually  decrease  to  zero  (foroig  < 1)  due  to  absorption  by  the  cloud  crystals. 
AFGL  measured  the  transmitted  spectral  radiance  at  the  location  that  corresponded 
to  maximum  visual  brightness.  The  reflected  radiance  is  predicted  by  our  model  to 
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increase  monotonically  with  optical  thickness,  and  to  approach  an  asymptote  (which 
is  effectively  reached  at  t a 4).  It  is  not  unreasonable  to  expect  that  the 
average  optical  thickness  traversed  by  the  reflected  radiation,  when  it  has 
reached  its  asymtote,  is  the  same  as  the  optical  thickness  traversed  by  the 
transmitted  radiation  when  it  is  a maximum,  if  the  net  scattering  angles  are  the 
same  in  both  cases.  This  assumption  implies  that  the  maximum  transmitted  and 
asymptotic  reflected  radiances  are  equal.  This,  in  any  case,  was  the  basis  for 
choosing  an  "equivalent"  reflection  geometry  for  the  scattering  calculation.  We 
assumed  an  optically  thick  cloud  (tq  ■ 10)  and  a mean  ice  particle  radius  equal  to 
32  micrometers,  which  is  a reasonable  value  for  cirrus  clouds. 

Figure  11  shows  the  measured  and  calculated  spectral  radiances.  The 
agreement  is  excellent,  but  it  may  be  fortuitous,  since  we  cannot  prove  the 
equivalence  of  the  transmission  and  reflection  geometries,  and  because  the  mean 
radius  was  selected  arbitrarily. 

3.5.3  Mission  817,  Runs  1,  2 and  3 

In  mission  817  the  spectral  radiance  of  a thick,  uniform  stratocirrus 
cloud  deck  was  measured  between  late  afternoon  and  evening,  corresponding  to  solar 
elevations  ranging  from  8 degrees  to— 1.5  degrees.  The  selected  spectra,  runs  1, 
2 and  3 represent  the  range  of  elevations  8 to  2.6  degrees,  and  four  different 
azimuth  directions  for  each  elevation  angle. 

Figure  12  shows  the  measured  and  calculated  radiance  spectra  for 
Run  1.  The  brightest  AFGL  spectrum  is  record  365,  corresponding  to  zero  relative 
azimuth;  the  second  brightest  is  record  367  (A$  * 180  degrees).  The  other  two 
AFGL  spectra,  for  relative  azimuths  of  90  and  270  degrees,  are  nearly  equal,  as 
one  would  expect . 

Agreement  between  the  calculated  and  measured  spectra  is  generally 
poor.  Disagreement  in  the  1800  to  2200  cm"”1  region  for  the  A$  ■ 90  and  270  degree 
cases  suggests  the  cloud  was  warmer  than  assumed.  Disagreement  in  the  2400  to 
3400  cm  * region  for  the  A$  ■ 0 and  A$  ■ 180  cases  and  the  general  trends 
predicted  by  the  model  suggest  that  the  cloud  crystals  may  be  significantly  larger 
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CLOUD  RADIANCE  SPECTRA- 


Figure  12.  Measured  and  Calculated  Spectra  for  AFGL  Mission  817/1 
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than  the  assumed  value,  32  micrometers.  Also,  there  may  have  been  some  structure 
in  the  cloud  top  which  cannot  be  accounted  for  in  the  model:  note  that  at  2700  cm 
the  measured  radiance  is  the  same  for  A$  - 90  degrees  and  180  degrees,  and  is  less 
than  that  measured  for  a4>  ■ 270  degrees.  Generally,  the  agreement  was  worse  when 
smaller  optical  thicknesses  and/or  particle  sizes  were  used.  We  could  not  perform 
calculations  for  larger  cloud  particle  sizes  because  the  largest  radius 
represented  in  the  scattering  parameter  data  set  is  32  micrometers.  Figure  13 
reproduces  the  calculated  spectra  shown  in  Figure  12. 

Figures  14  and  15  show  the  measured  and  computed  spectra  for 
mission  817,  run  2.  Note  that  there  are  four  distinct  spectra  for  this  case 
because  the  relative  azimuths  are  0,  90,  180  and  300  (rather  than  270)  degrees. 
Figures  16  and  17  show  the  results  for  mission  817,  run  3.  As  for  run  1,  the 
agreement  between  measured  and  calculated  spectra  is  poor,  and  suggests  that  the 
assumed  cloud  temperature  and  particle  size  are  both  too  small.  For  run  3, 
however,  the  brightest  calculated  specL 'urn  reproduces  the  spectral  variation  of 
the  corresponding  measured  radiance  quite  well. 
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CLOUD  RADIANCE  SPECTRA-  817/2 


Figure  14.  Measured  and  Calculated  Radiance  Spectra  for  AFGL  Mission  817/2 


CLOUD  RflDIRNCE  SPECTRA- 


Figure  16.  Measured  and  Calculated  Radiance  Spectra  for  AFGL  Mission  317/3 


SECTION  IV 

MODELS  FOR  SUNGLITTER  AND  SPECULAR  REFLECTION  BY  WATER 


4.1  SUNGLITTER  MODEL 

The  reflection  of  sunlight  from  a rough  water  surface  is  visually 
observed  to  consist  of  a broad  intense  glitter  pattern.  This  pattern  is  brightest 
when  looking  along  the  azimuth  toward  the  sun  and  diminishes  symmetrically  on 
either  side  of  this  direction.  At  any  given  instant  in  time  the  glitter  pattern  is 
seen  to  consist  of  many  highlights  These  highlights  are  single  small  reflecting 
elements  or  facets  of  the  surface  which  are  momentarily  oriented  to  specularly 
reflect  the  image  of  the  sun  in  the  direction  of  the  observei . The  number  of 
specularly  reflecting  highlights  in  a fixed  IFOV  is  a random  function  of  the 
surface  wind  direction  and  velocity.  Consequently,  the  apparent  brightness  within 
the  IFOV  results  from  spatial  averaging  over  many  unresolved  highlights  in  the 
scene. 


The  problem  of  calculating  the  average  sung] itter  radiance  has 
attracted  the  attention  of  many  investigators  who  were  motivated  by  the  need  for 
an  analytical  model  to  aid  in  assessing  the  performance  of  shipboard  infrared  and 
radar  search  and  t.rack  systems  operating  against  a sunglint  ocean  background.  The 
radar  system  analysts  have  traditionally  treated  sunglitter  as  a physical  optics 
problem  involving  a vector  form  of  the  Kirchoff  integral.  On  the  other  hand,  the 
infrared  systems  community  have  preferred  the  geometric  optics  approach  which 
consists  of  spatially  averaging  the  contributions  from  many  randomly  oriented, 
unresolved  highlights  about  the  point  of  interest. 

I 

The  validity  of  both  the  geometric  optics  and  the  physical  optics 
solutions  is  the  subject  of  current  research.  The  geometric  optics  approach  of 
Cox  and  Munk^  is  employed  in  this  study  because  of  its  simple  and  physically 
intuitive  interpretation  of  the  sunglitter  phenomenon. 
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The  basic  expression  for  Che  average  reflected  spectral  radiance  for 
surface  wind  speed  W,  local  angle  of  incidence  to,  and  reflection  geometry  angles 
<*.8.V,y>  is: 

<Nsun(X)>  " [p(w>X)  H0(X)  T(X)]  V,<r)P(Zu,Zv,W)/n  (27) 

The  angular  behavior  of  the  sunglitter  is  given  by  the  products  of  the  surface 
slope  shadowing  factor  (S),  a geometric  factor  (g),  and  a probability  factor  (P) 
which  expresses  the  probability  of  occurrence  of  wave  slopes  having  components  Z 

i U 

and  Z . The  spectral  distribution  of  the  radiance  is  determined  by  the  product  of 
the  atmospheric  spectral  transmission  (T)  along  the  sun-surface-observer  path, 
the  solar  spectral  irradiance  (Hq)  at  the  top  of  the  atmosphere,  and  the  spectral 
reflectivity  (P)  of  the  ocean  surface. 

It  is  assumed  that  warning  sensors  are  insensitive  to  the  polarization 
properties  of  the  surface,  so  that  p(X)  is  the  unpolarized  surface  reflectivity: 

2p«D,X)  - pv  (to,X)  + ph  (to,X)  (28) 

where  py  and  are  the  Fresnel  spectral  reflection  coefficients  for  vertical  and 
horizontal  polarization,  respectively. 

Each  of  the  factors  in  Eq.  (27)  is  a complicated  function  of  source- 

receiver  geometry,  wind  speed  and  direction,  surface  reflectivity,  and 

wavelength.  Explicit  expressions  for  each  of  these  factors  - p,S,g,  and  P — are 

given  in  Appendix  C.  These  expressions  are  sufficiently  complicated  to  require  a 

computer  code  to  calculate  <^N  C X ) *>. 

sun 

A search  of  the  available  literature  revealed  that  the  only  available 
validations  of  the  Cox-Munk  model  involved  broad-band  IR  measurements  of  the 
angular  variation  of  ^ with  wind  speed  and  source-receiver  geometry. 

Furthermore,  all  of  the  previous  validations  included  other  components  of  the  sea 
background,  such  as  reflected  skylight,  upweiling  ocean  radiance,  and  thermal  sea 
surface  emission. 
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The  strong  liquid  water  absorption  in  the  infrared  allows  us  to  safely 
ignore  the  upwelling  ocean  radiation  caused  by  volume  scattering  of  the  incident 
solar  radiation  from  beneath  the  ocean  surface. 


Thermal  emission  is  a significant  component  of  the  ocean  background 

8 

radiance  and  can  be  evaluated  using  the  Cox-Munk  technique.  Krishnan  and  Peppers 
have  found  that  average  thermal  emission  for  strong  ocean  absorption  is  given  by: 

<N(X)>th  “ > f(0,  Wu,  Wv)]n0(X,T)  (29) 

where  Nq(X,T)  is  the  blackbody  radiance  for  sea  temperature  T,  and  f(0,  Wu»  Wy) 
represents  an  average  surface  reflectance: 


:(B,  Wu,  Wy)  * cos  B P (8)  ]l+ 


P"(B)  - 2p*(0)  tan  (P) 

\ 

fa  2 cos2  a + a ^ sin2  a] 

Lu  v J 

+ 2p"(g~)  ^p'(B)  ctn  <0)j*£cru2  Sin2  a + ay2  cos2  aj 


(30) 


The  quantities  p'  and  p"  are  the  first  and  second  derivations  of  the  reflection 
coefficient  with  respect  to  surface  slopes,  o and  p are  the  elevation  aud  azimuth 
of  the  observed  radiance,  and  a , a are  the  rras  variation  in  surface  slopes  in 
directions  u and  v.  These  quantities  are  all  defined  in  greater  detail  in 
Appendix  C. 


The  other  significant  component  of  the  ocean  background  spectral 
radiance  is  reflected  skylight.  Reflected  skylight  can  be  calculated  from  the 
expression  used  for  sunglitter  (Eq.  27)  if  one  replaces  the  solar  irradiance 
term  Hq(X)  T(X)  by  the  irradiance  at  the  sea  surface  due  to  a small  segment  of  the 
sky  and  then  integrates  the  resulting  expression  over  the  sky  dome.  Krishnan  and 
Peppers  assumed  the  sky  spectral  radiance  Nfl(X)  was  uniform  over  the  entire  sky 
dome  and  derived  an  expression  for  the  average  reflected  skylight  spectral 
radiance  in  terms  of  Ng(X)  and  the  average  ocean  reflectance  f(0,  Wu,  Wy): 
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<N(X)>sky  “ Ns(X)  f(8>  Wu>  Wv)  (31) 

Thermal  emission  and  reflected  skylight  backgrounds  are  always  present 
in  sunglitter  measurements,  adding  an  additional  complication  to  the  task  of 
verifying  a sunglitter  model. 

Fortunately  the  AFGL  measurements*  of  5 December,  1974  are  available  in 
a form  which  can  be  used  to  obtain  only  the  sea  glitter  component.  This  data  was 
obtained  from  an  altitude  of  15,000  feet  over  the  Atlantic  Ocean  off  the  coast  of 
New  Hampshire  (pertinent  mission  parameters  are  given  in  Table  2.  Two  sets  of 
measurements  were  made  while  the  aircraft  was  banked  at  an  angle  of  45  degrees  to 
the  horizon.  The  first  measurement  (Figure  18)  was  made  looking  directly  into 
the  glitter  patch.  The  second  measurement  (Figure  19)  was  made  at  an  azimuth  of 
30  degrees  away  from  the  visible  sunglint  path.  A computer  calculation  for  the 
second  case  indicated  that  the  sunglint  should  be  reduced  by  a factor  of  1000  when 
compared  to  the  first  measurement.  Consequently,  the  second  spectrum  contains 
only  background  and  atmospheric  emissions  which  may  be  removed  from  the  first 
spectrum. 


TABLE  2 

MISSION  PARAMETERS  FOR  AFGL  FLIGHT  TR3 , RUN  1 


Date: 

5 December  1974 

Universal  Time: 

1639 

Latitude  (deg): 

43.26 

Longitude  (deg): 

70.08 

Wind  Direction  (deg): 

320-Concord,  N.H. 
260-Pease  AFB 

Wind  Speed  (knots 

):  5.0 

Visibility  (nmi): 

30.0 

Radiosonde  Data: 

1200  ITT 

Portland, 

Maine 

Figure  20  shows  the  difference  between  Figures  14  and  A very 
large  sunglint  component  appears  in  the  data  between  3000  and  6000  cm  . The  deep 
atmospheric  absorption  bands  from  5100  to  5650  cm  * and  from  3000  to  4000  cm  * are 
caused  by  the  1.9  and  2 . 7-micrometer  water  vapor  absorption  bands.  Two  narrow  CO,^ 
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absorption  bands  are  also  present  between  4800  and  5000  cm  AFGL  suspects  the 
validity  of  their  data  in  the  blank  area  from  2800  to  3500  cm  * in  Figure  19. 

The  sunglint  spectrum,  unlike  other  natural  backgrounds  is  decidedly 
non-Planckian  because  of  the  abrupt  change  in  the  ocean  surface  reflectivity  at 
2.94  micrometer  This  feature  is  evident  in  the  curve  labeled  CM  in  Figure  20. 
The  CM  curve  is  the  computed  result  using  the  Cox-Munk  expression  for  sunglint 
with  100%  atmospheric  transmission  for  the  geometry  and  wind  conditions  given  in 
Table  2.  This  curve  also  ignores  shadowing  (S)  which  only  becomes  significant 
at  grazing  angles. 

The  values  displayed  as  crosses  in  Figure  20  represent  the  reflected 
spectral  radiance  when  atmospheric  attenuation  is  included.  The  atmospheric 
spectral  transmittance  of  the  sun-surface-aircraft  path  was  computed  using  the 
modified  LOWTRAN  4 code  described  in  Section  III. 

In  general  the  model  predictions  agree  well  with  the  AFGL  measurements. 
The  differences  between  the  AFGL  data  and  the  transmittance-corrected  theoretical 
spectrum  are  probably  due  to  an  inaccurate  estimate  of  the  aerosol  optical 
thickness  along  the  transmission  path. 

4.2  SPECULAR  REFLECTION  MODEL 


simply 


The  spectral  radiance  reflected  specularly  by  a smooth  water  surface  is 


h(W,X)  = p((i),X)T(X)N0U),  (32) 

where  p(to,X)  is  the  surface  reflectivity  for  incidence  angle  uj,  Nq(X)  is  the 
extraterrestrial  spectral  radiance  of  the  sun,  and  T(X)  is  the  spectral  trans- 
mittance of  the  path  from  the  sun  to  the  water  surface  and  then  to  the  receiver. 
The  reflectivity  is  given  by  Eq . (28)  (see  Appendix  C also). 

AFGL  spectra  of  specular  reflections  were  compared  to  the  predictions  of 
Eq.  (32).  In  calculating  (li),X)  we  used  the  water  indices  of  refraction  tabulated 
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3 

by  Hale  and  Querry  . These  data  apply  to  pure  water,  but  the  solutes  typical  in 

ocear  water  will  not  appreciably  change  the  refractive  indices.  However,  the 

presence  of  films  such  as  oil  slicks  or  other  pollutants  can  have  a large  effect  on 

0(to,X).  One  of  the  objectives  in  analysing  the  AFGL  spectra  was  to  determine 

whether  they  differ  significantly  from  predictions  based  on  the  assumption  of  pure 

water.  The  Nn(X)  values  used  were  ba^cc  on  the  solar  irradiance  spectrum  of 
9 u -5 

Thekaekara  and  the  value  1.9ft  x 10  sr  for  the  average  solid  angle  subtended  by 
the  sun. 

Validation  of  the  specular  sunglint  model  for  smooth  water  surfaces  is 
complicated  by  other  backgrounds  in  the  receiver  field  of  view.  Thermal  emission 
and  reflected  skylight  are  always  present  in  measurements.  Fortunately,  the  AFGL 
flight  data  consists  of  2 sets  of  spectra.  One  set  contains  glint  plus  background 
and  the  other  set  consists  of  background  only.  The  difference  between  the  two 
sets  of  data  allows  one  to  isolate  the  glint  data  for  model  verification. 

Five  spectra  of  specular  sunglints  and  associated  backgrounds  were 
selected  from  AFGL  missions  729  and  803.  These  are  identified  in  Table  3.  The 
atmospheric  profiles  used  in  the  modified  LOWTRAN  4 calculations  were  based  on 
AFGL-supplied  radiosondes,  and  arc  given  in  Appendix  B.  The  differential  radiance 
spectra  (glint-plus-background  minus  background  only)  were  calculated  and  are 


displayed  as  curves  labeled  "AFGL"  in  Figures  21  through  25.  Included  in  each  of 

these  spectra  is  a multiplication  factor  of  16  which  converts  "apparent" 
spectral  radiance  to  true  spectral  radiance  within  the  solar  image;  i.e.,  the 
instrument  IFOV  diameter  was  two  degrees  and  the  sun's  image  subtends  approxi- 
mately 0.5  degrees. 


Figures  21  through  25  also  show  the  predicted  spectral  radiances  with 
and  without  atmospheric  attenuation  included  (the  cross  symbols  and  the  low 
resolution  solid  curves,  respectively).  Each  pair  of  predicted  spectra  have  been 
multiplied  by  a scaling  factor  a that  forces  agreement  between  the  predicted 
values  (those  with  atmospheric  attenuation  included)  and  the  measured  line 
radiance  values  in  the  window  regions  2400-3000  cm  * and  4000--4800  cm  ^ . These 
scale  factors,  given  in  Table  4,  may  be  regarded  as  the  fraction  of  the  solar  disc 
observed;  i.e.,  in  some  cases  a broken  image  may  have  been  observed  due  to  surface 
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roughness  (wind),  and  in  ac  least  one  case  (303/15)  the  image  was  not  whole 
because  mountains  blocked  some  of  the  low-angle  incident  sunlight.  Scale  factors 
were  determined  separately  for  each  of  the  two  window  regions  (see  Table  4)}  a 
weighted  average  was  used  for  the  plots  shown  in  Figures  21  through  25. 

There  appears  to  be  excellent  spectral  agreement  between  theory  and 
observation  in  all  of  the  plots.  The  slight  spectral  discrepancies  are  probably 
caused  by  the  inability  of  LOWTRAN  4 to  accurately  model  the  atmospheric  trans- 
mission over  the  long  slant  ranges  corresponding  to  the  low  sun  angles.  (There  is 
at  least  2.27  air  masses  between  the  (sun  and  the  water  surface  for  the  highest 
solar  elevation  — 26  degrees  and  at  least  8.90  air  masses  for  the  lowest  solar 
elevation  - 6 degrees.  The  only  suggestion  of  a spectral  reflectivity  variation 
different  from  that  of  pure  water  is  in  the  mission  803/4  case  for  Mono  Lake. 
However  the  predicted  spectrum  for  the  same  lake  observed  at  lower  sun  elevation 
(803/9)  shows  excellent  spectral  agreement  with  the  AFGL  data,  and  has  a scale 
factor  of  0.93. 

The  only  disturbing  aspect  of  the  results  is  the  case-to-case  variation 
in  the  predicted  radiance  levels.  The  best  theoretical  agreement  with  data  occurs 
for  the  Mono  Lake  spectrogram  at  a low  sun  angle  where  a ■ 0.93.  Yet  the  same 
scene  recorded  at  a higher  solar  elevation  has  0 * 0.02.  In  other  words,  a nearly 
a fifty-fold  reduction  in  scaling  factor  is  needed  to  explain  the  AFGL  data  for 
the  same  scene  observed  at  slightly  different  sun  angles. 
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SECTION  V 

BLUE  SPIKE  SPECTRAL  EMISSION/ ABSORPTION  MODEL 


An  efficient  model  and  computer  code  were  developed  to  calculate  the 
fully  resolved  spectrum  of  the  blue  spike  of  the  4.3-micrometer  CO^  band,  a well- 
understood  phenomenon  resulting  from  the  emission  of  hot  CO^  transmitted  through  a 
cooler  layer  (atmospheric  path)  also  containing  CO^.  The  model  avoids  using  the 
lengthy  line-by-line  procedure  by  representing  the  hot-and  cold-gas  spectral 
absorption  coefficients  as  two  different  Elsasser  bands  with  spectrally  varying 
line  intensity  S and  spacing  d.  Theoretical  formulas  for  the  intensities, 
halfwidchs  and  spacings  of  the  rotational  lines  of  the  00°1— 00°0  transition  are 
used  with  empirical  modifications  that  force  agreement  with  measured  spectra  and 
computed  line-by-line  spectra.  The  model  places  the  Elsasser  spectral  absorption 
maxima  at  the  exact  theoretical  line  positions,  which  preserves  the  observed 
spectrum  features  and,  in  particular,  the  correlation  of  emission  and  absorption 
lines.  The  model  sacrifices  very  little  accuracy  and  is  potentially  30  to  60 
times  faster  than  a conventional  line-by-line  calculation  of  the  blue  spike 
spectrum. 


Figure  26  shows  a measured  absorptance  spectrum  of  the  4 . 3-micrometer 
blue  spike  produced  by  a flame  at  2330  K.  A notable  characteristic  of  the 
spectrum  is  its  resemblance  to  the  periodic  Elsasser^  absorption  coefficient 

k(V)  * d cosh'  £8-"cos  ITmT  6 ' «3> 

which  represents  exactly  an  infinite  periodic  succession  of  equal  Lorentz  lines  of 
intensity  S,  halfwidth  Y and  spacing  d.  In  this  equation  v represents  the 
frequency  measured  from  the  center  of  any  line.  Of  course,  the  intensities  and 
apacings  are  not  constant  in  the  actual  (measured)  spectrum,  but  one  might  assume 
that  its  absorption  coefficient  can  be  approximated  by  Eq.  (33)  if  S and  d are 
allowed  to  be  slowly  varying  functions  of  frequency  v. 
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In  fact,  the  Elaasser  model  cannot  represent  the  actual  spectrum,  even 
when  variable  S and  d are  used,  because  the  model  assumes  Lorents  lines  whereas 
measurements  show  that  the  true  CO.  line  shape  iB  "sub-Lorentrian,"  i.e.,  is 

1 t ii  12 

described  by  the  well-known  Benedict  line  shape  or  those  propounded  by  Burch 
or  Yamamoto.  However,  it  is  still  possible  to  use  a modified  form  of  the 
Elaasser  model  to  mimic  the  essential  features  of  a measured  spectrum,  namely  the 
heights  of  the  absorption  maxima  and  the  levels  of  the  minima  between  absorption 
lines.  The  only  differences  then  remaining  are  the  values  of  absorption  coeffi- 
cient between  maxima  and  minima,  which  for  the  purposes  of  this  study  can  surely 
be  considered  a second-order  effect. 

The  heights  of  the  maxima  and  minima  in  the  Elaasser  absorption  coef- 
ficient are  obtained  by  setting  the  denominator  term  cos(2uv/d)  equal  to  i 1 . We 
find  that: 


k 


max 


S 

d 


k . ■ tanh  (0/2). 

mm  d 


(34) 


Hence,  the  product  of  maximum  and  minimum  absorption  coefficients  is: 


k 


k . 


max  min 


(S/d)2, 


(35) 


which  is  equivalent  to  the  equation 

log  (S/d)  - | [dog  kmax  ♦ log  ka.n)^j  (36) 

This  shows  that  the  median  ordinate  of  the  upper  and  lower  envelopes  of  the 
absorption  coefficient  plotted  on  a logarithmic  scale  is  equal  to  log(S/d).  We 
can  also  show  that: 


0 • 2 tanh 


-1 
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1/2 


(lo8I0kmin 
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1 

>J 


(37) 


In  other  words,  the  average  of  the  logarithmic  absorption  coefficient  and  the 
total  logarithmic  excursion  uniquely  determine  S/d  and  0,  the  two  parameters  of 
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Che  Elsasser  model.  It  is  easy  to  see  how  S(v)  and  Y(v)  [which  determines  3(b)  ■ 
2irY(v)/d(v)]  can  be  empirically  adjusted  so  that  the  Elsasser  model  fits  the 
maxima  and  minima  of  a measured  blue  spike  spectrum,  and  also  retains  the  actual 
line  positions  and  line  spacing  d(v). 

Our  calculations  represent  the  00°1— 00°0  transition,  which  is  the  only 

one  of  importance  in  the  blue  spike  region.  The  rotational  lines  included  have  J 

values  from  50  to  122  (corresponding  to  the  bandhead),  and  from  122  to  172  in  the 

return  line  spectrum.  Beyond  the  bandhead  we  represent  the  line  component  of  the 

spectrum  as  one  side  of  a Benedict  spectral  line  centered  at  the  bandhead  and 

matching  the  absorption  coefficient  at  the  center  of  the  J ■ 122  line.  However, 

in  the  region  near  and  beyond  the  bandhead,  in  the  case  of  the  cool  gas 

(atmospheric)  absorption  coefficient,  we  add  a continuum  contribution  represented 

by  the  sum  of  the  wings  of  two  Benedict  lines  placed  near  the  "centers  of  gravity" 

of  the  P and  R branches  of  the  00°l-00°0  band;  this  follows  a procedure  defined  by 
14 

Carpenter  to  represent  the  band  continuum.  This  continuum  correction  is 
unnecessary  for  the  hot  gas  absorption  coefficient  at  gas  temperatures  above 
500  K. 


After  computing  he  hot-gas  and  cool-gas  spectral  absorption  coef- 
ficients in  this  manner,  the  code  computes  the  spectral  radiant  emission  of  a hot- 
gas  slab  or  cylinder  of  specified  properties,  and  the  spectral  transmittance  of  a 
specified  cool-gas  (atmospheric)  path.  Finally,  it  multiplies  these  two  spectra 
together  on  a point-by-point  basis  to  obtain  the  received  blue-spike  spectral 
radiance. 

5.1  DETAILS  OF  THE  MODEL 

We  will  first  define  the  basic  equations  of  the  model,  then  the  empir- 
ical adjustments.  The  00°1-00°0  line  positions  are  obtained  using  the  molecular 
constants  in  the  AFGL  Line  Parameters  Compilation.  The  line  corresponding  to 
ordinal  number  ra  = J + 1 is  centered  at  wavenumber 
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V (cm"1)  - 2349.146  + mio.777358-m[o. 003078  + 


ra(5 . 312  x 10"7  - 6 x 10-1°m)]j, 


m - 51(2)173 

Note  that  m assumes  only  odd  values,  corresponding  to  even  J. 

The  approximate  line  spacing  d is  obtained  by  differentiating  Eq. 
(38)  with  respect  to  ra  and  then  multiplying  by  2: 

dm(cra_1)  - 2^0 .777358  - m[0. 006156  + 

ra(l. 5936  x 10~6  - 2.4  x 10~9  m)]|,  (39) 

ra  < 121 

This  equation  is  not  sufficiently  accurate  at  the  band  head  (m  ■ 123);  Eq. 
(38)  is  used  to  obtain  the  distance  between  the  lines  ra  ■ 121  and  ra  ■ 123. 


The  R-branch  return  lines  corresponding  to  ra  >123  are  assumed  to  fall 
precisely  on  the  outgoing  lines;  i.e.,  the  lines  ra  ■ 125,  127,  ...  are  assumed  to 
fall  on  top  of  the  lines  ra  * 121,  119,  ...,  respectively,  so  that  the  corre- 
sponding line  intensities  can  be  added.  This  is  actually  a very  good  approxi- 
mation for  the  first  few  return  lines,  and  beyend  these  the  line  intensities  have 
decreased  to  such  a low  level  that  their  contributions  are  insignificant. 


For  the  individual  line  intensities  we  used  the  well-known  approximate 


formula 


Sm  - So(300  K)  m exp(-l  .439m(ra-l)  •—)  (40) 
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together  with  the  values  (McClatchey 


15 


and  Gray^) 


S (300  K)  - 2443.0  cm 
o 


- 0.39022 


The  temperature  dependence  of  the  vibrational  partition  function  Q(T)  is  based  on 

14 

values  given  in  Table  5,  which  were  obtained  from  Figure  5 of  Carpenter.  The 
code  uses  table  lookup  and  interpolation  to  obtain  Q for  a given  temperature. 

TABLE  5 

APPROXIMATE  VIBRATIONAL  PARTITION  FUNCTION  FOR  C02  UP  TO  1500  K 


T (K) 

Q(T) 

100 

1.00 

200 

1.02 

300 

1.15 

400 

1.27 

500 

1.45 

600 

1.75 

700 

’ .9 

800 

- . 2 

900 

2.6 

1000 

3.0 

1200 

4.3 

1500 

6.7 

The  Lorentz  line  halfwidths  Y are  calculated  from: 


Y(  cm  1)  - 0.07(p  + 0.3  prn  ) 

co2 

where  p is  the  gas  total  pressure  in  atmospheres,  T 
and  p . is  the  partial  pressure  of  carbon  dioxide. 


(41) 


is  its  absolute  temperature, 


The  absorption  coefficient  in  the  vicinity  of  the  m^  line  is  assumed  to 


vary  as: 


■■“’■ft) 

e’  - (2TrVd)' 
M m m 


sinh  0' 


m 

cosh  0'o  - cos~£2TT(v-Vn)/d,m} 


(42) 
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where  the  primed  quantities  denote  linear  interpolates  of  S/d,  £ and  d between  the 
line  and  its  nearest  neighbors.  The  absorption  coefficient  is  calculated  by 


the  code  at  20  equal-spaced  points  over  the  interval  d centered  at  v * v . 

m m 


Near 


the  bandhead  is  taken  as  the  sum  of  the  intensities  of  the  outgoing  and  return 


spectrum  lines.  Beyond  the  bandhead  the  absorption  coefficient  is  expressed  as: 


k(v) 


S1 23  ' 


d123  tanh  (8^/2) 


fi-0. 135  <V-Vm> 


0.7 


(v-v123)2 . 


(43) 


which  represents  a Benedict  line  whose  line  center  absorption  coefficient  matches 


that  of  Eq.  (42)  at  V ■ V 


123 


the  bandhead. 


The  continuum  (wing  effect)  for  the  cool  gas  path  is  defined  by  the 


14 

absorption  coefficient  (Carpenter  ) : 


S y 

_E_ 

TT 


<v-v 


(v-v,)- 


-0.135  (V~V  ) 
e o 


0.7 


(44) 


where  v 


2349.146  cm 


-1 


the  00°1-00°0  band  center,  and  S , V_  and  \L  are 

p R P 


intensities  and  wavenumbers  that  depend  on  temperature.  For  T ■ 288  K, 


Sp  a 6.61  x 10  cm/gm 
■ 1.30  x 103  cm"2. 


2362.5  cm 


-1 


(45) 


2331.5  cm 


-1 


14s 


(see  Figure  11  of  Carpenter  ).  The  absorption  coefficient  given  by  Eq.  (44) 
is  added  to  the  net  "line  absorption  coefficient"  given  by  the  previous  equations. 


co2) 


-1 


The  absorption  coefficients  k,  as  defined,  have  the  units  (cm  STP  of 


They  may  be  multiplied  by  273  p/T  times  the  C02  mole  fraction  (cm  STP/kin) 


62 


i ~ - ' . i:-j  Jkifl.*' ir' ti . l. 


! to 


i 

& 


{ 


to  obtain  a corresponding  coefficient  K in  units  of  km  * (per  km  of  physical  path 
length),  if  p is  the  path  total  pressure  in  atmospheres  and  T is  its  absolute 
temperature. 


After  the  spectral  transmittances  of  the  hot-gas  slab  and  cool-gas  p-th 
( t^E  and  t ^)  are  calculated,  the  observed  spectral  radiance  is  computed  from: 


V BV^TB^tVEtVA+  S)A  ^1-tVE^BV^TE^ 


(46) 


+ “"WV* 

where  t^E  ■ exp(— K^Lg)  and  KgLg  is  the  product  of  absorption  coefficient  and  path 
length  (optical  thickness)  for  the  exhaust;  similarly,  t^  ■ exp(— K^L^) , BV(T) 
denotes  Planck's  function,  and  T , T , T denote  the  temperatures  of  the  hot  gas, 
the  cooler  path,  and  a background  assumed  to  be  behind  the  hot  slab  and  to  have 
unit  emissivity,  Note  that  the  cool  path  emission  is  included  in  the  computed 
radiance. 


Eq.  (46)  gives  the  observed  spectral  radiance  of  a hot  slab  source 
! that  fills  the  IFOV.  A more  realistic  geometry  for  a gaseous  threat  or  false 

threat  source  is  a cylinder  of  arbitrary  radius  and  length.  Then,  of  course,  the 
| appropriate  quantity  to  calculate  is  the  spectral  intensity  (watts/sr-cm  *) 

f received  at  some  distant  point. 

j Consider  a cylinder  of  radius  R and  length  SL  which  has  been  sliced  into 

\ 

\ a series  of  slabs  of  thickness  dz  (Figure  27).  The  spectral  radiance  emitted  in 

the  x-direction  by  a typical  slab  is 


[ 1 " exp(-2xKE)jBv(TE), 

and  the  total  emergent  intensity  of  the  entire  cylinder  i3 


I =•  2£B  (T  ) 

V Ci 


f [,  - S™'*] 


dz 


(47) 


63 


2WBv(Te) 


. -f 

Jo  °-t  > 


172 


-2RK_t 
e E dt 


(48) 


2 2 1/2 

where  the  second  form  is  obtained  by  substituting  z ■ (R  — x ) and  t • x/R  in 
the  first  form.  The  integral  in  (48)  can  be  solved  (cf.  Gradshteyn  and 
Ryzhik17,  page  316)  in  terms  of  known  functions  to  obtain 


I - RtByCTg)  [i,(2RKe)  - L1(2RRe)], 


(49) 


RADIANCE  PERPENDICULAR 
TO  CYLINDER  AXIS 


Figure  27.  Geometry  for  Calculating  the  Radiant  Intensity  of 
a Hot  Cylinder  of  Gas 

18 

The  Struve  function  I^(x)  — L^(x)  is  tabulated  in  Abramowitz  and  Stegun  . The 
spectral  intensity  observed  perpendicular  to  the  cylinder  axis  at  distance  is 
equal  to  the  product  of  (5.17)  and  the  cool-path  spectral  transmittance  tvA  ■ 
exp(-KALA) . 
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5.2  EMPIRICAL  ADJUSTMENTS 

The  model  described  above  vas  written  as  a computer  code  that  includes 
the  option  of  selecting  either  the  slab  or  cylinder  geometry  for  the  hot  source. 
Note  that  the  effects  of  atmospheric  and  background  emission  are  included  in  the 
apparent  slab  radiance  calculation  but  not  in  the  apparent  cylinder  intensity 
calculation.  Spectra  obtained  by  the  code  were  compared  to  measurements  and 
independent  line-by-line  calculations  to  determine  the  kinds  of  empirical  adjust- 
ments that  would  improve  the  model  predictions.  The  line-by-line  spectra  used  in 

19 

the  comparisons  were  supplied  by  Horace  Ory  of  R and  D Associates  ; these 
previously  unpublished  spectra  represent  the  cylindrical  source  geometry.  The 
comparisons  resulted  in  the  following  modifications  to  the  equations  of  Sec- 
tion 5.1. 

5.2.1  Variation  of  Line  Width  with  Temperature  and  Rotational  Quantum  Number 

The  Lorentz  halfwidth  at  STP  (the  factor  0.07  in  Eq . (41)  was  made 

a linear  function  of  ordinal  number  m “ J+l  by  fitting  a straight  line  to 

20 

theoretical  results  of  Yamamoto  et  al.  We  also  included  a halfwidth-temperature 
relationship  derived  by  Ory  from  Yamamoto’s  results  (which  he  used  in  the  line-by- 
line calculations).  The  result  of  these  corrections  is  to  c.iange  Equation  5.9  to 

Y ( cm~ 1 ) « (0.077  - 2.7  x 10"4  m)  (p  + 0.3  pPA  ) (300/T)n 

co2 

n * 0. 75  - 0.004(m-l)  for  m < 101  (50) 

a 0 . 35  for  ra  > 101 

This  halfwidth  appears  in  Eqs.  (42)  through  (44). 

5.2.2  Depth  of  Absorption  Minima 

The  absorption  coefficient  minima  between  line  pairs  will  be  too  high 
without  some  kind  of  correction,  because  the  Elsasser  band  represents  Lorentz 


lines.  It  was  found  that  the  trough  depth  could  be  matched  approximately  to  the 
line-by-line  (Ory's)  spectra  by  multiplying  the  in  Eq.  (42)  by  the  factor 

EMPIR  - ^0.726  + 0.00417(m-51)J  1/2  (51) 

A compensating  correction  is  applied  to  the  line  intensities  S Eq.  (40)  so 

m 

that  the  absorption  coefficient  at  the  line  centers  will  not  change;  i.e.,  S is 

m 

multiplied  by 

tanh(B  /2)/tanh|e  /(2EMPIR)]  (52) 

m (_  m J 

where  is  the  new  (empirically  corrected)  value  of  8 (8ee  Eq.  (34). 

5.2.3  End  Correction  to  Elsasser  Band 


This  correction  is  equivalent  to  replacing  an  infinite  array  of  lines 
(the  Elsasser  band)  by  a semi-infinite  array  of  lines  that  ends  at  the  band  head. 
It  consists  of  first  computing  the  quantity 


. S sinh  8 

_ra m 

2 d cosh  8 -1 

m m 


123-m 


1 

( 2nj)2+0  2 ’ 
J m 


m = 85(2)  123, 


(53) 


and  evaluating  the  linear  interpolate  of  (53)  at  the  same  frequencies  used  in 

evaluating  Eq.  (42).  These  values  are  then  subtracted  from  the  absorption 

coefficients  determined  from  Eq.  (42).  Note  that  the  correction  affects  only 

the  spectral  region  spanned  by  the  20  lines  nearest  the  bandhead;  it  is  performed 

after  the  corrections  on  B and  S . 

m m 
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5.2.4  Modification  of  Line  Intensities  and  Cool-Gas  Continuum 

The  intensities  of  the  three  lines  nearest  the  bendhead  (n>  ■ 119,  121 

and  123)  are  multiplied  by  0.94,  0.70  and  0.30,  respectively.  The  cool-path 

continuum  coefficient  given  by  Ec.  (44)  is  increased  uniformly  by  the  factor 

1.74.  During  validations  of  the  model  the  continuum  correction  factor  was  revised 

a second  time  (see  Section  5.3).  Note  that  this  factor  and  the  overall  shape  of 

the  continuum  correction  are  strongly  determined  by  the  extreme  wings  of  strong 

CO  lines;  the  far-wing  line  shape  and  its  temperature  dependence  are  still 
2 21 
uncertain  (see  Burch  and  Gryvnak  ). 

5.3  VALIDATION  OF  THE  BLUE  SPIKE  MODEL 

Table  6 lists  some  of  the  cases  for  which  predicted  spectral  intensity 
or  absorptance  spectra  were  compared  to  Ory's  line-by-line  calculations  and 
measured  spectra.  The  effective  broadening  pressures  given  in  the  table  are 
defined  as  the  sum  of  the  total  pressure  and  0.3  times  the  CO^  partial  pressure 
(cf.  Eqs.  (41)  and  (50)). 

The  results  of  the  validations  are  presented  as  a series  of  figures 

which  each  have  a part  (a)  and  part  (b).  Part  (a)  always  shows  a spectral 

quantity  computed  using  the  blue  spike  model,  while  part  (b)  shows  the  same 

19 

quantity  determined  from  measurements  or  Ory's  line-by-line  calculations. 
Table  6 identifies  the  figure  applicable  to  each  case. 

Figures  28  through  37  show  excellent  agreement  betv  n the  model 
predictions  and  the  independently  coraputed/measured  spectra.  The  je<.sured  spec- 
trum of  Figure  34b  contains  emission  lines  of  the  isotope  C 0 ly,  which  has  a 

— 1 ^ 

bandhead  near  2385  cm  . These  lines  are  not  included  in  the  blue  spike  model 

because  the  atmospheric  paths  corresponding  to  typical  AAM  and  SAM  detection 

geometries  will  usually  be  opaque  at  wavenumbers  less  than  2385  cm  Figure  5—10 

shows  that  our  model  predicts  substantially  lower  radiance  than  SAI's  line-by-line 

code  when  an  optically  thick  slab  is  observed  at  long  ranges  (10  km  in  this  case). 

This  is  due  entirely  to  the  atmospheric  continuum  (wing  effect)  used  in  our  model, 
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Figure  28.  Validation  of  Blue  Spike  Model  for  500  K Source 
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SPECTRAL  RADIANCE 


SPECTRAL  RADIA 


Figure  36.  Same  as  Figure  35a,  Except  Based  on  a ised  Correction  | 

Factor  for  the  Atmospheric  Continuum  j 
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Figure  37. 


Same  as  Figure  32a,  Except  Based  on  a Revised  Correction 
Factor  for  the  Atmospheric  Continuum 


for  if  we  reduce  the  continuum  coefficient  by  a factor  of  3 (after  the  empirical 
1.74  increase  referred  to  in  Section  5.2.4)  we  obtain  the  spectrum  shown  in 
Figure  36,  which  is  in  much  better  agreement  with  Figure  35b.  Figure  37 
shows  how  the  predicted  spectrum  of  Figure  32a  would  be  affected  by  the  same 
factor-of-3  change  in  the  atmospheric  continuum. 

The  blue  spike  model  was  also  validated  against  ERIM  measurements  of  the 
spectral  absorptance  of  hot,  undiluted  CO^  at  a pressure  of  20  Torr,  and  at 
temperatures  of  500  R,  700  K and  900  K.  These  results,  reported  in  the  11th  R&D 
Status  Report  of  the  study,  show  excellent  agreement  between  the  calculations  and 
measurements . 

5.4  UTILITY  OF  THE  MODEL 

The  computer  code  written  to  perform  the  blue  spike  spectral  calcula- 
tions has  the  capability  to  represent  both  emission  and  absorption  paths  by  up  to 
20  segments  of  different  pressure,  temperature  and  C0^  fraction.  It  can  also 
degrade  spectral  resolution  according  to  a triangular  slit  function  of  specified 
width.  The  blue  spike  model/code  should  be  useful  for: 

1.  Calculating  the  blue  spike  spectral  signatures  of 
generic  threat  models 

2.  Calculating  the  same  for  hot-gas  type  false  threats 
(e.g.,  industrial  emissions) 

3.  Calculating  the  spectral  signatures  of  continuum-like 
false  threats  in  the  region  of  the  blue  spike  (i.e., 
use  the  model  to  calculate  the  fully  resolved  atmos- 
pheric spectral  transmittance) 

4.  Analyzing  lower  resolution  false  threat  spectral  sig- 
nature data. 
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SECTION  VI 

INDUSTRIAL  SITE  BACKGROUND  MODELS 


Procedures  have  been  developed  for  representing  industrial  site 
processes  that  result  in  intense  infrared  spectral  signatures  in  the  3.5  to 
5. 0-micrometer  region.  In  the  majority  of  this  work  we  have  relied  on  high 
resolution  Fourier  Transform  Spectrometer  (FTS)  measurements  obtained  by  the  Air 
Force  Geophysics  Laboratory  (AFGL)^  and  by  the  Environmental  Protection  Agency 
(EPA) . While  both  of  these  laboratories  provided  excellent  calibrated  spectra  of 
industrial  processes,  comparable  data  on  physical  parameters  (e.g.,  exact  sizes, 
temperatures,  compositions)  has  generally  not  been  available.  Thus  in  most  cases 
it  has  been  necessary  to  rely  on  estimates  or  perform  parametric  calculations  to 
obtain  physical  parameters  in  modeling  the  spectra.  However,  no  unreasonable 
assumptions  were  required,  and  we  feel  that  the  procedures  developed  and  site 
parameters  inferred  in  this  program  have  been  adequately  validated  for  use  in 
analyzing  the  discrimination  capabilities  of  IR  threat  warning  receivers. 

6.1  RADIATIVE  MODELING  OF  HOT  GAS  FALSE  THREAT  SOURCES 


Our  primary  data  base  for  analyzing  and  validating  hot  gas  emission 

s 

sources  has  been  the  high  resolution  FTS  spectra  supplied  by  William  Herget  of 

22  —1 
EPA  . These  measurements  consist  of  high  resolution  (0.25  — 1.0  cm  ),  narrow 

field-of-view  (0.5  degree),  short  range  (0.5  — 1.0  km)  spectra  cr  various  hot  gas 

emission  sources.  Included  in  the  measurements  are  gas  burnoff  flares,  boiler 

stack  plumes,  refinery  tower  emission  plumes,  aircraft  plumes  and  other  similar 

sources.  Since  the  EPA  interferometer  was  ground  based,  and  the  instrument  field 

of  view  is  comparable  to  the  plume  sizes  of  these  sources,  the  measurements 

represent  the  source  signatures  with  little  or  no  background  interference.  We 

have  used  this  data  as  the  primary  resource  in  validating  the  hot  gas  emission 

models. 
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The  computer  codes  used  in  the  analysis  of  the  hot  gas  emission  spectra 

are  ATLES  and  BLUESPIKE.  The  ATLES  code  is  a moderate  resolution  band  model  code 

23 

developed  by  Stephen  Young  of  Aerospace  Corporation  . This  is  probably  the  most 
sophisticated  band  model  code  (and  associated  parameter  set)  currently  available. 
In  the  course  of  the  present  study  a number  of  minor  modifications  and 
improvements  were  made  to  the  ATLES  code. 

Four  modifications  have  been  made  to  ATLES  which  />■  Jtantially  increase 
its  capability  and  computational  efficiency: 

1.  The  program  flow  has  been  altered  such  that  PTC  (pres- 
sure, temperatures  and  composition)  variation  along  an 
emission  path  need  only  be  calculated  once  rather  than 
recalculated  at  each  spectral  position,  as  is  done  in 
the  original  code. 

2.  Pointers  to  desired  band  model  parameters  are 
calculated  by  subroutine  SPCLOP  at  the  beginning  of 
each  run  rather  than  recalculated  on  each  call  to 
BPARAM. 

3.  Up  to  five  species  may  be  used  in  the  absorption  and/or 
emission  path  segments. 

A.  Transmittance  and  transmittance  derivatives  can 
optionally  be  calculated  for  water  and/or  nitrogen 
continuum  for  the  absorption  path  only.  Coupled  path 
calculations  assume  that  the  continuum  transmittance 
derivatives  go  to  zero  upon  entering  the  emission  path. 

The  BLUESPIKE  code  is  described  in  Section  V.  This  model  predicts  line- 
by-line  CO^  emission  and  atmospheric  absorption  in  the  blue  spike  region  between 
2375  end  2400  cm 
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6.2 


STACK  PLUMES 


The  most  commonly  encountered  hot  gas  emission  source  will  be  various 
types  of  8 tack  plumes.  Some  examples  of  these  are  the  EPA  spectra  shown  in 
Figures  38,  39,  40  and  41.  The  first  two  spectra  present  the  emission  of  refinery 
stacks.  In  these  spectra  nearly  all  of  the  spectral  structure  is  the  result  of  CO2 
emission.  The  only  exceptions  seen  are  isolated  H„0  emission  lines  beyond 
2000  cm  and,  in  the  case  of  the  sulfur  recovery  unit  exhaust,  the  SO,,  emission 
centered  at  2500  cm  (Figure  39).  Both  of  these  spectra  represent  relatively 
small  plumes  estimated  to  be  less  than  one-half  meter  in  diameter. 

Spectra  40  and  41  are  typical  of  boiler  exhaust  plumes.  Figure  40  is 
the  exhaust  from  a power  plant  boiler  at  Duke  University.  The  familiar  blue  spike 
and  red  spike  plume  structures  are  readily  apparent.  Also  visible  is  the  CO 

fundamental  band  centered  at  2140  cm  and  a C0„  band  with  its  Q-branch  at 

-1  1 
2175  cm  . This  red  wing  COj  band  is  prominent  in  large  optical  depth  COj  plumes 

and  appears  to  be  a good  discriminant  between  boiler  plumes  and  aircraft  or 

missile  plumes. 

The  CO  boiler  plume  spectra  in  Figures  41  and  42  result  when  CO  (which 
is  a byproduct  of  another  refinery  process)  is  burned  to  heat  a boiler.  Clearly, 
this  spectra  is  similar  to  the  other  three  even  though  it  results  from  purely  C0o 

Am 

emission.  Superimposed  on  Figure  42  is  an  ATLES  computation  of  the  predicted 
spectra.  Except  for  the  obvious  difference  in  resolution,  the  ATLES  spectra  is 
generally  a good  representation  of  the  measured  result.  There  is  a discrepancy  in 
the  weaker  CO^  bands  between  2040  and  2140  cm  *,  the  reason  for  this  discrepancy 
is  unclear.  Two  possibilities  are  that  it  either  results  from  inaccurate 
assumptions  of  the  plume  optical  thickness  or  from  inaccurate  band  model 
parameters  in  this  region. 

To  get  a more  general  representation  of  the  variation  which  is  to  be 
expected  with  stack  plumes  we  have  made  parametric  calculations  of  the  peak 
radiances  predicted  for  the  blue  spike  (near  2390  cm  ),  the  red  spike  (near 
2210  cm  *)  and  the  red  wing  (near  2050  cm  *).  These  results  are  shown  in 
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Figure  41.  Measured  Emission  Spectra  of  the  Exhaust  Plume  of  a CO  Fueled  Boiler  Stack 
(Lower  Curve)  and  the  External  Stack  Wall  (Upper  Curve) 


Figures  43,  44  and  45.  As  is  seen  from  those  figures,  the  radiance  values  of  all 
three  bands  increase  similarly  with  temperature.  This  prediction  is  in  agreement 
with  the  behavior  exhibited  by  all  other  spectra  examined  in  this  study.  Simi- 
larly, for  the  temperatures  considered  here,  the  blue  spike  and  red  spike  peaks 
are  relatively  independent  of  CC^  optical  thickness.  Thus  the  relative  peak 
intensities  in  the  red  spike  and  blue  spike  cannot  be  used  to  infer  either  plume 
temperature  or  plume  optical  thickness.  If,  however  we  examine  Figure  45  it  is 
seen  that  the  red  wing  peak  radiance  is  approximately  linearly  related  to  the 
plume  CO2  optical  depth  for  all  three  temperatures  considered. 

Since  aircraft  or  missile  plumes  are  typically  much  thinner  optically 
than  power  plant  stack  plumes,  the  ratio  of  the  red  wing  to  red  spike  radiation  may 
be  a viable  discriminant  between  threats  and  stack  plumes.  An  indication  of  the 
validity  of  this  type  of  discrimination  can  be  9een  by  examining  the  measured 
aircraft  plume  shown  in  Figure  46.  Compare  the  noticeable  lack  of  a red  wing 
signature  in  that  figure  with  the  prominent  red  wing  in  Figures  39  , 40  and  41. 

6 . 3 FLARES 

Flares  (open  gas-fueled  flames)  are  widely  used  in  the  chemical  and 
steel  industries  to  dispose  of  flammable,  unwanted  gaseous  waste  products.  They 
are  also  used  in  many  locations  to  burn  off  unusable  natural  gas  from  oil  fields. 
EPA  has  provided  us  with  two  examples  showing  the  spectral  structure  of  flares; 
these  are  presented  in  Figures  47  and  48.  The  flare  represented  in  Figure  47  is 
used  to  burn  off  waste  gases  at  a refinery,  and  that  represented  in  Figure  48  is 
used  to  burn  off  vinyl  chloride  related  wastes  at  a chemical  plant.  The 
similarities  between  these  two  spectra  is  obvious,  particularly  with  respect  to 
the  shape  and  relative  magnitudes  of  the  red  spike  and  blue-spike.  Generally,  the 
width  of  the  redspike  and  the  pinpoint  extending  from  the  blue  spike  are  very 
characteristic  of  these  very  high  temperature  sources. 

The  spectra  shown  in  Figure  48  also  contains  other  molecular  emission 
besides  the  expected  HjO  and  CC^.  The  overtone  band  of  CO  is  clearly  visible  as 
uniformly  spaced  emission  lines  between  2040  and  2120  cm  *.  Also  visible  is  HCL 
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Figure  43.  Computed  Radiance  of  the  Blue  Spike  Peak  for  an  Aysutned 
Stack  Plume  Viewed  through  a One  Kilometer  Atmospheric 
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Figure  44.  Computed  Radiance  or  the  Red  Spike  Peak  for  an  Assumed 
Stack  Plume  Viewed  Through  a One  Kilometer  Atmospheric 


emission  lines  with  regular  spacing  of  approximately  20  cm  ^ between  2500  and 
2900  cm  ^ . This  molecule  results  from  the  combustion  of  the  chlorine  compounds  by 
the  flare. 


The  peak  intensities  of  the  HCL  lines  were  used  to  deduce  a rotational 
temperature  of  2400  K for  the  flame.  Using  this  temperature  and  the  emission  and 
absorption  path  quantities  given  in  Table  7,  calculations  were  made  of  the 
predicted  emission  spectra  using  ATLES . As  shown  in  Figure  48  the  absorption  from 
atmospheric  N^0  must  be  included  to  obtain  shape  agreement  between  the  prediction 
and  the  measurement.  In  general,  the  agreement  between  the  measurement  and  the 
prediction  is  good.  Of  course,  the  prediction  does  not  show  the  line  structure  of 
either  the  CO  or  the  HC1  bands.  The  sharp  depression  at  2500  cm  * results  from  a 
gap  in  the  H„0  band  model  parameters  caused  by  the  merging  of  two  not-quite- 
overlapping  band  model  parameter  sets. 

It  is  also  apparent  from  the  comparison  that  bend  model  parameters 

currently  available  for  ATLES  are  inadequate  for  modeling  the  bluespike  at  high 

temperature.  The  prediction  of  the  BLUESPIKE  computer  code  when  applied  to  this 

measurement  is  shown  in  Figure  49.  While  this  result  is  significantly  better  than 

given  by  ATLES  it  is  still  somewhat  inadequate.  In  fact,  for  most  situations  the 

BLUESPIKE  predictions  will  be  better  than  indicated  by  Figure  49,  since  with 

__  1 

longer  atmospheric  paths  the  lines  below  2390  cm  * will  be  strongly  attenuated  by 
CO2  absorption. 

Figure  50a,  b,  c and  d show  the  contributions  of  the  various  active 
sp  :ies  in  the  flare  as  computed  by  ATLES.  These  calculations  do  not  contain  NjO 
in  the  atmospheric  path.  One  fact  to  notice  is  the  comparison  between  the 
measured  high  resolution  spectra  (Figure  48a)  and  the  computed  ATLES  spectra.  In 
particular,  note  the  strikingly  simple  manner  with  which  the  constituent  molecules 
can  be  identified  from  the  measured  spectra  as  compared  with  the  computed  low 
resolution  spectra  where  everything  merges  together  into  a continuum. 
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TABLE  7 

PHYSICAL  CONDITIONS  USED  FOR  CHEMICAL  BURNOFF  FLARE  MODEL 


EMISSION  PATH 


Thickness : 
Temperature: 
Total  Pressure: 


1.0  atm 


Partial  Pressures  of  Active  Species 
H20: 


C02: 
CO: 
HCL : 


ABSORPTION  PATH 


50  cm 
2400  K 


0.015  at.a 
0.010  atm 
0.002  atm 
0.001  atm 


Length: 

Temperature : 

Total  Pressure:  1.0  atm 

Partial  Pressure  of  Active  Species 


H20: 

C°2  ‘ 
CO: 

N20: 


1 km 
272.2  K 


4.344  E-03  atm 
0.00033  atm 
7.5  E-08  atm 
2.8  E-07  or  0 atm 


Water  and  Nitrogen  Continuum  is  included. 
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Figure  49.  Blue  Spike  Radiance  Calculated  by  the  BLUESPIKE 
Computer  Code  Compared  With  Flare  Measurements 
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HEATED  GRAV’30DY  RADIATIVE  SOURCES 


The  measurements  available  for  graybody  sources  were  obtained  by  AFGL 
over  industrial  areas  of  the  Southeast  and  Midwest  in  June  1978.  The  measurements 
were  taken  from  the  AFGL  KC-135  aircraft  primarily  with  the  number  103  FTS  system, 
whi.ch  has  a one— degree  f ield-of-view.  The  measurements  were  made  wV  i le  the 
aircraft  circled  the  target  at  a range  of  approximately  3500  meters.  Thus,  the 
f ield-of-view  projects  a 60— meter  diameter  spot  on  the  ground.  Since  this  is 
significantly  larger  than  the  diameters  of  most  sources,  the  resulting  spectra 
contain  contributions  from  the  backgrounds  as  well  as  localized  hot  sources. 
Auxiliary  data  was  available  for  some  of  the  measurements  in  the  form  of  thermal 
imagery,  movies  and  still  photographs.  However,  this  data  was  of  limited  utility 
in  analysis  of  the  spectral  data.  The  radiometric  imagery  was  used  to  locate  the 
sources  and  approximate  areas  of  the  hot  spots,  and  the  photography  was  useful  for 
tho  same  purpose.  However,  because  of  blooming  in  the  thermal  images  and  the 
difficulty  in  identifying  hot  spots  from  the  photography,  the  source  areas  could 
only  be  estimated.  The  radiometric  imagery  was  uncalibrated  and  therefore  was  not 
useful  for  estimating  source  radiance  or  intensity  values. 

The  complete  set  of  industrial  site  spectra  obtained  from  AFGL  are  given 
in  Appendix  D.  (AFGL  obtained  these  spectra  in  support  of  the  IR  Signature  Study, 
and  does  not  plan  to  publish  them  as  an  AFGL  report.)  The  corresponding 
radiosonde  data  are  included  in  Appendix  B.  The  spectra  cover  the  wavelength 
range  measured  by  AFGL;  however  the  present  study  only  examined  the  35  to 
5.0-micrometer  region.  -Although  significant  solar  scattering  was  observed  at  the 
shorter  wavelengths,  no  solar  contributions  were  observed  in  the  3.5  to 
5.0-micrometer  region.  Therefore  the  present  analysis  ignored  solar  scattering. 

Figures  51,  52  and  53  are  typical  examples  of  the  measured  spectra 
of  various  graybody  sources.  The  most  conspicuous  structure  is  the  dip  in  the 
radiance  due  to  the  CO^  absorption  band  between  2220  and  2380  cm  *.  Another 
prominent  feature  observable  in  all  of  the  graybody  spectra  is  the  small  emission 
peak  near  2220  cm  ^ resulting  from  the  higher  transmittance  between  the  P and  R 
branches  of  the  N^O  atmospheric  absorption  band.  The  appearance  of  this  feature 
makes  many  of  the  graybody  spectra  look  similar  to  the  red  spike  structure  seen  in 
molecular  spectra  in  this  wavelength  region. 
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We  carefully  analyzed  a number  of  the  AFGL  spectra  in  two  respects:  we 
estimated  the  temperatures  and  emissivity-areas  of  the  sources,  and  we  also 
obtained  the  apparent  intensity  of  the  sources. 

6.4.1  Temperatures  and  Emissivity-Areas  of  Graybody  Sources 

4 

LOWTRAN  4 was  used  in  the  emission  mode  to  analyze  the  measured  graybody 
spectra.  The  following  procedure  was  used:  first,  the  measured  data  was  plotted 
on  semi-log  ordinste  scales.  Measured  spectra  which  showed  molecular  emission 
were  put  aside  for  later  analysis.  Next,  a number  of  LOWTRAN  radiance' 
calculations  were  made  assuming  source  temperatures  between  400  and  1000  K and 
unit  emissivity.  These  calculations  were  also  plotted  in  semi-log  form.  An 
example  computation  with  a source  temperature  of  600  K is  shown  in  Figure  54, 
With  the  calculations  in  this  form  a vertical  offset  in  the  graphs  is 
approximately  equivalent  to  changing  the  emiss ivity-area  of  the  source.  (This 
would  be  rigorously  true  except  for  the  contributions  of  path  radiance  which  are 
negligible  in  the  infrared  window  regions.)  Finally,  by  overlaying  the  LOWTRAN 
calculations  on  the  measured  spectra,  the  approximate  temperature  and  emissivity- 
area  of  the  graybody  source  was  obtained.  In  this  process  the  temperature  was 
determined  by  selecting  the  LOWTRAN  spectra  which  best  fit  the  shape  of  the 
measured  spectra  in  the  3.5  to  4.1-micrometer  window  region  regardless  of  the 
level  of  the  curve.  The  temperature  of  the  best-fit  LOWTRAN  spectrum  was  then 
taken  as  the  approximate  temperature  of  the  source.  We  estimate  the  accuracy  of 
this  procedure  is  approximately  J^lOO  K.  The  offset  of  the  best-fit  LOWTRAN 
calculation  when  overlayed  on  the  measured  spectra  gives  the  source  emissivity- 
area  as  a ratio  of  the  total  area  within  the  f ield-of-view  of  the  instrument.  A 
similar  procedure  has  also  bsen  defined  based  on  the  computer  algorithm  described 
in  Appendix- E.  The  emissivity-areas  inferred  by  either  procedure  are  critically 
dependent  on  the  inferred  temperatures  of  the  sources;  however,  if  the  emissivity- 


area  and  temperature  pairs  are  used  together,  results  consistent  with  the  j 

; S 

measurement  will  be  obtained.  ! J 

i 


The  temperatures  and  emissivity-areas  obtained  are  given  in  Table  8. 
Not  all  of  the  spectra  were  analyzed  in  this  manner,  but  the  values  given  in  the 
tabl?  encompass  the  range  of  sources  measured.  Of  course  the  measurements 
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represent  complex  scenes  and  thus  the  temperatures  must  be  considered  effective 
temperatures,  which  will  be  near  the  peak  temperature  observed  within  the  scene. 
In  general,  similar  physical  sources  are  inferred  to  have  similar  temperatures. 
Any  scene  which  was  identified  as  having  molten  or  "orange"  metal  visible  was 
anlayzed  to  have  a temperature  of  800  K or  slightly  higher.  Hot  coke  ovens  and 
other  hot  buildings  have  temperatures  ranging  from  500  to  600  K.  The  emissivity- 
areas  are  typically  20—200  square  meters;  these  are  small  compared  to  the  instru- 
ment f ield-of-view  or  the  physical  sizes  of  the  buildings.  These  values  are  also 
in  agreement  with  areas  obtained  from  the  infrared  imagery  for  cases  where  imagery 
was  available,  The  small  effective  sizes  indicate  the  difficulties  in  obtaining 
the  emitting  areas  from  visible  photography,  which  does  not  distinguish  the  hot 
areas  from  the  remainder  of  the  structure. 

6.4.2  Graybody  Source  Intensities 

The  apparent  radiant  intensities  of  sources  were  also  obtained  directly 
from  the  measured  spectra.  These  results  were  obtained  for  the  mid-IR  window 
(3.5  — 4.15  micrometer),  the  blue  spike  (4.15  - 4.20  micrometer")  and  the  red 
spike  and  wing  (4.4  — 5.0  micrometer  ) spectral  regions;  these  results  are  shown 
in  Table  9.  These  values  were  obtained  by  integrating  the  measured  spectra  over 
the  wavelength  bands  of  interest  and  subtracting  similarly  integrated  background 
intensities.  Since  these  values  were  obtained  directly  from  the  measured  spectra 
they  are  independent  of  the  temperature  and  emissivity-area  results  and  the 
uncertainties  therein.  Because  the  ranges  and  background  radiances  are  somewhat 
uncertain,  these  values  have  associated  uncertainties  of  perhaps  +25  percent,  but 
they  do  show  the  typical  intensities  expected  for  these  types  of  sources. 


The  intensity  ratios  given  in  Tables  8 and  9 were  computed  in  the 
hope  that  a unique  indication  of  a graybody  source  could  be  obtained  from  these 
ratios.  Unfortunately,  it  was  found  that  the  ratios  varied  considerably  from 
source  to  source;  probably  because  of  the  sensitivity  of  the  blackbody  function  to 
temperature  within  the  region  examined. 
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6,4.3 


Combined  Graybody  and  Gaseous  Emission  Sources 


A number  of  Che  AFGL  measurements  include  both  graybody  emission  and 
gaseous  emission  sources  within  the  f ield-'of-view.  Examples  of  this  type  of 
source  are  shown  in  Figures  55,  56  and  57.  The  measured  stack  flare 
emission  spectrum  shown  in  Figure  55  was  one  of  the  strongest  sources  of  gaseous 
emission  as  well  as  the  brightest  radiator  of  all  sources  measured.  As  can  be  seen 
from  Figure  55,  the  "red  spike"  region  is  very  similar  in  structure  to  some  of 
the  hotter  graybody  sources.  (Compare  with  Figure  53.)  This  is  because  the 
primary  structure  in  the  red  spike  region  for  high  temperature  or  graybody  sources 
results  from  the  absorption  structure  in  the  atmosphere.  Some  of  the  molecular 
line  complexes  in  the  region  show  a "W"  shape  which  may  be  characteristic  of  a 
gaseous  emission  source.  This  type  of  structure  is  not  evident  in  the  graybody 
measurements . 

The  blue  spike  is  clearly  visible  in  all  of  the  molecular  emis  ion 
spectra  at  2395  cm  . This  feature  was  clearly  distinguishable  in  every  spectrum 
where  molecular  emission  was  expected  and  thus  appears  to  be  an  unambiguous 
indication  of  molecular  emission. 

These  composite  spectra  were  also  analyzed  in  a manner  similar  to  that 
described  for  the  graybody  spectra.  The  graybody  component  was  extracted  by 
concentrating  on  the  3.5  to  4 . 15-micrometer  window  region  where  little  or  no 
molecular  emitting  areas  were  inferred  from  the  expected  temperatures  of  the 
sources  and  the  peak  radiances  of  the  red  spike  and  blue  spike  above  the  graybody 
levels.  These  areas  must  be  considered  crude  estimates;  however  in  the  case  of 
the  "flame  from  stack"  the  inferred  area  agrees  closely  with  the  value  obtained 
from  the  thermal  imagery.  The  results  are  given  in  Table  9. 

The  separation  of  the  emission  into  graybody  and  gaseous  emission 
components  appears  to  be  valid  except  for  the  case  of  the  stack  flare.  For  that 
case,  the  continuum  emission  appears  to  have  a slope  much  too  high  to  be  explained 
as  heated  metal  emission.  It  is  our  hypothesis  that  this  results  either  from 
particulates  in  the  flame  or  perhaps  HjO  emission  from  weak  hot  band  lines.  In 
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Figure  55.  Measured  Apparent  Radiance  of  Stack  Flare  (Gas  Burnoff)  at  Gary,  IN 


wpwww.5w.Tt  ~f”r- 


w W— ,T/>vu;,-i,^-»-»-T7 »».— rwr.-.s;  ■■^wr-  rv/  *^7»TV!W«egF!'1'*'r 

r 


general,  the  roolecular  emission  from  the  flare  or  flame  sources  dominates  the 
graybody  emission;  however,  except  for  the  blue  spike,  the  molecular  emission  from 
the  stack  sources  is  dominated  by  the  graybody  emission  from  the  hot  stack  or 
associated  structures. 


SECTION  VII 

REFUSE  FIRE  FALSE  THREATS 


Burning  refuse  material  represents  a serious  and  common  false  threat 
source  over  battlefield  environment?.  The  only  direct  spectral  radiance  measure- 
ment located  for  this  type  of  threat  is  a spectrum  of  rice-stubble  burnoff 
measured  by  AFGL  during  a flight  over  the  San  Joaquin  Valley  in  California.  nis 
spectrum  is  shown  in  Figure  58.  With  the  limited  resolution  of  the  data  we  have 
not  been  able  to  identify  any  molecular  emission  except  a very  weak  blue  spike  CO^ 
emission.  Other  molecular  emission  may  be  present;  however,  a careful  computer 
analysis  would  probably  be  required  to  identify  it. 

The  window  region  between  3.5  and  4.2  micrometers  in  the  measured 
spectrum  shows  evidence  of  graybody  emission  similar  to  that  seen  in  the  hot 
industrial  sources.  We  have  analyzed  the  spectrum  in  a manner  similar  to  that 
described  in  the  previous  section.  The  results  indicate  an  effective  graybody 
temperature  of  575  ± 100  K and  an  effective  emissivity-area  of  7.5  percent  of  the 
instrument  f ield-of-view.  Given  the  AFGL  estimate  that  the  visible  flame  filled  1 
to  2 percent  of  the  f ield-of-view,  this  result  for  the  effective  emissivity-area 
of  the  source  seems  reasonable. 

A set  of  measurements  of  the  absorption  of  the  products  of  refuse  fires 
was  located.  This  data  is  summarized  in  Table  10  and  the  measured  spectra  are 
shown  in  Figures  59  to  63.  These  measurements  were  made  by  burning  the 
substance  in  question  in  the  laboratory  and  collecting  the  gas.s  and  particulates 
produced  in  cells  or  on  filters. 

While  it  is  doubtful  that  accurate  source  emission  intensities  could  be 
obtained  from  this  set  of  data,  it  does  provide  useful  qualitative  information. 
For  most  materials  and  spectral  regions,  the  particulate  absorptivities  are  much 
larger  than  the  molecular  absorptivities;  but,  there  seems  to  be  a window  in  the 
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Figure  62.  Absorptivity  of  the  Gaseous  Products  of  Burning 
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particulate  absorption  between  1800  and  2800  cm  * in  the  important  molecular 
signature  bands.  All  of  the  refuse  material  spectra  show  strong  hydrocarbon  and 
CO  bands  in  contrast  to  the  industrial  site  spectra.  Thus,  it  seems  that  the 
presence  of  unburned  hydrocarbons  and  particulates  might  be  a good  discriminant 
for  refuse  fires. 

In  order  to  calculate  the  spectral  radiance  of  these  sources  it  is 
necessary  to  obtain  estimates  of  flame  temperature  and  the  size  of  the  radiating 
volume,  along  with  the  size  and  spectral  extinction  coefficient  of  the  smoke 
produced  by  the  fire. 

Several  estimates  of  fire  temperatures  are  available  from  USDA  Forest 

2 5 26 

Service  Measurement  Programs.  Both  Palmer  and  ..irsch  et  al  report  typical 
temperatures  of  1000  K for  large  free-burning  fires. 

The  size  of  the  radiating  flame  volume  can  be  determined  semi-empiri- 
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cally  from  measurements  by  Thomas  . The  basic  relationship  between  the  flame 
length  (L)  and  the  flame  width  (D)  can  be  related  to  the  material  properties  of  the 
blaze  by  the  expression: 

L/D  - 42  |ro  " (PQ  /gD)_I 


where  ra"  » mass  fuel  consumption  rate  per  unit  area,  pQ  ■ combustion  gas  density, 
and  g ■ acceleration  of  gravity.  The  magnitude  of  the  exponent  T)  is  a function  of 
burning  material.  A typical  value  for  n is  0,61.  Normally  L/D  will  vary  between  2 
and  2.5  for  kerosene  pool  fires  that  are  several  meters  in  diameter. 

The  size  of  large  fire  plumes  ht,s  been  studied  extensively  in  the  fluid 

mechanics  literature.  A useful  relationship  between  the  height  (H)  of  the  smoke 
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and  the  rate  of  heat  production  (Q)  by  the  fire  was  derived  by  Morton  et  al  for  a 
stably  stratified  atmosphere  with  a known  temperature  gradient: 
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H - 31  (l+m)“3/8  Q1/4, 


where  Q is  the  rate  of  heat  production  of  the  source  in  kilowatts  and  m is  the 
ratio  of  the  vertical  gradient  of  the  absolute  temperature  to  the  atmospheric 
lapse  rate.  Applying  this  expression  to  an  ICAO  standard  atmosphere,  Morton  finds 
that  H varies  from  450  meters  for  a small  bonfire  to  as  much  as  2200  meters  for  a 
burning  forest. 

The  attenuation  of  the  flame  radiation  by  the  smoke  can  be  estimated 
from  Mie  scattering  theory.  The  primary  material  in  the  smoke  cloud  is  partially 
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combusted  hydrocarbon  and  carbon  in  a wide  variety  of  particle  sixes.  PalmeriJ 
indicates  that  smoke  particles  from  free-burning  cellulose  fires  are  generally 
less  than  1 micrometer  in  radius.  The  maximum  in  the  number  density  size  dis- 
tribution occurs  around  0.3  micrometer  near  the  bottom  of  the  smoke  column.  The 
concentration  of  smoke  particles  varies,  but  a reasonable  order  of  magnitude 
estimate  is  10  partioles/craJ. 

The  exact  Lorentz-Mie  calculation  of  the  smoke  spectral  extinction 

coefficient  may  often  be  replaced  by  approximate  analytic  solutions.  A very 

29 

useful  result  is  due  to  Casperson  who  generalized  the  anomalous  diffraction 
method  Van  de  Hulst  to  include  absorbing  particles.  His  expression  for  the 
exponential  extinction  coefficient  is: 


Y(X)  * (fTa* 


2 N/n2)  I" 

max  o i_ 


2(n+2)(n+l)  + Yj-Y2  sin  Y^-Y^  cos  Y^ 


where  the  Y^(A)  are  given  in  Appendix  F , Nq  is  the  particle  number  density,  amax 
is  the  radius  at  which  the  particle  distribution  peaks,  and  n is  the  inde..  for  a 

power- law-exponent ial  particle  size  distribution  function.  A plausible  value  for 

■ . . 29 

n that  is  typical  of  many  smoke  ensembles  is  n - 4.6.  Casperson  gives  several 

specific  examples  of  Y(A)  for  absorbing  water  fogs  and  for  non-absorbing  smoke 

ensembles  with  an  index  of  refraction  equal  to  1.5.  The  calculation  of  Y(A)  for 

30 

carbon  smoke  involves  utilizing  the  optical  constants  of  soot  which  are  quite 
different  from  fogs. 


- A -t;  ..Wlfctw.  .,1 
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SECTION  VIII 


MUNITION  AND  SMOKE  FALSE  THREATS 

The  muzzle  flash  of  artillery,  explosive  shell  bursts  and  screening 
smoke  dispersions  may  also  cause  false  threat  events  in  warning  receivers. 
Table  11  summarizes  the  available  information  on  these  types  of  events. 

A substantial  amount  of  radiometric  data  and  some  moderate  resolution 
(32  cm  ) spectrometric  data  is  available  on  muzzle  flash  signatures.  Typical 
muzzle  flash  signatures  have  a blue  and  red  spike  CO^  signature  similar  to  a 
rocket  exhaust  but  with  an  underlying  continuum  emission  attributable  to 
particulates  in  the  plume.  Typical  duration  of  the  flash  is  5 to  100  milli- 
seconds. Because  of  classification,  none  of  the  measurement  results  are  included 

31  32 

here;  however,  they  are  available  in  the  Pardes  and  Robertson  references.  ’ 

No  shell  burst  data  has  been  located  for  the  spectral  region  considered 
in  the  present  study.  Radiometric  measurements  in  the  visible  and  near-IR 
indicate  that  shell  bursts  intensities  are  somewhat  less  than  muzzle  flashes  of 
the  related  weapon  and  have  pulse  widths  of  from  1 to  5 milliseconds. 

Screening  smokes  may  be  widely  used  in  future  ground  engagements. 
Hexachlcroe thane  (HC)  and  white  phosphorous  (WP)  are  generally  disseminated  by 
artillery  fire  and  produce  smoke  through  burning.  Thus,  the  initial  smoke 
temperatures  are  fairly  high.  As  the  smoke  cloud  diffuses  in  the  wind,  it  is 
quickly  cooled  by  mixing  with  the  air  and  within  20  to  100  feet  of  the  wind,  it  has 
cooled  to  within  a few  degrees  of  ambient.  Sulphur  trioxide  (FS)  and  Fog  Oil  smoke 
are  generally  disseminated  by  a spray  so  they  typically  have  lower  temperatures 
than  the  pyrotechnic  smokes.  In  the  visible  and  near-IR,  these  smokes  screen 
through  aerosol  scattering.  In  the  mid-IR  region  of  interest  here,  scattering  and 
aerosol  absorption  are  comparable  effects  (see  Table  12).  Thus,  when  observed  in 
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TABLE  11 

POTENTIAL  MUNITION  FALSE  THREAT  SOURCES 


TABLE  11 


f 


-PWT 


i 

\ 

i 

5 
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TABLE  12 

CALCULATED  FS  SMOKE  ATTENUATION  FOR  THE  NOMINAL. 

SIZE  DISTRIBUTION  AND  50X  ACID 

DROPLET  CONCENTRATION^ 

Total 


Laser 

Wavelength 

Scattering 

Coefficient 

Absorption 

Coefficient 

Attenuation 

Coefficient 

(|jni) 

(m2/g) 

(m2/g) 

(m2/g) 

0.  55 

3.  526 

0.  000 

3.  526 

1,06 

2.  104 

0.  000 

2.  104 

3.  39 

0.237 

0.  216 

0.  453 

3.  80 

0.  187 

0.  162 

0.  349 

10.60 

0.  058 

0.  146 

0.  204 

EXTINCTION  COEFFICIENT  (MVGM) 


r»wgX!^'»r  "Ofr*-  -?r&n-rwa  - ■nr^i,«  ■-»>  r i 


EXPERIMENTAL  DATA:  RP 

CL  = .945  GM/M2 
MMD  = 1.13  um 
on  = 1.40 

Drop  Cone.  = 95.2%  H3PO4 
RH  = 58% 


* 1 1 I | ! I i 1 | 1 1 i I [ 1 " I ""I  "T j I I T '■  I ’ ’1 

5 3.0  3.5  4.0  4.5  5.0 

WAVELENGTH  (ym) 
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Figure  66.  Measured  Extinction  of  WP  (RP)  Smoke 


Absorptivity  or  standard  deviation 


this  region,  the  screening  smokes  will  appear  primarily  as  continuous  emitters  At 
the  temperature  of  the  smoke.  However,  as  shown  by  Figure  58,  there  is  a 
significant  molecular  content  to  the  smokes,  and  thus  molecular  absorption  or 
emission  may  be  superimposed  on  the  graybody  signatures. 

Colored  signaling  smokes  are  also  included  in  Table  11  and  one  spectral 
example  is  shown  in  Figure  69.  While  this  type  of  smoke  is  not  expected  to  be  a 
large  component  in  the  battlefield,  the  spectra  is  interesting  because  it  shows 
the  large  variety  of  molecular  constituents  which  can  be  present  in  these  sources. 
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APPENDIX  A 

SCATTERING  PARAMETERS  FOR  CLOUD  RADIANCE 


These  scattering  parameters  were  computed  by  Eric  Shettle  of  AFGL  fop 
log-normal  distributions  of  water  and  ice  spheres  (see  Section  3).  Each  of  the 
10  tables  represents  a particular  type  of  particle  (water  or  ice)  and  a particular 
mean  radius  for  the  log-normal  distribution.  The  type  and  radius  (in  micrometers) 
is  given  by  the  first  line. 

The  remainder  of  each  table  consists  of  8 columns,  which  give  the 
wavelength  in  micrometers;  the  real  and  imaginary  indices  of  refraction;  the 
extinction,  scattering  and  absorption  coefficients  cIq,  $q,  Kq  (each  in  units  of 
cm  1 and  representing  a number  density  of  1 per  cubic  centimeter);  the  single 
HJscattering  albedo  and  the  anisotropy  parameter  g.  Only  Kq,  and  g are 

required  by  the  cloud  radiance  model  described  in  Section  3.  The  water  and  ice 

. . . 3 
indices  of  refraction  were  taken  from  Kale  and  Querry. 


Table  A-l 


WATER 


wavlgh 

REFR- 

IMOICES 

EXT9C-CDEF 

SCATT-COEF 

ABSRP-COFF 

4LRDO 

G 

?.600o 

1.242 

-.0032 

8.3452E-07 

7 . 7 8 5 1 E-07 

5.6009^-08 

.933 

.8721 

2.6500 

1.219 

-.0067 

8.2564E-07 

7.2538E-07 

1 • 0026E-0  7 

.879 

.8970 

2.7000 

1.188 

-.0190 

7,(.865E-07 

5.9340E-07 

1 .954SE-07 

.752 

.9315 

2.7500 

1.157 

-.0590 

7.0340E-07 

3.9904E-07 

3.0434E-07 

.567 

.9575 

2.0000 

1.142 

-.1150 

6.6156E-07 

3.1793E-07 

3.4363E-07 

.481 

.9606 

2.8500 

1.149 

-.1550 

6.7170E-07 

3.1073E-07 

3.6097E-07 

• 463 

.9533 

2.9000 

1.201 

-.2680 

7,061 6E-07 

3.3340E-07 

3.7276E-07 

.472 

.9393 

2. 0500 

1.292 

-.2980 

7.3636F-07 

3.5571E-07 

3.8113E-07 

• 483 

.9281 

T.nOOO 

1.  371 

-.2720 

7.5B45E-07 

3.7164E-07 

3.8681E-07 

.490 

.9277 

3.0500 

1.426 

-.2400 

7 , 7337E-07 

3.9345E-07 

3.8992E-07 

.496 

.9150 

3 • i 0 0 0 

1.467 

-.1920 

7.8771E-07 

3.9704E-07 

3.9065E-07 

.504 

.9086 

3.1500 

1.483 

-.1350 

0.O196E-O7 

4. 1740E-07 

3.8456E-0? 

.520 

.897? 

3.0000 

1.478 

-.0924 

8,1 382E-07 

4.4591E-07 

3.6790E-07 

.548 

.8809 

3.2500 

1.467 

-.0610 

8.2463E-07 

4. 9794F-07 

3.3669E-07 

.592 

.8585 

3.3000 

1*450 

-.0368 

8.3517F-07 

5.5204E-07 

2.B312E-07 

.661 

.831? 

3.3500 

1.432 

-.0261 

8.420BE-07 

6.0195E-07 

2.4014E-07 

.715 

.8173 

3.4000 

1.420 

-.0195 

8.4742F-07 

6.4439F-07 

2.0303E-07 

.760 

.807! 

3.4500 

1.410 

-.0132 

8.5267E-07 

6.9559E-07 

1.5707E-07 

.816 

.7953 

3.5000 

1.400 

-.0094 

8.5679E-07 

7 • 350 1 F-07 

1.2179E-07 

.858 

.7890 

3.6000 

1.385 

-.0052 

8.6296E-07 

7.B958E-07 

7.3376E-00 

.915 

.7831 

3.7000 

1.374 

-.0036 

3.6654E-07 

9.1424E-07 

5.2275E-09 

.940 

.7839 

3.4000 

1.364 

-.0034 

8.6B24E-07 

8.2006E-07 

4.8183E-08 

.9*5 

.7894 

3.9000 

1.157 

-.0038 

8.6916E-07 

8.1737E-07 

5.1 788E-08 

.940 

.7951 

4. 0000 

1.351 

-.0044 

8.6903E-07 

8.0929E-07 

5.9734E-08 

.911 

.0013 

4.1000 

1 . 346 

-.0056 

H.6783E-07 

7.9861E-07 

6.9222E-09 

.920 

.8073 

4.2000 

1.342 

- . 0069 

8.6579E-07 

7.8574E-07 

8.0048E-08 

.900 

.8131 

4.3000 

1.338 

-.0085 

8.6294F-07 

7.7035E-07 

9.2587E-08 

.893 

.8192 

4. 4000 

1.334 

-.0103 

8.5919E-07 

7.S31SE-07 

1 .0604E-07 

.877 

.8254 

4.5000 

1.332 

-.0134 

6.5442E-07 

7.2734F-07 

1.2708E-07 

.851 

.832? 

4.4000 

1.330 

-.0147 

8.5096E-07 

7.1739E-07 

1 , 3363E-0  7 

.043 

.835? 

4.^000 

1.330 

-.0157 

8.4840F-O? 

7.1047E-07 

1 >3800E-07 

.037 

.0364 

4.4000 

1.330 

-.0150 

8.4743E-07 

7. 1580 E -07 

1 .3155E-07 

.045 

.8349 

4.Q000 

1 . 328 

-.0137 

S.4557E-07 

7.2459E-07 

1 • 2098E-0  7 

.057 

,8336 

5. 0000 

1.325 

-.0124 

8.4274F-07 

7.3258E-07 

1 .1016E-07 

.869 

.0331 

5.1000 

1.322 

-.01 1 1 

8.3951E-07 

7.4031 F-07 

9.9199E-08 

.082 

.8325 

5.?000 

1.317 

-.0101 

8.3400E-07 

7.43B9E-07 

9.0108E-08 

.892 

.0338 

5.7000 

1.312 

-.0098 

8.2706E-07 

7.4095F-07 

8.61 16E-08 

.096 

.8363 

5.4000 

1.305 

-.0103 

8,  1 677E-07 

7.2909E-07 

8.7683E-08 

.893 

.8413 

5.5000 

1.298 

-.0116 

8.0465F-07 

7. 10255-07 

9.4402E-08 

.883 

.8*73 

5.4000 

1.289 

-.0142 

7.8796E-07 

6. 79855-07 

1.081  IE-07 

.863 

,8558 

5.7000 

1.277 

'••0203 

7.6288E-07 

6.2538F-07 

1.3750E-07 

.820 

,8685 

5.5000 

1.262 

-.0330 

7.2692E-07 

5.4195E-07 

1 .8497E-07 

.746 

.0845 

5.5000 

1.248 

-.0622 

6.8467E-07 

4.3026E-07 

2.5441E-07 

.628 

.9002 

6.8000 

1.26S 

-.1070 

6.8877E-07 

3.7273E-07 

3.1604E-07 

.541 

.9004 

- I • 0 


Table  A- 2 


WATF»  4. 


WAVLGH 

refr-iwices 

EXTNC-CDEF 

SCAT  T-COEF 

A8SRP-C0EF 

ALRDO 

6 

2. 6000 

1 

.242  -.0032 

3.0467F -06 

2.66&5E-06 

3.8616E-07 

.873 

.8780 

2.6500 

1 

.219  -.0067 

3.0668F-06 

2.4278E-06 

6.3902E-07 

.792 

.9034 

2.7000 

1 

.188  -.0190 

3.0677E-06 

2.0058E-06 

1 .0619E-06 

.654 

.9426 

2.7500 

1 

.157  -.0590 

2.9293E-06 

1.567?E-06 

1 • 3621E-06 

.535 

.9716 

2.0000 

1 

.142  -.1150 

2 , 7926E-06 

1 .3Bl*E-06 

1.4109E-06 

.495 

.9748 

2. 0500 

1 

.149  -.1850 

2.7746E-06 

1.3608E-06 

1.4138E-06 

.490 

.9634 

2.4000 

1 

.201  -.2680 

2.8245E-06 

1.4137E-06 

1 .41 08E-06 

.501 

.9558 

2. 4500 

1 

.292  -.2980 

2 . 8769E-06 

1.4595E-06 

1 * 4 1 74E-06 

.50? 

.9462 

3.O000 

1 

.371  -.2720 

2.9132006 

1.484  J E-06 

1 .4292E-06 

.509 

.9419 

3.0500 

1 

.426  -.2400 

2.9369E-06 

1 .499QE-06 

1 .4370E-0S 

.511 

.9387 

3.1000 

1 

.467  -.1920 

2 • 9593006 

1.5136E-06 

1 .4447E-06 

.512 

.936? 

3.1500 

1 

.483  -.1350 

2.9786006 

1.5302E-O6 

1.4484E-06 

.514 

,9336 

3. »000 

1 

.478  -.0924 

2.9957006 

1 .5587E-06 

1 .4371E-06 

.520 

.928? 

3.  ?50O 

1 

.467  -.0610 

3.01 16E-06 

1 .61 73E-06 

1 .3943E-06 

.537 

.9165 

3.7000 

1 

.450  -.0368 

3.0282006 

1 .7437^-06 

1 • 2849E-06 

.576 

.8950 

3.7500 

1 

.432  -.0261 

3.0415E-06 

1 .8708E-06 

1.1 706E-06 

.615 

.8783 

3.4000 

1 

.420  -.oic 5 

3. 0532006 

2*001  IE-06 

1 .0522E-06 

.655 

.8633 

3.4500 

1 

.410  -.0132 

3.0663E-06 

2. 19S3r-06 

8 . 7846F-0  7 

.714 

.8448 

3.400,0 

1 

.400  -.0094 

3.0755E-06 

2.3528F-06 

7.2275E-07 

.765 

.8316 

3,  #,000 

1 

.3e5  -.0052 

3.0954E-06 

2.51  96006 

4.7584E-07 

.846 

.8142 

3.7000 

1 

.174  -.0036 

3.1 143006 

2.759SE-06 

3 , 5448E-0  7 

,886 

• 807l' 

3.4000 

1 

.364  -.0034 

3,1 318E-06 

2.8013E-06 

3 ■ 3044E-0  7 

.894 

.809? 

3.4000 

1 

.357  -.0038 

3.1 50 1 006 

2.7959F-06 

3.5423E-07 

.888 

.8126 

4,0000 

1 

.351  -.0046 

3.166&OQ6 

2.7633E-06 

4 • 0330E-0  7 

.873 

.8178 

4.1000 

l 

.146  -.0056 

3.1  825E-06 

2.7218E-06 

4.6074E-07 

.855 

.8234 

4,3000 

1 

.3*2  -.0069 

3.1 968E-06 

2.6720E-06 

5 . 2483E-0  7 

.836 

.8292 

4.7000 

l 

.338  -.0085 

3.2090E-06 

2.6137E-06 

5.9531 F-0  7 

.8]4 

.8358 

4. <*000 

1 

.334  -.0103 

3.2199E-06 

2. 5522E-06 

6.6773E-07 

,793 

.8429 

4.4000 

l 

.332  -.0134 

3.2264E-06 

2.4528F-06 

7.7359E-07 

.760 

.8520 

4, #,000 

1 

.330  • 0 1 47 

3.2356E-06 

2.4294E-06 

8 • 061 4E-0  7 

.751 

.8552 

4.7000 

1 

.330  -.0157 

3.2442E-06 

2.4164E-06 

8.2781F-07 

.745 

.8567 

4.7000 

1 

.330  -.0150 

3.2562E-06 

2.457>E-06 

7.9903E-07 

.755 

..  8540 

4.Q000 

1 

.326  -.0137 

3,270  IE-06 

2.5204E-06 

7.4962E-07 

.771 

.8507 

5.O000 

l 

•325  -.0124 

3 ,284000  6 

2.5871E-06 

6.9683E-07 

.788 

• 8480 

5.1000 

1 

.322  -.0111 

3.2974E-06 

2.5566E-06 

6.4083E-07 

.806 

.8453 

5,3000 

1 

.317  -.0101 

3.3105F-06 

2.71  75E-06 

5.9299E-07 

.821 

.8446 

5.3000 

1 

.312  -.0098 

3. 321 2E-0b 

2.7490E-06 

5.7221E-07 

.828 

.8461 

5,4000 

1 

.305  -.0103 

3,3293006 

2.7465006 

5.82T6E-07 

.825 

.8509 

5.4000 

1 

•298  -.0116 

3.3331F-06 

2.71  UE-06 

6.2168E-07 

,813 

.0577 

5.4000 

l 

.289  -.014? 

3. 3299006 

2.6328E-06 

6.971  IE-07 

.791 

.8681 

5,7000 

l 

. ?77  -.0203 

3. 3083F-06 

2.4629E-a6 

8.4544E-07 

.744 

.8850 

5.4000 

1 

.262  -.0330 

3.2504E-06 

2. 1977E-06 

1 .0527E-06 

.676 

.9071 

5.4000 

1 

.248  -.0622 

3.1268F-06 

1.9374E-06 

1 .2894E-06 

.588 

.9297 

6.6000 

1 

.265  -.1070 

3.0520E-06 

1.61 1 8E-06 

1 .4402E-06 

.528 

.9355 

-1.0 
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Table  A- 3 


MATES  8. 


WAvLGH 

refr-inoices 

EXT^C-COEE 

SCATT-COEF 

ABSRP-COEF 

ALBOO 

6 

2.6000 

1.242 

-.0032 

1.1 357E-05 

8. 

9470E-06 

2.41 00E-06 

.788 

.9072 

2. 6500 

1.219 

-.0067 

1 . 1 373E-05 

7. 

7$5r>E-06 

3.5785E-06 

.685 

.9304 

2. 7000 

1.188 

-.0190 

1.1 408E-05 

6. 

4456E~06 

4.9622E-06 

.565 

.9626 

2.7500 

1.157 

-.0590 

1.1 28  IE-05 

5. 

7621 F -06 

5.51 88E-06 

.511 

.9814 

2.8000 

1.142 

-.1150 

1.1 056E-05 

5. 

5550E-06 

5.501  IE-06 

.502 

.9817 

?.«500 

1.149 

-.1850 

1.1001E-05 

5. 

5835E-06 

5.41 70E-06 

.508 

.9750 

2.QOOO 

1.201 

-.2680 

1.1 079E-05 

5. 

7481E-06 

5.3314E-06 

.519 

.9627 

2.O500 

1.292 

-.2980 

1.1 1 72E-05 

5. 

8649E-06 

5.3071E-06 

.525 

.95^5 

3.0000 

1.371 

-.2720 

1.1 235F-05 

5. 

9135E-06 

5. 321 9E-06 

.526 

.9497 

3=0500 

1.426 

-.2400 

1.1 277E-05 

5. 

9446E-06 

5.3325E-06 

.527 

.9471 

3.1000 

1.467 

-.1920 

1.131 4E-05 

5. 

963?E-06 

5.3503E-06 

.527 

,9460 

3.1500 

1.483 

-.1350 

1 .1346E-05 

5. 

969?E-06 

5. 3773E-06 

.526 

.9465 

3. ?000 

1.478 

-.0924 

1.1 373E-05 

5. 

9778E-06 

5. 3957E-06 

.526 

.9472 

3.S500 

1.467 

-.0610 

1.1 398E-05 

6. 

020SE-06 

5.3775E-06 

.528 

.9459 

3.1000 

1.450 

-.0368 

1.1 422E-05 

6. 

1725F-06 

5.2496E-06 

.540 

.9394 

3.1500 

1.432 

-.0261 

1 .1443F-05 

6. 

3B61E-06 

5.0570E-06 

.558 

.9314 

3.4000 

1.420 

-.0195 

1 .1455E-05 

6. 

6565F-06 

4.7990E-06 

.581 

.9215 

3.4500 

1.410 

-.0132 

1 .1483E-05 

7. 

156OE-06 

4.3266E-06 

.623 

.9054 

3.0000 

1.400 

-.0094 

1.1493E-05 

7. 

6867E-06 

3.8067E-06 

.669 

.8911 

3.6000 

1.385 

-.0052 

1.1538F.-05 

8. 

7459F-06 

2.7921E-06 

.758 

.8675 

3.7000 

1.374 

-.0036 

1.1576E-05 

9. 

374SE-06 

2.201 9E-06 

.810 

.8563 

1.4000 

1 » 364 

-.0034 

1.161 3E-05 

9. 

5271 E -06 

2.0859E-06 

• 8?0 

.8550 

3. sGGo 

1.357 

-.0038 

1 • 1 626E-05 

9. 

41 9nE”06 

2 . 2070E-06 

.810 

.8586 

4.0000 

1.351 

-.0046 

1 . 1631E-05 

9. 

1797E-06 

2.451 8E-06 

.789 

.8643 

4.1 000 

1.346 

-.0056 

1 • 1651E-05 

8. 

91665-06 

2.7344E-06 

.765 

.8698 

4,?000 

1.342 

-.0069 

1.1680r>05 

8. 

6415E-06 

3.03R3E-06 

.740 

.8754 

4.1000 

1.138 

-i 0085 

1 .1710E-05 

8. 

3557F-06 

3.3545E-06 

.714 

.8818 

4.4000 

1.334 

-.0103 

1 .1740E-05 

B • 

0831 5-06 

3.6573E-06 

.688 

,8885 

4.4000 

1.332 

-.0134 

1 .1769E-05 

7. 

7055E-06 

4. 0630E-06 

.655 

.8979 

4 , AO  00 

1.130 

-.0147 

1.1 797E-05 

7. 

6129E-06 

4.1 840E-06 

. 645 

.9005 

4. 7000 

1.330 

-.0157 

1 . 1 820E-05 

7. 

5600E-06 

4. 2603E-06 

.640 

.9018 

4.4000 

1.130 

-.0150 

1.1845E-05 

7. 

69255-06 

4.1622E-06 

.649 

.8983 

4.Q00O 

1.328 

-.0137 

1.1  874E-05 

7. 

985?E-06 

3.9890F-06 

.664 

.8933 

S.nOOO 

1.325 

-.0124 

1 . 1 905E-05 

8. 

110?E-06 

3.7951E-06 

.681 

.8884 

5.1000 

1.322 

-.0111 

1.1939E-05 

8. 

36095-06 

3.5784E-06 

.700 

.8833 

5.7000 

1.317 

-.0101 

1 . 1 976E-05 

8. 

590?f-06 

3.3854E-06 

.71  7 

.8796 

5.7000 

1.312 

-.0098 

1 .2014E-05 

9. 

71 0SE-06 

3.3035E-06 

.725 

.8786 

5.4000 

1.105 

-.0103 

1.2053E-05 

8. 

69505-06 

3.3577E-06 

.721 

.881? 

5.5000 

1.298 

-.0116 

1 .2092F-05 

9. 

55835-06 

3.5333E-06 

.708 

.8868 

5.4000 

1.289 

-.0142 

1 .2130E-05 

B • 

27Bi c-06 

3.8516E-06 

.682 

.8958 

5.7000 

1.277 

-.0203 

1 .2158E-05 

7. 

7474E-06 

4.4104E-06 

.637 

.912? 

5.4000 

1.262 

-.0330 

1 .2145E-05 

7. 

088] 5-06 

5.0574E-06 

.584 

.9335 

5.Q000 

1.248 

-.0622 

1 .2003E-05 

6. 

4062E-06 

5.5971E-06 

.534 

.9538 

ft.nooo 

1.265 

-.1070 

1.1  81 5E-05 

6. 

0201 r-06 

5.7953E-06 

.510 

.9585 

1.0 
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Table  A-4 


WATEP 

16. 

wavlgh 

9EFR-INDICES 

2.5000 

1.242  -.0032 

2.5500 

1.219  -.0067 

2.7000 

1.188  -.0190 

2.7500 

1.157  -.0590 

2.0000 

1*142  -.1150 

2.8500 

1.149  -.1850 

2.0000 

1.201  -.2580 

2. 0500 

1.292  -.2980 

3.0000 

1.371  -.2720 

3. 0500 

1.426  -.2400 

3.1 000 

1.467  -.1920 

3.1500 

1.483  -.1350 

3.?000 

1 -478  -.0924 

3 . ?500 

< .467  -.0610 

3.1000 

1 . *.50  - .0368 

3.1500 

1.432  -.0261 

3.4000 

1.420  -.0195 

3.4500 

1.410  -.0132 

3.0000 

1.400  -.0094 

3.5000 

1.385  -.0052 

3.7000 

1.374  -.0036 

30000 

1.364  -.0034 

3.Q00O 

1.357  -.0038 

4.0000 

1.35V  -.0046 

4.1 000 

1.346  -.0056 

4. ?000 

1.342  -.0069 

4.1000 

1.338  -.0085 

4.4000 

1.334  -.0103 

4. *000 

1.332  -.0134 

4.5000 

1.330  -.0147 

4.7000 

1.330  -.0157 

4.8000 

1.330  -.0150 

4.0000 

1.128  -.0137 

5.0000 

1.325  -.0124 

5.1000 

1.322  -.0111 

5.POOO 

1.317  -.0101 

5.1000 

1.312  -.0098 

5.4000 

1.305  -.0103 

5. soon 

1.298  -.0116 

5.5000 

l.?89  -.0142 

5.7000 

1 .277  -.0203 

5.4000 

1.262  -.0330 

5,0000 

1 .248  -.062? 

6. n000 

1.265  -.1070 

-1.0 

EXTMC-CDEF  SCfcTT-COEF 
4.3990E-05  3.0324E-C5 
4.3947E-05  2.629?E-0S 
4.3943E-05  2.3141E-05 
4.3747E-05  2.2269E-05 
4.3432E-05  2.2l3oE-05 
4.3356E-05  2.2405E-O5 
4.3511L-05  2.3055F-05 
4.3691E-05  2.3413E-05 
4.3815E-05  2.3541E-05 
4.3899E-05  2.3626E-05 
4.3970E-05  2.3665E-05 
4.4034E-05  2.3647005 
4.4087F-05  2.3606E-05 
4.4136E-05  2.35B7E-05 
4.41R4E-05  2.366OE-05 
4.4220E-O5  2.3B5OE-05 
4.4252E-05  2.42UE-05 
4.4289E-05  2.510RE-05 
4.4327E-05  2.S343F-05 
4.4397E-05  2.9560E-05 
4.4491E-05  3.1973E-05 
4.4524F-05  3.2514E-05 
4.459‘»F“C'"  3.2037F-05 
4.4580F-05  3.097q£-05 
4.4S74E-05  2,99595-05 
4.4679E-05  2.99S6E-05 
4.4784E-05  2.7948E-05 
4.4786E-05  2.7043E-05 
4.4857E-05  2.502(SF-0S 
4.4898F-05  2.574QE-05 
4.4930F-05  2.5604E-05 
4. 50165-05  2.58865-05 
4.50 1 6F-05  2.6329E-05 
4.5133E-05  2.5927E-05 
4.5125F-05  2.7545F-05 
4.5250E-05  2.9233E-05 
4 . 5244E -0 5 2.9493E-05 
4.5364E-05  2.339*5-05 
4.5360E-05  2.79285-05 
4.5476F-05  2.S96PE-05 
4.5475E-05  2.553JF-05 
4.551  IE-05  2.41  975-05 
4.5375E-05  2.323*5-05 
4,51 45E-05  2. 30395-05. 


A3S8P-C0EF 

ALB00 

G 

1 * 3S65E-05 

.691 

.9381 

1.7655E-05 

.598 

.9591 

2.0802E-05 

.527 

.9792 

2.  J47BE-05 

.509 

.9863 

2.1252E-05 

.51  1 

.9846 

2. 087  IE-05 

.519 

.9778 

2. 0456E-05 

.530 

.9655 

2.0278E-05 

.536 

.9562 

2.0274E-05 

.537 

.9526 

2.0272E-05 

.538 

.9501 

2.0305E-05 

.538 

.9491 

2.0388E-05 

.537 

.9501 

2.0481E-05 

.535 

.9519 

2.0549E-05 

.534 

.9535 

2.0524E-05 

.535 

.9541 

2.0371E-05 

.539 

.9533 

2.0038E-05 

.547 

.9506 

1.91812-05 

.567 

.9430 

1.7985E-05 

.594 

.9334 

1.4837E-05 

.666 

.9127 

1.2518E-05 

.719 

.8994 

1.2006E-05 

.730 

.8986 

1.2557E-05 

.718 

.9016 

1 . 3602E-05 

.695 

.9090 

1.4715E-05 

.671 

.9147 

1.5793E-05 

.64  7 

.9217 

1 .6836E-05 

.824 

.9280 

1 .7743E-05 

.604 

.9343 

1. 88315-05 

.580 

.9418 

1.9147E-05 

.574 

.9437 

1 .9326E-05 

.570 

.9449 

1.91 30E-05 

.575 

.942? 

1.8687E-05 

.585 

.9390 

1.8206E-05 

.59  7 

.9347 

1.7579E-05 

.610 

.9305 

1.701 7E-05 

.624 

.9266 

1.6762E-05 

.630 

.9252 

1.696BE-05 

i 626 

.9268 

1.7532E-05 

.613 

.9310 

1 .8506E-05 

.593 

.9378 

1.9945E-05 

.561 

.9495 

2.1315E-05 

.532 

.9619 

2.2090E-05 

.51  3 

.9710 

2.2106E-05 

.510 

.970  3 
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Table  A- 5 


WATFO  32. 


WAVLOH 

REFR-INOICES 

2. <,000 

1 

.242 

-.0032 

2,*500 

1 

.219 

-.0067 

2.7000 

1 

.188 

-.0190 

2.7500 

1 

.157 

-.0590 

2.8000 

1 

.142 

-.1150 

2.8500 

1 

• 149 

-.1850 

2.Q000 

1 

.201 

-.2680 

2.9500 

1 

.292 

-.2980 

3.0000 

1 

.371 

-.2720 

3. 0500 

1 

.426 

-.2400 

3.1 000 

1 

• 467 

-.1920 

3.1500 

1 

.483 

-.1350 

3.2000 

1 

.478 

-.0924 

3.^500 

1 

.467 

-.0610 

3.1000 

1 

.450 

-.0368 

3.3500 

1 

.432 

-.0261 

3.4000 

1 

.420 

-.0195 

3.4500 

1 

.410 

-.0132 

3.8000 

1 

.400 

-.0094 

3.8000 

1 

.385 

-.0052 

3.7000 

1 

.374 

-.0036 

3.4000 

1 

.364 

-.0034 

3.9000 

1 

.357 

-.0038 

4 . nOOO 

1 

.351 

-.0046 

4.1  000 

1 

.346 

-.0056 

4.2000 

1 

.342 

-.0069 

4.3000 

1 

.338 

-.0085 

4.4000 

1 

.334 

-.0)03 

4.8000 

1 

.332 

-.0134 

4.8000 

l 

.330 

-.0147 

4.7000 

1 

.330 

-.0157 

4.POOO 

1 

.330 

-.0150 

4.9000 

1 

.328 

-.0137 

S.nOOO 

1 

.325 

-.0124 

5.1000 

1 

.322 

-.0111 

5.7000 

1 

.317 

-.0101 

5.3000 

1 

.312 

-.0098 

5.4000 

1 

.305 

-.0103 

5.5000 

1 

.298 

-.0116 

5.8000 

1 

• ?89 

-.0142 

5.7000 

1 

.277 

-.0203 

5,4000 

1 

.262 

-.0330 

5.9000 

1 

.248 

-.0622 

6.n000 

1 

.265 

-.1070 

-1.0 


EX7MC-C3EF  SCAT  T-COEF 
1.7059E-04  1.034?E-04 
1.7067F-04  9.344QE-05 
1.7065E-04  B.9546E-05 
1.7026F-04  9.7516E-0S 
1.6975E-04  8.7754E-05 
1.6967E-04  8.906?E-05 
1.7000E-04  9.1066E-05 
1 « 7036E-04  9. 2269E-05 
1.7061F-04  9.2661E-05 
1.7079E-04  9.2924E-05 
1.7094F-04  9.3023E-05 
1.7107E-04  9.2907E-05 
1 • 7 1 1 8F-04  9.269QE-05 
1.7129F-04  9. 25Q6E-05 
1.7139F-04  9. 23?6E“05 
1.7147F-04  9.2253F-05 
1.71S5F-04  9.244(SE“05 
1.7163E-04  9.3374F-05 
1.7171F-04  9.5149E*05 
1.7183F-04  1.0189F-04 
1.7201F-04  1.0872F-04 
1.7208F-04  1.1041E-Q4 
1.7225E-04  l’.0873E“04 
1.7236E-04  1.0559E-04 
1.7248E-04  1.025OE-04 
1.7259E-04  9. 997RE-05 
1.7272F-04  9.770*F-05 
1.7281F-04  9.5923F-05 
1.7294E-04  9.407SF-05 
1.7303E-04  9.364OE-05 
1.7315E-04  9.3431E-05 
1.7328F-04  9.3B44E-05 
1.7338E-04  9.4575E-05 
1.7352E-04  9.554OF-05 
1 • 7365E-04  9.676«E-05 
1.7375E-04  9.B013E-05 
1.7390F-04  9.9631E-05 
1.7400E-04  9.920BE-05 
1.7411E-04  9.704SE-05 
1.7422F-04  9.5177E-05 
1.7429E-04  9.2715E-05 
1.7428E-04  9.0946E-05 
K7405F-04  8.9905E-05 
1.7376F-04  9.020AE-05 


ABSRP-COEF 

ALRDO 

G 

6.71 73E-05 

.606 

.9606 

7.7218E-05 

.548 

.9757 

8.2104E-05 

.519 

• 98£*1 

8.?748E-05 

.514 

.9878 

8.2000E-05 

.517 

.9857 

8.0607E-05 

.525 

.9789 

7.B929E-05 

.536 

.9665 

7.R0O2E-05 

.542 

.9573 

7.7953E-05 

.543 

.9536 

7.7865E-05 

.544 

.9511 

7.79) 7E-05 

.544 

.9501 

7.R166E-05 

.543 

.951? 

7.8485E-05 

.542 

.953? 

7.8781E-05 

.540 

.9552 

7.9062F-05 

.539 

.9573 

7.9218E-05 

.5  38 

.9590 

7.9105E-05 

.539 

.9595 

7.R257E-05 

.544 

.9581 

7.6557E-05 

.554 

.9549 

6 . 9940E-05 

.593 

.9429 

6.3288E-05 

.632 

.9325 

6.1669E-05 

.642 

.9313 

6.3515E-05 

.631 

.9347 

6.6772E-05 

.613 

.9406 

6.9890E-05 

.595 

.9461 

7.P60BE-05 

.579 

.951? 

7,501 7E-05 

.566 

.9556 

7.6891E-05 

.555 

.9594 

7.8861E-05 

.544 

.9632 

7.9393E-05 

.541 

.9643 

7.9718E-05 

.540 

.9647 

7.9434E-05 

.542 

.9638 

7.8807E-05 

.545 

.9624 

7 , 7980E-05 

.551 

.9607 

7,688  IE-05 

.557 

.9585 

7.5734E-05 

.564 

.9567 

7.5266E-05 

.56? 

.9560 

7.5788E-05 

.564 

.9574 

7.7069E-05 

.557 

.960? 

7.9042E-C5 

• 546 

.9643 

8, 1563E-05 

.532 

.9701 

8.3435E-05 

.521 

.9750 

8.4145E-05 

.517 

.9774 

8. 3554E-05 

.519 

.9743 
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Table  A-6 


ICE  ?. 


XAvLGM 

REFR-I90IC6S 

EXT9C-C3EF 

SC0TT-COE6 

ABSRP-C06F 

albdo 

G ! 

2.5000 

225  -.0000 

8..2926F-07 

7.59836-07 

6.90356-08 

.916 

.8865 

2.5600 

212  -.0000 

8,26696-07 

7.59566-07 

6.70336*08 

.919 

.8960 

2.6320 

192  -.0050 

8 , 1 51 3E-07 

7.37BoE-07 

7.7298E-08 

.905 

.9123 

2.7030 

162  -.0080 

7 • 7792F-07 

6.71795-07 

1.0613E-07 

.860 

.9358 

2.7780 

108  -.0200 

6.28806-07 

O.080O5-07 

1 .BuBOE-O  7 

.712 

.9667 

2.«170 

065  -.0380 

O.6772E-07 

2.33036-07 

2.3069E-07 

.098 

.9776  1 

2.0250 

056  -.0000 

0.0553E-07 

1.99376-07 

2.O616E-07 

.007 

.9780  1 

2.0330 

006  -.0520 

0.32116-07 

1 .73So6-07 

2.5857E-07 

,002 

.9785 

2.0010 

037  -.0600 

O.2951E-07 

1.612(16-07 

2.6831E-07 

.375 

.9761 

2.0090 

027  -.0680 

0.3130E-07 

1.5577F-07 

2 . 7553E-0  7 

.361 

.9773 

2.0570 

016  -.0790 

0.0051E-07 

1.60686-07 

2.8383E-07 

.361 

.9757 

2. 0700 

995  -.1100 

0.93896-07 

1,93366-07 

3.00S3E-07 

.391 

.9698 

2.0900 

979  -.1060 

5.0006E-07 

2.27516-07 

3, 1255E-07 

.021 

.9623  I 

2.9070 

960  -.1930 

5.82156-07 

2.61066-07 

3.21096-07 

.008 

.9522 

2.Q2O0 

967  -.2570 

6 . 2277F-0  7 

2.99396-07 

3.3338E-07 

.065 

<■9812 

2.0330 

978  -.2760 

6.333OE-07 

2.35656-07 

3.3769E-07 

.067 

.9386 

2.0010 

987  -.2950 

6.O301F-07 

3.01976-07 

3.01 0OE-07 

.070 

.9358 

2.0590 

000  -.3310 

6.59O3E-07 

3,132(16-07 

3.O623E-07 

.075 

.9302 

2.0760 

021  -.3670 

6.7375F-07 

3. 23506-07 

3.S021E-07 

.080 

.9200  > 

2.0000 

035  -.0110 

6.8787F-07 

3.350?6-07 

3.5285E-07 

.087 

.9169 

3.0120 

060  -.0660 

7 . 0036F-0  7 

3,08056-07 

3. 56316-07 

.090 

.9080 

3.0300 

100  -.5220 

7 . 1 975F-0  7 

3.50186-07 

3.5957E-07 

.500 

.8991  i 

3.0090 

170  -.5800 

7 . 3635E-0  7 

3.72966-07 

3.63396-07 

.507 

.8890 

3. 0580 

217  -.6120 

7.O056E -07 

3.78976-07 

3.  #>5596-07 

.509 

.8809  1 

3.0670 

273  -.6270 

7.S162E-07 

3.83386-07 

3.68206-07 

.510 

.8821  ) 

3.0770 

330  -.6250 

7.570  IE-07 

3.359?c-07 

3.7109E-07 

.510 

.881?  ; 

3.0860 

380  -.6100 

7.6107F-07 

3.8778F-07 

3.7362E-07 

.509 

.8809  j 

3.0960 

032  -.5930 

7 .6502E-07 

3.98836-07 

3.76196-07 

.500 

.8815  i 

3.1060 

069  -.5600 

7.6760E-07 

3.990OF-07 

3.78576-07 

.507 

.8031  i 

3.1150 

099  -.5390 

7.6995F-07 

3.B95OF-07 

3.8002E-07 

.506 

.8800 

3.1250 

527  -,'150 

7.72366-07 

3.903(16-07 

3.R206E-07 

• 5u5 

.8605  j 

3.1350 

553  -.0910 

7 . 70766-07 

3.91216-07 

3.8355E-07 

.505 

.8808  j 

3.1050 

576  -.0650 

7.772OF-07 

3.92236-07 

3.8500E-07 

.505 

.8850  ! 

3.1550 

603  -.0380 

7 ,79866-07 

3.935(16-07 

3.R536E-07 

.505 

.8850  ' 

3.1650 

623  -.0060 

7.82O3F-07 

3.905nc--o7 

1.8783E-07 

.500 

.885  3 

3.1750 

600  -.3720 

7 . 851 QF-07 

3.95866-07 

3.8920E-07 

.500 

.8855  , 

3.1850 

650  -.3360 

7.8775E-07 

3,97126-07 

3.9053E-07 

.500 

.8860  < 

3. ?050 

653  -.2670 

7.93UF-07 

0.00036-07 

3. 92696-07 

.505 

,8865  1 

3.7?60 

601  -.2100 

7.98086-07 

0.05566-07 

3.9292E-07 

.508 

.8852 

3.5790 

593  -.1210 

8.10976-07 

0.28916-07 

3.8216E-07 

.529 

.8725  ; 

3.73^0 

505  -.0700 

8 . 2238F-0  7 

0.5B52F-07 

3.5386E-07 

.570 

.8509  ; 

3.3900 

511  -.0010 

8.3013F-07 

5.3798F-07 

2.95206-07 

.605 

.8173  ! 

3. <.080 

076  -.0200 

8.O392F-07 

6.13756-07 

2.30176-07 

.727 

.7902  j 

3,5090 

009  -.0150 

8.5157F-07 

6.79536-07 

1.7200E-07 

.798 

.7810  j 

3.57.0 

028  -.0090 

8. 57906-07 

7.39816-07 

1.18106-07 

.862 

.7728  i 

3.5360 

010  -.0080 

«. 61526-07 

7«559i 6*07 

1 .05616-07 

.877 

. T7B3  1 

3.7000 

398  -.0070 

8.60596-07 

7. 713?6-07 

9,32706-00 

.892 

.7811  1 

3.S060 

376  -.0070 

8.67506-07 

7.78156-07 

8.93O86-(<0 

.897 

.7931 

O.nOOO 

361  -.0100 

8.65976-07 

7 . 52396-07 

1.13586-07 

.869 

.8090  j 

0.1670 

308  -.0120 

8.62786-07 

7.37086-07 

1 » 2530F-0  7 

.855 

.8210 

0.5550 

303  -.0150 

8.58696-07 

7.10206-07 

1.00506-07 

.832 

.8295  i 

0,7080 

337  -.0200 

8.51686-07 

6.7905E-07 

1.72236-07 

.798 

.8009 

0.  i860 

337  -.0230 

8.OB73E-07 

6.51626-07 

1.871 16-07 

.780 

.0050 

0.4000 

338  -.0260 

0. 06056-07 

6.05986-07 

2.00066-07 

.760 

.8081  j 

0.5050 

300  -.0290 

8.03796-07 

6.32026-07 

2.11766-07 

.709 

.8500  ( 

0.5050 

300  -.0300 

8. 00376-07 

6.2899‘„-07 

2.15386-07 

.705 

.809?  j 

0.6080 

307  -.0250 

8.0939E-07 

6.5678F-07 

1.92616-07 

.773 

.8008 

0.6510 

307  -.0210 

0.52096-07 

6.80796-07 

1.71696-07 

.799 

.8309  ) 

0.766? 

301  -.0150 

8,53786-07 

7.20376-07 

1.33016-07 

.800 

.8280  j 

5.0000 

327  -.01-20 

8.O0O06rI07 

7.  36856-07 

1 .1)7606-07 

.873 

.8311 

5.5630 

312  -.0100 

8.23956-07 

7,091 2r~07 

1.1080E-07 

.861 

.8027 

5.5560 

296  -.0200 

7.88176-07 

6.28076-07 

1.60106-07 

.797 

.8620 

5.7100 

287  -.0350 

7.61956-07 

5.60106-07 

1 .97866-07 

.700 

.8701 

5.7070 

285  -.0380 

7.55806-07 

5.09336-07 

2.06076-07 

.727 

.8767  ] 

5.7800 

280  -.0020 

7.501 9F-07 

5,32726-07 

2.17076-07 

.710 

.8791 

5.B1O0 

283  -.0060 

7.00596-07 

5.17266-07 

2.27036-07 

.695 

.8813 

5.8080 

283  -.0520 

7.39286-07 

0.97886-07 

2.O1O0E-07 

.673 

.8837 

5.8820 

285  -.0570 

7 . 37 1 9E-0  7 

0,95106-07 

2.52056-07 

.658 

.8805 

6.0610 

301  -.0690 

7.01796-07 

0.67596-07 

2.70206-07 

.630 

.8808 

1 


j 


•; 


■4 

■j 

i 
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Table  A-7 


t TT  *W . ,•  > '.v^r  ■■  -y  -V 


ICE 


WAvLOH 

9EFH- 

M01CES 

Exm-coer 

SCATT-COEF 

absrp-coff 

ALRDO 

G 

2.5000 

1.225 

-.0040 

3.040  3F -06 

2.5  r’54E-06 

4.67R9E-07 

.846 

.8914 

2.5640 

1.212 

-.0040 

3.0540F-06 

2.6086E-06 

4.5533E-07 

.851 

.8980 

2,6320 

1.192 

-.0050 

3.0961E-06 

2.5768E-06 

5, 1 928E-0  7 

.832 

.9124 

2.7030 

1.162 

-.0080 

3.1 206E-06 

2.4418F-06 

6.78735-07 

.782 

.9376 

2.7780 

1 • 1 OB 

-.0200 

2.9771F-06 

1.9667^-06 

1.0105E-06 

.661 

.9742 

2.8170 

1.065 

-.0380 

2.4870E-06 

1.3124E-06 

1 • 1 746E-06 

.528 

.9868 

2.8250 

1.056 

-.0440 

2. 370  7F-06 

1.16935-06 

1.201SE-06 

.493 

.9878 

2,8330 

i .046 

-.0520 

2.2741F-06 

1.04665*06 

1 .2274E-06 

.460 

.9882 

2.8010 

1.037 

-.0600 

2.2259F-06 

9. 80B2E-07 

l .24515-06 

.441 

.9881 

2.84,90 

1.027 

-.0680 

2.2024E-06 

9.4743F-07 

1.2550E-06 

.430 

.9877 

2.8570 

1.016 

-.0790 

2.2230E-06 

9. 57605-07 

5.26545-06 

.431 

.9865 

2.8740 

.995 

-.1100 

2.3490F-O6 

1.0656E-06 

1 ,2834E-0b 

.454 

.9817 

2.8900 

. .979 

-.1460 

2. 4631E -06 

1 , 1 70  35-06 

1.2927E-06 

.475 

.9754 

2.9070 

.960 

-.1930 

2.5585E-06 

1.26615-06 

1.2923E-06 

.495 

.9664 

2.9240 

.967 

-.2570 

2.6459E-06 

l . 340  3E-06 

1.3057E-06 

.507 

.9566 

2.9330 

.978 

-.2760 

2.6685E-06 

1.35545*06 

1.31 31E-0b 

.508 

.9543 

2.9410 

.987 

-.2950 

2.688BE-06 

1 .371OF-06 

1.3178E-06 

.510 

.9517 

2.9590 

1.004 

-.3310 

2. 7228E-06 

1 . 399OE-06 

1.3238E-06 

.514 

.9466 

2.9760 

1.021 

-.3670 

2.7S18E-06 

1 .4251F-06 

1 . 3268E-06 

.518 

.9413 

2.9940 

1.035 

-.4110 

2.7799E-06 

1 .45595-06 

1.3246E-06 

.523 

.9342 

3.0120 

1.064 

-.4660 

2.8125E-06 

1 .46B7E-06 

1 .3237E-06 

.529 

.9256 

3.0300 

1.104 

-.5220 

2.8427F-05 

1.51965-06 

1.3231E-06 

.535 

.9170 

3.0490 

1.170 

- « 5840 

2.8752F-06 

1.55125-06 

1.3240E-06 

.539 

.9076 

3.0580 

1.217 

-.6120 

2.8911F-06 

1.56515-06 

1.3260E-06 

.541 

.9032 

3.0670 

1.273 

-.6270 

2.9046F-06 

1.57365-06 

1 ,33lOE-Ob 

.542 

.9006 

3.0770 

1.330 

-.6250 

2.9147F-06 

1 .5763^-06 

1.33B4E-0b 

.541 

.8999 

3.0860 

1 .384 

-.6140 

2.9227E-06 

1.57705-06 

1.3456E-06 

.540 

.8999 

3.0960 

1.432 

-.5930 

2.9287F-06 

1.57545-06 

1.3533E-0b 

.538 

.9008 

3.1060 

1 .*69 

-.5640 

2.9328F-06 

1.57145-06 

1.3614E-06 

.536 

.9026 

3.1150 

1.499 

-.5390 

2.9364F-06 

1 .5690^-06 

1.3674E-06 

.534 

.9038 

3.1250 

1.527 

-.5150 

2 . 940 1 E-06 

1 . 567*5-06 

1.37P6E-06 

.533 

.9046 

3.1350 

1.553 

-.4910 

2.9437E-06 

1 . 56655-06 

1 .3772E-06 

.532 

.9053 

3.1450 

1.578 

-.4650 

2.9474E-06 

1.56575-06 

1.3817E-06 

.581 

.9059 

3.1550 

1.603 

-.4380 

2.951 1F-06 

1.56545-06 

1 . 3857E-06 

.530 

.9064 

3 • 1650 

1.623 

-.4060 

2.9546F-06 

1. 56405-06 

1.3906E-06 

.529 

.9074 

3.1750 

1.640 

-.3720 

2.95817 -06 

1.56265-06 

1.3956E-06 

.528 

.9084 

3.1850 

1.650 

-.3360 

2 • 96 1 4F-0b 

1.56025-06 

1.401 2E -06 

.527 

.9099 

3.9050 

1.653 

-.2670 

2.9680F-06 

1.55525-06 

1.4128E-06 

.524 

.9134 

3,9260 

1.641 

-.2100 

2.9748F-06 

1.55175*06 

1.42315-06 

.522 

.9166 

3.9790 

1.593 

-.1210 

2.9914F-06 

1.55815-06 

1 .43375-06 

.521 

.9192 

3.7330 

1.545 

-.0740 

3.0079F-06 

1 .598oc;-06 

1.4099E-0S 

.531 

.9128 

3.3900 

1.511 

-.0810 

3.0255F-06 

1.71925-06 

1 .3064E-06 

.568 

.8903 

3.4480 

1.476 

-.0240 

3.0431F-06 

1.91025-06 

1 . 1 328E-0b 

.628 

.8637 

3.5090 

1.449 

-.0150 

3.0584F-06 

2.12695-06 

9,31 4BE-07 

.695 

.8413 

3.5710 

1.428 

-.0090 

3.0721F-06 

2.37205-06 

7.0018E-07 

.772 

.8211 

3.6360 

1.410 

-.0080 

3.0894E-06 

2.44585-06 

6.4358E-07 

.792 

.8193 

3.7040 

1.398 

-.0070 

3.1  062F-06 

2.52325-06 

5.R303E-07 

.8)2 

.8163 

3.8460 

1076 

-.0070 

3.1  300E-06 

2.56425-06 

5.6582E-07 

.819 

.8218 

4.0000 

1061 

-.0100 

3.1573E-06 

2.46095-06 

6.9636E-07 

.779 

.8360 

4.1670 

1048 

-.0120 

3.191 2E-06 

2.4236E-06 

7.5761E-07 

.762 

.8454 

4.9550 

1043 

-.0150 

3.1898F-06 

2.34115-06 

8.4B71E-07 

.734 

.8548 

4.7480 

1037 

-.0200 

3,1 963E-06 

2.22595-06 

9.7038E-07 

.696 

.8680 

4.7860 

1037 

-.0230 

3, 1 959F-06 

2.16575-06 

1.0302E-06 

.678 

.8737 

4.8440 

1038 

-.0260 

3.1971E-06 

2.117)5-06 

1 .0800E-06 

.662 

.8782 

4.5050 

1040 

-.0290 

3.  1 978E-06 

2. 07505-06 

1.1 228E-06 

.649 

.8819 

4.5450 

1 044 

-.0300 

3.1991F-06 

2.06375-06 

1.1 354E-06 

.645 

• 8B16 

4.6080 

1047 

-.0250 

3.21 23E-06 

2 • 15935-06 

1 .05305-06 

.672 

.8711 

4.6510 

1047 

-.0210 

3.2228E-06 

2.25215-06 

9, 7077E-07 

.699 

.8625 

4.7662 

1041 

-.0150 

3.2475E-06 

2.44255-06 

B.O510F-07 

.752 

.8495 

5 . ft  0 0 0 

1027 

-.oi-eo 

3.2840F-06 

2.60085-06 

6.R328E-07 

.792 

.8459 

5.2630 

1.312 

-.0140 

3.3066F-06 

2.58155-06 

7.P503E-07 

.781 

.8569 

5.5560 

1.296 

-.0240 

3.2981F-06 

2.35455-06 

9.4366E-07 

.714 

. 8R 1 9 

5.7140 

1.2B7 

-.0350 

3.2669F-06 

2.167SE-06 

1 .0994E-06 

.663 

.8988 

5.7470 

1.285 

-.0380 

3.2571E-06 

2.12545-06 

1.131 7E-06 

.653 

.9024 

5.7800 

1.284 

-.0420 

3.2448E-06 

2.07375-06 

1 . 1 71 2E-06 

.639 

.9062 

5.8140 

1.283 

-.0460 

3 . 2326F-06 

2.027?5-06 

1 .2054E-06 

.627 

.9095 

5.8480 

1.283 

-.0520 

3.2156E-06 

1 . 96495-06 

l .2507E-06 

.611 

.9134 

5.8820 

1.285 

-.0570 

3.2045E-06 

1.921ic-06 

1 .2834E-06 

.600 

.9155 

6,0610 

1001 

-.0690 

3.1961F-06 

1.84885-06 

1 .34765-06 

.078 

.9155 

-1.0 
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Table  A-8 


ICE  8. 


WAV LGH 

RE58-I901CES 

EXTVC-C3EE 

SCAT  T-COEF 

ABSRP-COEF 

AL800 

0 

2.5000 

1.225  -.0040 

1,1 330E-05 

8.5100F-06 

2.8201E-06 

.751 

.9182 

2.5640 

1.212  -.0040 

1 .1 374E-05 

b.5996E-06 

2.7746E-06 

.756 

.9210 

2.6320 

1.192  -.0050 

t • 1 422E-05 

8.3377F-06 

3.0844E-06 

.730 

.9307 

2.7030 

1.162  -.0080 

1.1 505F-05 

7. 7204E-06 

3.7851E-06 

.671 

.9492 

2.7700 

1.108  -.0200 

1 . 1 595P-05 

6.7033E-06 

4.891 6E-06 

.578 

.9801 

2.8170 

1.065  -.0780 

1 . 0904E-05 

5.6834F-06 

5.2209E-06 

.521 

.9913 

2.0250 

1.056  -.0440 

1 .0600E-05 

5.3492F-06 

5.2S07E-06 

.505 

.9922 

2.0320 

1.046  -.0520 

1 .0286E-05 

5.0165F-06 

5.2696E-06 

.488 

.9927 

2. 0410 

1.037  -.0600 

1 .0089E-05 

4. 81595-06 

5.2728E-06 

.477 

.9926 

2.0490 

1 .027  -.0680 

9.9662E-06 

4.7061F-06 

5.2600E-06 

.472 

.9922 

2.0570 

1.016  -.0790 

9.9622F-06 

4.7198E-06 

5.2425E-06 

.474 

.9912 

2.0740 

.995  -.1100 

1.0  1 91F-05 

4.$9l2E-06 

5. 1997E-06 

.490 

.9869 

2. 0900 

.979  -.1460 

1 . 0406E  *05 

S.2509E-06 

5.1551E-06 

.505 

.9811 

2. 0070 

.960  -.1930 

1.0581E-05 

5.4936E-06 

5.0877E-06 

.519 

.9726 

2.0240 

.967  -.2570 

1 .0747E-05 

5.6790F-06 

5.0678E-06 

.528 

.9632 

2.0330 

.978  -.2760 

1.0791E-05 

5.7143F-06 

5.0762E-06 

.530 

.9610 

2.0410 

.987  -.2950 

1 .0829E-05 

5. 752®F-06 

5.0766F-06 

.531 

.9586 

2.0590 

1.004  - - 3310 

1.0B94F-05 

5.8239E-06 

5.0702E-06 

.535 

,9536 

2.0760 

1.021  -.3670 

1 .0949E-05 

5.S926E-06 

5.0568E-06 

.538 

.9483 

2.0040 

1.035  -.4110 

1.1003E-05 

S.9777F-06 

5.0250E-06 

.543 

.9413 

3.0120 

1.064  -.4660 

1 • 1 065E-05 

6. 071 6E-06 

4.9936E-06 

.549 

.9328 

3.0300 

1.104  -.5220 

1 .1  124E-05 

6.1S99E-06 

4.9648E-06 

.554 

.9242 

3.0490 

1.170  -.5840 

1.1  187E-05 

6.247>E-06 

4.9399E-06 

.558 

.9148 

3.0580 

1 .217  -.6120 

1.1 21 8E-05 

6.2845F-06 

4.9337E-06 

.560 

.9104 

3.0670 

1.273  -.6270 

1.1245F-05 

6.30375-06 

4.9409E-06 

.561 

.9077 

3. 0770 

1.330  -.6250 

1.1264E-05 

6.3035E-06 

4.9604E-06 

.560 

.9070 

3.0860 

1.384  -.6140 

1.1 279F-05 

6.2975F-06 

4.9812E-06 

.558 

.9071 

3.0960 

1.432  -.5930 

1 • 1 290E-05 

6.2841F-06 

5.0054E-06 

.557 

.9080 

3.1060 

1.469  -.5640 

1.1 296F-05 

6.26345-06 

5.0324E-06 

.555 

.9099 

3.1150 

1.499  -.5390 

1.1 302F-O5 

6.2498E-06 

5.0522E-06 

.553 

.9112 

3.1250 

1.527  -.5150 

1.1 308E-05 

6.2395E-06 

5.0687E-06 

.552 

.9121 

3.1350 

1.553  -.4910 

1.131 4F-05 

6.231 2F-06 

5.0829E-06 

.551 

.9129 

3.1450 

1.578  -.4650 

1.1320F-05 

6.2233E-06 

5.0966E-06 

.550 

.9136 

3.1550 

1.603  -.4380 

1.1 326F-05 

6.2174F-06 

5.1 084E-06 

.549 

.9141 

3. 1 65u 

1.623  -.4060 

1.1331F-05 

6,20745-06 

5. 1 236E-06 

.548 

.9152 

3.1750 

1.640  -.3720 

1 . 1 336F-05 

6.  19645*06 

5.1 392E-06 

.547 

.9164 

3.1850 

1.6S0  -.3360 

1 .1341F-05 

651831 E-06 

5. 1 575E-06 

.545 

.9181 

3.7050 

1.653  -.2670 

1 . 1349F-05 

6.1521F-06 

5. 1 972E-06 

.542 

.9222 

3.7260 

1.841  -.2100 

1.1 359F-05 

6. 12245-06 

5.2362E-06 

.539 

.9264 

3.7790 

1.593  -.1210 

1.1 382F-05 

6.07125-06 

5.3113E-06 

.533 

.9343 

3.7330 

1*545  -.0740 

1 .1  406E-O5 

6.06245-06 

5, 3439E-06 

.531 

.9384 

7.7900 

1.511  -.0410 

1.1 4J3F-05 

6.18285-06 

5.2497E-06 

.541 

.9342 

3.4480 

1.476  -.0240 

1.1459F-05 

6.5017F-06 

4.9572E-06 

.567 

.9220 

3.5090 

1 .449  -.0150 

1. 1485E-05 

7.0193F-06 

4.4652E-06 

.61  1 

.9047 

3.5710 

1.428  -.0090 

1. 1500F-05 

7.7899E-06 

3,71 04E-06 

.677 

.8832 

3.4360 

1.410  -.0080 

1.151 OF-05 

8.00985-06 

3.5009E-06 

.696 

.8796 

3.7040 

1.398  -.0070 

1.1524E-05 

8.26765-0 

3.2566E-06 

.717 

.8742 

3.«460 

1 .376  -.0070 

1.1 602F-05 

8.39715-06 

3.2050E-06 

.724 

.8757 

4.O000 

1061  -.0100 

1.1 61 4F-05 

7.9791E-06 

3.7345E-06 

.678 

,8896 

4.1670 

1046  -.0120 

1.1 677F-05 

7.69325-06 

3.9835E-06 

.659 

.8962 

4.7550 

1 .343  -.0150 

1 . 1 700E-05 

7.3963E-06 

4.3032E-06 

.632 

,9049 

4.7480 

1037  -.0200 

1.171 7F-05 

7.03305-06 

4.6838E-06 

.600 

.9164 

4,7860 

1037  -.0230 

1.1721F-05 

6.9689F-06 

4.8524E-06 

.586 

.9215 

4,4440 

10  38  -.0260 

l • 1 729F-05 

6.7465E-06 

4.9823E-06 

.575 

.9253 

4.5050 

1040  -.0290 

1.1 739F-05 

6.6485F-06 

5.0900E-06 

.566 

.92S4 

4.5450 

1044  -.0300 

1.1 743F-05 

6.6246F-06 

5.11 83E-06 

.564 

.9285 

4.4080 

1047  -.0250 

1.1 763E-05 

6. B4B75-06 

4.9147E-06 

.582 

.9200 

4.4510 

1.347  -.0210 

1.1782E-05 

7.08395*06 

4.6927E-06 

.602 

.9120 

4.7662 

1041  -.0150 

1 , 1 826E-05 

7.64845-06 

4.1772E-06 

.647 

.8966 

5.O00O 

1.327  -.0720 

1 . 190.3F-05 

9. 1623E-06 

3.7407E-06 

.686 

.8866 

5.7630 

1012  -.01<*0 

1.1 997E-05 

8.0675F-06 

3.9292E-06 

.672 

.8928 

5.5560 

1.296  -.0240 

1 .2070E-05 

7.3673E-06 

4.7O25E-06 

.610 

.9154 

5.7140 

1.287  -.0350 

1.2081F-05 

6.92945-06 

5.  1527E-06 

.573 

.9306 

5.7470 

1.285  -.0380 

1 .2079F-05 

6.94305-06 

5.2360E-06 

.567 

.9337 

5.7800 

1.284  -.0420 

1 .2071F-05 

6.73985-06 

5. 3308E-06 

.558 

.9371 

5.«140 

1.283  -.0460 

1 .2062F-05 

6.65295-06 

5.4089E-06 

.552 

.9400 

5.B480 

1.283  -.0520 

1.2043E-05 

6.53895-06 

5.5038E-06 

.543 

.9434 

5.0820 

1.285  -.0570 

1.2027E-05 

6. 46075-06 

5.5664E-06 

.537 

.9454 

6. 0610 

1001  -.0690 

1.201 1F-05 

6.33405-06 

5.677  IE-06 

.527 

.9464 

-1.0 
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Table  A-9 


ICE  16. 


HAyLOH 

REFft- 

INDICES 

EXTNC-COEF 

SCATt-COEF 

ABS9P-C0EF 

AL0DO 

0 

2.8000 

1.225 

-.0040 

4 . 3B34E-0S 

2. 

8673E-05 

1.5161E-05 

.654 

.9464 

2.5640 

1.212 

-.0040 

4.3913E-05 

2. 

888?E-05 

1.5031E-05 

,658 

.9494 

2.6320 

1.192 

-.0050 

4.3942E-05 

2. 

7792E-05 

1 .6150E-05 

.632 

.9571 

2.7030 

1.162 

-.0080 

4.4033E-05 

2. 

5656E-05 

1.8376E-05 

.583 

.9698 

2.7780 

1.108 

-.0200 

4.4153E-05 

2. 

31 5?E-05 

2.1 001E-05 

.524 

.9876 

2.8170 

1.065 

-.0380 

4. 3533E-05 

2. 

2093E-05 

2.1440E-05 

.507 

.9940 

2.8250 

1 .056 

-.0440 

4.3063E-05 

2. 

1650E-OS 

2.1414E-05 

.503 

.9946 

2.8330 

1.046 

-.0520 

4.2470E-05 

2. 

L126E-05 

2.1 344E-05 

.497 

.9948 

2.8410 

1.037 

-.0600 

4.2028E-05 

2. 

0775E-OS 

2.1 253E-05 

.494 

.9947 

2.8490 

1.027 

-.0680 

4.1716E-05 

2. 

0!>3?~-05 

2.1 134E-05 

.493 

.9943 

2.RS70 

1.016 

-.0790 

4.  1625E-05 

2. 

064)t-&J 

2.0984E-05 

.496 

.9933 

2.8740 

.995 

-.1100 

4.1941E-05 

2* 

1292E-05 

2.0649E-05 

.508 

.969? 

2.8900 

.979 

-.1460 

4 , 2292E-05 

2. 

1941 E-05 

2.0351C-05 

.519 

.9836 

2.9070 

.960 

-.1930 

4.2593F-05 

2. 

2606E-05 

1.9987E-05 

-.531 

.9753 

2.9240 

.967 

-.2570 

4.2900E-05 

2. 

3124E-05 

1.9775E-0S 

.539 

.9661 

2.9330 

.978 

-.2760 

4.2984E-05 

2. 

3218E-05 

1.9766E-05 

.540 

.9639 

2.9410 

.987 

-.2950 

4 , 3058E-05 

2. 

3325E-05 

1 .9732E-0S 

.542 

.9614 

2.9590 

1.004 

-.3310 

4 . 3182E-05 

2. 

3532E-05 

1 .9650E-05 

.545 

.9565 

2,9760 

1.021 

-.3670 

4 • 3288E-05 

2. 

3739E-05 

1.9549E-05 

.548 

.9512 

2. 99''0 

1.035 

-.4110 

4.3392E-05 

2. 

40 1 0F-05 

1.9382E-05 

.553 

.9442 

3.nl  -‘O 

1.064 

-.4660 

4.3514E-05 

2. 

43095*05 

1.9205E-05 

.559 

.9356 

3 . 0 3 0 

1.104 

-.5220 

4 . 3630E-05 

2. 

459i.F-05 

1 .9040E-05 

.564 

.9269 

3.0490 

1.170 

-.5840 

4 • 3757E-05 

2. 

4873E-05 

1.8B83E-05 

.568 

.9173 

3.0580 

1.217 

-.6120 

4.3819E-05 

2. 

4991 E-05 

1 .8829E-05 

.570 

.9129 

3.0670 

1.273 

-.6270 

4.3873E-05 

2. 

504*5-05 

1 .8829E-05 

.571 

.910? 

3.0770 

1.330 

-.6250 

4.391 1 E-05 

2. 

5029E-05 

1 .88B3E-05 

.570 

.9095 

3.0860 

1.384 

-.6140 

4 . 394 1 E-05 

2. 

49955-05 

1.B946E-05 

.569 

.9096 

3.0960 

1.43? 

-.5930 

4.3963F-05 

2. 

4936F-05 

1.9027E-05 

.567 

.9105 

3.1060 

1.469 

-.5640 

4.3975K-05 

2. 

4853E-05 

1.9122E-05 

.565 

.9125 

3.1150 

1.499 

-.5390 

4 . 3986E-05 

2. 

4795F-05 

1.91 91E-05 

.564 

.9138 

3.1250 

1.527 

-.5150 

4 • 3998F-05 

2. 

4751 E-05 

1.9247E-05 

.563 

.9147 

3.1350 

1.553 

-.4910 

4.4009E-05 

2. 

47UE-05 

1.9295E-05 

.562 

.9155 

3.1450 

1.578 

- . 4650 

4 . 4020E-05 

2. 

467<?f-05 

1 .Q341E-05 

.561 

.9163 

3.1550 

1.603 

-.4380 

4. 4031 E-05 

2. 

4651 E"05 

1.9379E-05 

.560 

.9168 

3.1650 

1.623 

•«  4060 

4.4040E-05 

2. 

4608E-05 

1 .9432E-05 

.559 

.9180 

3.1750 

1.640 

-.3720 

4.4048E-Q5 

2. 

4563E-05 

1.9485E-05 

.558 

.919? 

3.(850 

1.650 

-.3360 

4.4055E-05 

2. 

4506E-C5 

1.9549E-05 

.556 

.9209 

3.i>05O 

1.653 

-.2670 

4 • 4069E-05 

2. 

437QE-05 

1 .9690E-05 

.553 

.9251 

3.H260 

1.641 

-.2100 

4.4084E-05 

2. 

4255^-05 

1.9829E-05 

.550 

.9294 

3.5790 

1.593 

-.1210 

4 . 41 26F-05 

2. 

4015E-05 

2.011  IE-05 

.544 

.9381 

3.9130 

1.545 

-.0740 

4.4171E-05 

2. 

3851 E-05 

2.031 9E -05 

.540 

.9445 

3.9900 

1.511 

-.0410 

4 . 422 1 E-05 

2. 

3341 E-05 

2.0380E-05 

.539 

.9474 

3.4480 

1.476 

-.0240 

4.4264E-05 

2. 

4115E-05 

2.0149E-05 

,545 

.9469 

3.6090 

1.449 

-.0150 

4.431  5E-0S 

2. 

4912E-05 

1 . 9403E-05 

.562 

.9410 

3.5710 

1.428 

-.0090 

4.4357E-05 

2. 

6657E-05 

1 .7700E-05 

.601 

.927? 

3.4360 

1.410 

-.0080 

4 . 4368E-05 

2. 

7241E-05 

1.71 25E-05 

.614 

.9250 

3.7040 

1.398 

-.0070 

4 • 4464E  *05 

2. 

3045E-05 

1.6419E-05 

.631 

.9203 

3 = *460 

1.376 

-.0070 

4.4520F-05 

2. 

924*5-05 

1.6277E-05 

.634 

.9219 

4.4000 

1.361 

-.0100 

4 • 457 1 E -05 

2. 

6675E-05 

1 .7894E-05 

.599 

.9347 

4.1670 

1.348 

-.0120 

4 • 4672F-05 

2. 

6100E-05 

1.R571E-05 

.584 

.9400 

4.9550 

1.343 

-.0150 

4.4710E-05 

2. 

5361 E-05 

1 • 9349E-05 

.567 

.9464 

4.  1480 

1.337 

-.0200 

4 . 4 76  3E-05 

2. 

4585E-05 

2.0178E-05 

.549 

.9531 

4.9860 

1.337 

-.0230 

4.4762E-05 

2. 

4280E-05 

2. 0482E-05 

.5*2 

.9558 

4.4440 

1.338 

-.0260 

4 . 4775F-05 

2. 

4075E-05 

2.0700E-05 

.538 

.9576 

4.5050 

1.340 

-.0290 

4 , 490  7E-05 

2. 

3933E-05 

2.0873E-05 

.534 

.9588 

4.5450 

1.344 

-.0300 

4 . 4828E-05 

2. 

391  IE-05 

2.0915E-05 

.533 

.9586 

4,4080 

1.347 

-.0250 

4.4873E-05 

2. 

4269E-05 

2. 0603E-05 

.541 

.9546 

4.4510 

1.347 

-.0210 

4.4892E-05 

2. 

4697E-05 

2.0205E-05 

.550 

.9508 

4.7662 

1.341 

-.0150 

4.4978F-05 

2. 

5963E-0S 

1.911 4E-05 

.575 

.9414 

S.nOOO 

1.327 

-.01-80 

4.5132E-05 

2. 

703*E-05 

1 .P047E-05 

.600 

.9333 

5.5630 

1.312 

-.0140 

4.5253E-05 

2. 

6634E-05 

1 . 861 9E-05 

.589 

.9375 

5.5560 

1.296 

-.0240 

4.5364F-05 

2. 

487?E-05 

2.0492E-05 

.548 

.9528 

5.7140 

1.207 

- • 0 3S0 

4 . 541 OE -05 

2. 

405) P-05 

2, 1 359E-05 

.530 

.9607 

5.7470 

1.285 

-.0380 

4 . 541 6E-05 

2. 

3917E-05 

2 . 1 499E-05 

.527 

.962? 

5.7800 

1.284 

-.0420 

4.541 7E-05 

2. 

37  7 1 f-05 

2.1 645E-05 

.523 

.9636 

5.8140 

1.283 

-.0460 

4 • 541 6E-05 

2. 

3656E-05 

2.1759E-05 

.521 

.9648 

5.8480 

1.283 

-.0520 

4.5400F-05 

2. 

3522E-05 

2.1 878E-05 

.51  B 

.9660 

5.8820 

1.285 

-.0570 

4.5388E-05 

2. 

3445E-05 

2.1 943E-05 

.517 

.9666 

6.0610 

1.301 

-.0690 

4 • 540  SF-05 

2. 

3389E-05 

2.2025E-05 

.515 

.9657 

1 .n 
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32. 

9EFR-IWICES 

Table  A- 10 

EXT9C-C0EF  SCATT-COEF 

absrp-coef 

ALBDO 

0 

2. soon 

1.225 

.0040 

1 .7046F-04 

9.B810E-05 

7,1 653E-05 

.580 

.9684 

2.5640 

1.212 

.0040 

T.7059E-04 

9.9205E-05 

7.1 3B8E-05 

.582 

.9695 

2. *.320 

1.192 

.0050 

1.70665-04 

9.6356E-05 

7.4302E-05 

.565 

.9750 

2. ’030 

1.162 

.0080 

1.7078E-04 

9.164?F-05 

7.9134^-05 

.537 

.9832 

2.1  “30 

5*138 

.0200 

1.7073E-04 

B.7162E-05 

B.3571E-05 

.511 

.9921 

? .**1  ?'•:■ 

! -»C£>5 

.0380 

1.6979E-04 

8.557ftF-0S 

8,421 4E-05 

.504 

. 9954 

i',fl25(.' 

1.056 

.0440 

1.6916E-04 

8.5100E-05 

8.4060E-05 

.503 

.9957 

2.«j3(f 

I.C46 

.0520 

1 .602BE-04 

B.4530E-05 

8.3755E-05 

.502 

.9958 

2.04io 

i .037 

.0600 

1 .6756E-04 

8.4172E-05 

8.3385E-0S 

.502 

.9956 

2.0490 

1.027 

.0680 

1 .6699E-04 

8.40505-05 

8.29415-05 

.503 

.9952 

2.8570 

1.016 

.0790 

1 .687SE-04 

B. 44045-05 

8.2346E-05 

.506 

.9942 

2.0740 

.995 

.1100 

1.671 5E“04 

8.6219E-05 

B.0930E-05 

.516 

.9903 

2.0900 

.979 

.1460 

1.6772F-04 

8.  3093E-05 

7.9630E-05 

.525 

.9847 

2.0070 

.960 

.1930 

1 .6826E-04 

9.0173E-05 

7.8088E-05 

.536 

.9765 

2.9240 

.967 

.2570 

1 .68e4E-04 

9.1B14E-05 

7.7030E-05 

.544 

.9672 

2.0330 

.978 

.2760 

1 .6901E-04 

9.2098E-05 

7.6910E-05 

.545 

.9650 

2.0410 

.987 

.2950 

1.6915E-04 

9.24415-05 

7.671  IE-05 

.546 

.9626 

2.0590 

1.004 

.3310 

1 .6940E-04 

9.3116E-05 

7.6280E-05 

.550 

.9576 

2.0760 

1.021 

.3670 

1.6961F-04 

9. 3809E-05 

7.5798E-05 

.553 

.9^23 

2.0040 

1.035 

.4110 

1 .69815-04 

9.47455-05 

7.5068E-05 

.558 

.9452 

3.0120 

1.064 

.4660 

1 .7006E-04 

9. 57935-05 

7.4274E-05 

.563 

.9366 

3.n300 

1.104 

.5220 

1 .7029E-04 

9.67665-05 

7 . 3525E-05 

.568 

.9278 

3.0490 

1.170 

.5840 

1 .7055F-04 

9.7761E-05 

7.2786E-05 

.573 

.9181 

3.0580 

1.217 

.6120 

1 .706BE-04 

9.8171E-05 

7 . 2504E-05 

.575 

.9136 

3.0670 

1*273 

.6270 

1.7078E-04 

9.93435-05 

7.2441E-05 

.576 

.9104 

3.0770 

1.730 

.6250 

1.7087F-04 

9.82685-05 

7.2597E-05 

.575 

.9102 

3.0860 

1.384 

.6140 

1.7093F-04 

9.9125F-05 

7.2802E-05 

.574 

.9102 

3.0960 

1.432 

.5930 

1 .7097F-04 

9.78915-05 

7.3081E-05 

.573 

.9112 

3.1060 

1.469 

.5640 

1.710 0F-04 

9.75715-05 

7.3428E-05 

.571 

.9132 

3.1150 

1.499 

.5390 

1.7102F-04 

9.7345E-05 

7.36775-05 

.569 

.9145 

3.1250 

1.527 

.5150 

1 .7105E-04 

9.717,c-05 

7 « 3B76E-05 

.568 

.9155 

3.1350 

1.553 

.4910 

1.71075-04 

9.  70255-05 

7 .40455-05 

.567 

.9163 

3.1450 

1.578 

.4650 

1.710  9F-04 

9.6996r“05 

7.4207E-05 

.566 

.9171 

3.1550 

1.603 

.4390 

1.7112F-04 

9.677.35-05 

7.4342E-05 

.566 

.9176 

3.1650 

1.623 

.4060 

1.711 3F-04 

9.66025-05 

7.4531E-05 

.564 

.9188 

3.1750 

1.640 

.3720 

1.711 5F-04 

9.64255-05 

7.4725E-05 

.563 

.9200 

3.1850 

1.650 

.3360 

1.71165-04 

9.62005-05 

7.4962E-0S 

.562 

.9218 

3 • ?050 

1.653 

.?670 

1 .7119E-04 

9.57025-05 

7.54R5E-05 

,559 

.9261 

3.?260 

1.641 

.7100 

1 .7121F-04 

9.52185-05 

7.5996E-05 

.556 

.4304 

3.3790 

1.593 

.1210 

1 .7130F-04 

9.42685-05 

7.7029E-05 

.550 

.9393 

3.7330 

1.545 

.0740 

1 .71 39F-04 

9.35675-05 

7 . 7822E-05 

.546 

.9459 

3.7000 

1.511 

.0410 

1.7 1 49F-04 

9.31445-05 

7.8348E-05 

.543 

.9502 

7.4480 

1*476 

.0740 

1.71595-04 

9.29555-05 

7 . H638E-05 

.542 

.9537 

7.O090 

1.449 

.0150 

1 .7169F-04 

9.34965-05 

7.81 92E-05 

• 545 

.9545 

3.0710 

1.428 

.0090 

1.71785-04 

9.59315*05 

7.58S1E-05 

.558 

.9505 

7.4360 

1.410 

.0080 

1 .7189F-04 

9.  6925c-05 

7 . 4963E-05 

.564 

.9501 

3.’040 

1.398 

.0070 

1 .7196F-04 

9.8396^-05 

7 . 3574E-05 

.572 

.4482 

7.0460 

1.376 

.0070 

1 ,721  7E-04 

9.87385-05 

7 . 3429E-05 

.574 

,9494 

4.O000 

1.361 

.0100 

1 .7232F-04 

9.54675-05 

7. 68565-05 

.554 

. 957R 

4.1670 

1.348 

.0120 

1 .7254F-04 

9.43625-05 

7.8178E-05 

.547 

.4613 

4.2550 

1.343 

.0150 

1.7263F-04 

9.31765-05 

7.9453E-05 

.540 

.9642 

4.7480 

1.337 

.0200 

1 .7271F-04 

9.21155-05 

8.0598E-05 

.533 

.9670 

4.7860 

1037 

.0230 

1.72755-04 

9.17915-05 

8 , 0954E-05 

.531 

.9678 

4 .4440 

1.338 

.0260 

1 .7280E-O4 

9.16075-05 

8,11  93E-05 

.530 

.9681 

4. 0050 

1.340 

.0290 

1.7286E-04 

9.15035-05 

8. 1 354E-05 

.529 

.968? 

4.0450 

1.344 

.0300 

1 .7290F-04 

9.15365-05 

8.1 362E-05 

.529 

.9678 

4.4080 

1.347 

.0250 

1 .7299E-04 

9.19155-05 

8, 1 075E-05 

.531 

.9667 

4.4510 

1.347 

.1)210 

1 .7305E-04 

9. 23S(.c-05 

8.0671E-05 

.534 

.9656 

4. 366? 

1.361 

.0150 

1 .7323E-04 

9.39195-05 

7.931 3E-05 

.542 

.9627 

5. 0000 

1 . 327 

.01-20 

l .7352F-04 

9. 58475*05 

7 . 7670E-05 

.552 

.9599 

5.2630 

1.312 

.0140 

1 .7382E-04 

9.4917E-05 

7.P910E-05 

.546 

.9630 

5.0560 

1.296 

.0240 

1 .7409E-04 

9.20786-05 

8.201  IE-05 

.529 

.9700 

5.7140 

1.287 

.0350 

1.7416F-04 

9.1059E-05 

8.3101E-05 

.523 

.9728 

5.7470 

1.285 

.0380 

1 .741 7E-04 

9.09H  5-05 

8. 32625-05 

.522 

.9733 

5.7800 

1.284 

.0420 

1 .7418E-04 

9. 07745-05 

8. 34015-05 

.521 

.9736 

5.»140 

1.283 

.0460 

1 .7418F-04 

9.06785-05 

8.35045-05 

.521 

.9738 

5.0480 

1.283  -.0520 

1 .7416F-04 

9.05805-05 

8.35785-05 

.520 

.9739 

5.0820 

1.285 

.0570 

1.7-155-04 

9.05656-05 

8.35855-05 

.520 

.9737 

6.0610 

.0 

1 . 30 1 -.0690 

1 .7426E-04 

9.07965-05 

8.34595-05 

.52! 

.9719 

145 


APPENDIX  B 

RADIOSONDE  DATA/ LOWTRAN  INPUTS 

This  appendix  gives  meteorological  profiles  used  in  computing  LOWTRAN 
transmittances  and  radiances  for  the  model  validations.  Tables  B—  1 through  B— 3 
correspond  to  AFGL  missions  that  provided  cloud  radiance  data  (see  Section  3). 
Tables  B— 4 and  Br-5  correspond  to  missions  that  observed  specular  sunglints  from 
lakes  and  ponds  (see  Section  4).  The  remaining  tables  correspond  to  missions  that 
observed  industrial  sites  (see  Section  6), 

Tables  B—l  through  B— 5 are  listings  of  the  input  cards  to  the  modified 
LOWTRAN  4 code  (the  listings  include  additional  cards  that  identify  columns  and 
fields,  as  well  as  a header  card).  The  pressure,  temperature  and  dew  point 
profiles  are  based  on  radiosonde  data  provided  by  AFGL. 

Tables  Br-6  through  B— 9 give  the  raw  radiosonde  data  provided  by  AFGL 
for  mission  31.r'  and  816. 


Table  B-l  INPUTS  TO  MODIFIED  LOWTRAN  4 FOR  MISSION  726 


MISSION  726 

12345678901234567890123456789012345678901234567890123456789 


2 21 


0.  12.6  8.88 

ALT  PRES  TEMP 

(KM)  (MB) <C> 


0. 13  1000.  17.2 

.40  972.  14.2 

.76  929.  12.4 

<9  914.  12.6 

1.3  873.  15.6 

1.5  850.  15. 

2.52  753.  10.2 

3*125  700.  9.6 

4.7  576.  “0*3 

5*83  500.  -9.5 

6*46  460*  “14.1 

6.58  453.  -13.7 

7.51  400.  -21. 

9.57  300.  -38. 

10.8  250.  -48. 

12.25  200.  -55. 

14.07  150.  -61. 

16.5  100.  -64. 

20.  54.  -64. 

30. U. -64. 

1234567890 1^34567890 123451 

1. R?Q  0.9144  95. 

1800.  ' 3600.  5. 


DEM  PT 
(C) 

i 3” 

12.5 

12.3 

11. 

6.9 

-14.4 

-15. 

-19.8 

-20.4 

-30.3 

-39.5 

-34.1 

-43.7 

-51. 

-68. 

-68. 

-74. 

-77. 

-77. 

-76. 

-72. 

5oT?34B6789o  123456789(^1 23456789 

20. 


Table  B-2  INPUTS  TO  MODIFIED  LOWTRAN  4 FOR  MISSION  727 


MISSION  727 

1234567890 1234567890 12 34567890 1234567890 1234567890 123456789 


7 2 

2 

1 2 

17  1 

0. 

-32. 

14.25 

ALT 

pR'cS 

TEMP 

DEW  PT 

< KM ) 

(MR) 

<C> 

(C) 

0. 

1010. 

23.2 

0.2 

0.6 

940. 

23.2 

0.2 

1.34 

864. 

23.2 

0.2 

1 .45 

850. 

22. 

1.0 

3.09 

700. 

6.0 

-6.0 

5.75 

500. 

-11.5 

-22.5 

6.95 

430. 

-19.3 

-49.3 

7.4? 

400. 

-24,0 

-40.0 

9.01 

319. 

-37.5 

-47.5 

9.45 

30  or 

-41.0 

-55. 

10.66 

250. 

-51. 

-61. 

10.95 

?3°.  • 

-54.1 

-66. 

12.1 

200. 

-54.3 

-70. 

13.93 

150. 

-58.5 

-74. 

16.47 

100. 

-61.5 

-77. 

20. 

54. 

-62. 

-77. 

30. 

11. 

-62. 

-72. 

1234567890 1234567890 123456^890 1234567890 1234567890’ 23456789 
6.000  8.40  75.75  20. 

1800.  3600.  5. 
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Table  B-3  INPUTS  TO  MODIFIED  LOWTRAN  4 FOR  MISSION  817 


MISSION  817 

l 234567890 1234567890 1234567890 1234567890 1234567890 123456789 


7 2 

2 

1 

2 26  1 

0. 

-49.3 

2.6 

ALT 

PRES 

TEMP 

DEW  PT 

(KM) 

<M9) 

<C> 

(C) 

0. 

1020. 

14.4 

10.7 

.104 

1000. 

14.4 

10.7 

.70 

931. 

14.4 

10.7 

1.05 

895. 

10.4 

8.7 

1.477 

850. 

7.8 

6.6 

? .6? 

737. 

2.6 

-0.8 

3.059 

700. 

-0.5 

-1.8 

3.55 

656. 

-4.3 

-5.6 

3.90 

622. 

-5.7 

-17.7 

4.48 

584. 

-8.5 

-17.5 

4.91 

552. 

-12.1 

-14.7 

5.3 

526. 

-13.9 

-18.9 

S.67 

500. 

-16.7 

-20. B 

6.13 

469. 

-20.5 

-26.5 

7.06 

430. 

-25.5 

-32.5 

7.31 

400. 

-29.5 

-39.5 

7.5? 

388. 

-31.1 

-40  . 

8.54 

334. 

-39.7 

-46. 

9.29 

300. 

-45.5 

-56. 

\0. ?S 

259. 

-53.3 

-62. 

10.48 

250. 

-50.9 

-63. 

10.71 

241. 

1-48.9 

-64.. 

11. 95 

200. 

-47.3 

-69. 

13.18 

167. 

-46.3 

-72. 

13.85 

150. 

-48.1 

-73. 

16.5 

100. 

-52.5 

-77. 

1234S67890 1234567890 1234567890 1234567890 1234567890 123456789 
10 • OSS  9.754  92.5  20. 

1800.  3600.  5. 
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Table  B-4 


INPUTS  TO  MODIFIED  LOWTRAN  4 FOR  MISSION  729 


MISSION  729 

1 2345678901234567890123456789012345678901234567890123436789 
7 2?  1 2 19 


0. 

3.5 

7.34 

ALT 

PRES 

TEMP 

DEW  PT 

( KM) 

(Mr) 

(C) 

<C> 

0. 

ion . 

23.3 

6.4 

0.09c? 

1000  . 

21.8 

1 1 .8 

O.S 

949. 

20.2 

1.2 

1 .4 

857. 

13.6 

-5.4 

1 .54 

842. 

15.0 

-15.0 

2.2 

780. 

11.2 

-18.8 

2.47 

755. 

1 1.8 

-18.2 

3.11 

700. 

7.8 

-22.2 

5.78 

500. 

-12. 

-52.2 

6.28 

470. 

-15.7 

-45 . 7 

6.9 

421  . 

-22,1 

-50. 

7.45 

400. 

-25.2 

-55. 

9.47 

300. 

-42,7 

-65. 

10.68 

250. 

-48,7 

-68. 

12.18 

200. 

-52.7 

-72. 

1 3.9f 

150. 

-58.9 

-75. 

16,47 

100. 

-63.7 

-77. 

20.00 

54. 

-64. 

-76. 

30. 00 

1 1 . 

-64. 

-72. 

1234567890 1204567890 1234567890 1234 5678*0 1234567890 123456789 
4.57?  0.  97.34  37. 
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Table  H-b 


INPUTS  TO  MO  Dll'  1 KH  i.OWTRAN  4 i'Ol,  *'s ' ON  HO' 


MISSION  803 

1 2340678901234567890123456789012345678901234567890123456789 
7 2?  1 3 18 


n . 

3.5 

5.85 

ALT 

PROS 

TEMP 

DEW  PT 

(KM) 

(MR) 

(C) 

<C) 

1 .2S 

876* 

18.3 

“4.7 

1 .500 

850. 

17.5 

-5.5 

?.0?4 

8C0. 

15.5 

-14.5 

1.137 

700. 

6,4 

-10.8 

3.50? 

662. 

3.0 

-13.4 

A .603 

56  3. 

•6.5 

-16*6 

A.  007 

561  . 

~8  • 6 

-15.4 

0.4? 

5?5. 

-7.6 

-19.3 

0.8  0 

500. 

-9.2 

-23.2 

6.44 

460. 

-12.5 

-42.5 

7.48 

400  . 

-22.7 

-52.7 

o.5? 

300. 

■*4 0 • 4 

-65.0 

]0.7'. 

250. 

-49.9 

-68. 

12,? 

200. 

-50.2 

-72. 

14.04 

1 50  • 

-58.2 

-75. 

16.0? 

100. 

-71.6 

-77. 

20. 

56. 

-61,2 

-7 e . 

30  . 

12. 

-55. 

-72. 

1 234067890 12 3456 78 90 123456 78 90 12 3456 7 89012 34567890 12 3456 789 

0.79 

1 .95 

113,63 

130  . 

1 800, 

5800. 

20. 
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Table 

B-6  RADIOSONDE  DATA  APPLICABLE 
MISSION  816,  RUNS  1-4 

TO 

MIDWAY 

23  1846Z  JUNE 

78  (Used  for  GARY 

data) 

Height* 

Pressure 

Temperature 

Dew  Point 

H(km) 

P(mb) 

T(C) 

(C) 

0. 173 

1000 





(surface) 

998 

17.8 

6.8 

988 

17.8 

4.8 

980 

21.0 

1.0 

968 

21.0 

2.0 

875 

13.0 

-1.0 

1.558 

350 

11.4 

5.4 

o!3 

8.4 

5.3 

769 

7.2 

-0.8 

743 

0.6 

-5.5 

3.156 

700 

2.0 

-1.8 

Mean  height  above  sea-level.  Heights  given  at  TTAA 
(mandatory)  levels.  Heights  at  TTBB  (significant) 
levels  may  be  obtained  by  interpolating  H vs.  In  P. 


Table  B-7  RADIOSONDE  DATA  APPLICABLE  TO 
MISSION  816,  RUNS  6-10 

FLINT  23  1849Z  JUNE  78  (Used  for  WYANDOTTE  data) 


H (km) 

V (mb) 

T (C) 

Dew  Point 

(C) 

0.178 

1000 

— 

— 

(Surface)* 

993 

10.6 

8.3 

980 

12.2 

5.2 

922 

12.2 

2.2 

1.536 

850 

9.0 

2.0 

797 

4 . 6 

-2.4 

781 

4.0 

-12.0 

3 . 1 1 9 

700 

2.2 

-27.8 

680 

2.2 

-27.8 

542 

-9.7 

-39.7 

? 

5.760 

500 

-14.9 

-44.9 

? 

451 

-21.7 

-38.7 

420 

-24.3 

-54.3 

7 

7.410 

400 

-27.1 

-57.1 

7 

9.41 

300 

-42.1 

— 

10.63 

250 

-52.5 

— 

12.04 

200 

-57.5 

— 

13.88 

150 

-55.3 

— 

16.45 

100 

-57.3 

— 

it 

Surface  data: 

DETROIT 

23 

1408Z  T = 69F, 

Dew  Point 

= 44F 

DETROIT 

23 

1611Z  T = 72F, 

Dew  Point 

= 5 OF 

Table  B-8  RADIOSONDE  DATA  APPLICABLE  TO  MISSION  815/2 


ASHLAND  22  1453Z  JUNE  78 


H (km) 

P (mb) 

T (C) 

Dew  Point  (C) 

0.160 

1000 

_ _ 

(Surface) 

990 

17.2 

13.9 

969 

16.0 

10.0 

950 

17.2 

8.2 

905 

15.8 

6.8 

1.544 

850 

12.4 

11.4 

792 

8.6 

7.5 

775 

8.6 

-0.4 

723 

3.0 

-1.2 

713 

2.0 

-8.0 

710 

4.4 

-25.6 

3.146 

700 

5.4 

-24.6 

5.81 

500 

-13.1 

-33.1  ? 

7.47 

400 

-25.3 

-55.3  ? 

9.48 

300 

-43.3 

— 

10.69 

250 

-49.9 

— 

12.13 

200 

-52.3 

— 

13.98 

150 

-55.1 

— 

16.55 

100 

-59.3 

— 

154 


it  ’4li.r_,yiA  ».l 


•TTT 


Table  B-9 

RADIOSONDE  DATA 

APPLICABLE 

TO  MISSION  815/4 

BIRMINGHAM 

22  1454Z  JUNE  78 

(Used  for 

FAIRFIELD  data) 

H (kin) 

P (mb) 

T (C) 

Dew  Point  (C) 

0.144 

1000 

— 

/ 

(Surface) * 

995 

22.0 

20.5 

984 

22.2 

19.8 

973 

23.2 

19.3 

962 

24.0 

19.3 

947 

22.2 

17.5 

915 

20.2 

16.8 

882 

18.8 

11.8 

877 

18.6 

11.8 

866 

17.8 

11.8 

1.552 

850 

17.0 

9.0 

832 

15.8 

6.8 

826 

15.8 

6.8 

812 

14.0 

9.4 

783 

13.0 

3.0 

763 

11.2 

5.2 

734 

9.4 

2.4 

728 

9.0 

0.0 

3.182 

700 

6.8 

-0.2 

Near  surface  use: 

HUNTINGTON  13Z 


T = 66F,  Dew  Point  = 55F 


APPENDIX  C 

DEVELOPMENT  OF  SUNGLITTER  MODEL 


A successful  model  for  predicting  the  IR  spectral  radiance  of  sunlight 
reflected  by  a rough  ocean  surface  or  other  body  of  water  has  been  developed  using 
the  Cox-Munk  (Reference  Cl)  geometric  optics  approximation.  Although  the  model 
has  been  verified  by  comparison  with  ocean  sunglitter  data  obtained  by  the  AFGL 
Background  Measurement  Program,  one  may  also  apply  the  results  to  computations 
involving  the  specular  reflectance  of  the  sun  from  smooth  or  convex  surfaces. 
This  appendix  outlines  a step-by-step  computation  procedure  that  serves  to  define 
the  model  as  well  as  describe  the  details  of  its  implementation. 

Reflection  Geometry 


A cartesian  coordinate  system  is  used  to  define  the  reflection  geometry 
with  the  u and  v axes  defining  the  horizontal  plane  wf  the  undisturbed  ocean 
surface  (Figure  C— 1).  The  Z axis  points  vertically  upward  and  the  u axis  is 
directed  upwind.  The  locations  of  the  sun  and  the  observer  are  given  by  the  solar 
zenith  angle  <J*  , the  solar  azimuthal  angle  V with  respect  to  the  upwind  direction 
and  the  corresponding  zenith  angle  $ and  azimuthal  angle  a of  the  observer.  The 
orientation  of  the  reflecting  facet  of  the  surface  is  given  by  the  tilt  angle  p and 
the  local  angle  of  incidence  u)  at  the  fa^et.  The  tilt  angle  points  in  the 
direction  of  the  local  surface  normal  and  is  related  to  the  scattering  geometry 
angles  (<J>,  V;  a, £5)  by  the  expression  (Reference  C2). 

cos  p = [cos  $ + cos  B]/  2 cos  co  (Cl) 

Likewise,  the  local  angle  of  incidence  can  be  expressed  in  terms  of  the  scattering 
geometry. 

cos  ( 2 0))  = cos  $ cos  0 - sin  $ sin  <p  cos  (v-  a)  (C2) 
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^OBSERVER 


/^V 


c Zz> 


REFLECTING  FACET 
WITH  SLOPES  Z , Z 

U 'f 


Figure  C-l  SEA  SURFACE  REFLECTION  COORDINATE  SYSTEM.  The  observer  is  located 
at  azimuthal  angle  a and  zenith  angle  3.  The  sun  orientation  is 
given  by  azimuthal  angle  V and  zenith  angle  <f>. 


\ "PWg JUJ ; ; ?ip ygj* 


"^FT’w’Tf'  .'TO'TryiT 


It  will  also  be  necessary  to  relate  the  wave  slope  components  for  specular 
reflection  to  the  scattering  geometry.  These  results  will  be  useful  when  deter- 
mining the  probability  of  occurrence  of  slope  components  Z^  and  Z^. 

Zu  “ 3Z/3u  = (sin  0 cos  a -sin$  cos  \>)(cos4>  + cosg)  * (C3) 

Zv°3Z/3v=(sin$sina-sin4isinV)(cos4>+co8g)  * (C4) 

Basic  Equations  for  Reflected  Components 


The  geometric  optics  expression  for  the  spectral  reflected  sea  radiance 
is  given  by: 

H 

N (X)  - P (w,A)  (X)  T (X)  S (U))  (a,B;y,<t>)  P (Z  , Z , W)  (C5) 

a a <1  U V 


The  parameters  in  this  expression  are: 


N - spectral  sea  radiance  for  polarization  a and  wavelength  X 

8 

a - polarization  state  of  the  reflected  radiation  which  is  defined  to  be 
either  vertical  (a  = v,  parallel  to  plane  of  incidence)  or  horizontal  (a 
= h,  perpendicular  to  plane  of  incidence) 

Hq  - spectral  solar  irradiance  outside  the  earth's  atmosphere 

T - atmospheric  transmission  along  with  path  defined  by  a line  extending 

from  the  sun  to  the  surface  and  then  to  the  receiver 

S - a slope  shadowing  factor  which  prevents  N from  becoming  infinite  at  the 

horizon  due  to  the  sec  g term  in  the  geometric  factor  g (S  = 0 when  sec  g 
= 0)  [Reference]c3  . 

g - a geometric  scattering  factor  which  is  defined  later 


P - probability  of  the  occurrence  of  a particular  wave  slope  having  com- 
ponents Z^  and  for  a wind  velocity  W 

- reflectance  of  the  surface  for  incident  polarization  state  a and  local 
angle  of  incidence  to. 

If  one  makes  the  assumption  that  the  receiver  response  is  insensitive  to 
the  polarization  of  the  source,  then  the  total  spectral  reflected  solar  radiation 
at  the  receiver  is  just  the  sum  of  the  two  possible  polarization  states  (v  = 
vertical  and  h = horizontal). 


N (X)  = N (X)  + N.  (X) 
sun  v h 


(C6) 


A complete  definition  of  Nfl(X)  requires  further  expressions  which 
completely  define  the  parameters  which  appear  in  the  Equation  (C5). 

A geometric  scattering  factor  is  used  to  account  for  the  angular 
distribution  of  the  scattered  sunlight  (Reference  C2). 


g(d, 3; v, = sec  0 


l + cos  0 cos  <{>  - sin  0 sxn  $ cos  (V  - a) 


(C7 ) 


-4 


(cos  0 + cos  ) 


Saunders  (Reference  C3)  derives  a simplified  expression  for  the  surface 
shadowing  factor: 


S = 2 1 + erf  (V)  + ( 


v/iT)  1 erf  (-V2)J 


-1 


(C8) 


where  V = o ^ tan  <}> , O i*»  the  rms  surface  slope,,  and  erf  (v)  denotes  the  error 
function. 


The  entire  expression  for  N (X)  requires  a knowledge  of  the  probability 

3 

of  occurrence  of  wave  slopes  with  the  correct  orientations.  Cox  and  Munk 
(Reference  Cl)  have  measured  this  function  for  surface  wind  velocities  up  to 
14  m/s  (at  12  m above  sea  level). 
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H 4 


niu.^.ll  . 


nmsxp--,r:-r-  r 


P (z  ,z  ,w) 

u v 


-1 


exp 


^-l/2(£2  + 02)j  l - c(w)j 


(C9) 


The  rms  values  of  surface  slopes  in  the  u and  the  v directions  are  given  by  0^ 
and  0v;  5 and  n are  normalized  slopes  that  describe  the  deviation  of  the  slopes  of 
the  waves  from  a Gaussian  distribution. 


c(w)  " i c2i  (S2  " + l c03  <n3  - 3n) 

■ h c40  (S4  - 6c2  + 3)  - £ c22  (i-2  - l)(n2  - 1)  (CIO) 


It  is  important  to  note  that  Equation  C9  is  only  valid  for  values  of  £ 
and  n that  are  less  than  2.5  in  magnitude  and  for  <J>  < 55°  because  of  experimental 
limitations  on  the  surface  measurements.  Furthermore  Cox  and  Munk  measured  values 
of  rms  slope  components  for  a clean  surface  of 

a 2 = 0.003  + 0.00192  W 

v 

a 2 - 0.000  + 0.00316  W (Cl 1 ) 

u 

a2  = 0.003  + C. 00512  w 


along  with  coefficients  determining  the  skewness  and  peakedness  of  the  surface 
slope  distribution. 

C21  = 0.01  - 0.0086  W 
C03  = 0.04  - 0.033  W 
C40  = 0.40 
C22  = 0.12 
C04  =0.23 

The  most  imoortant  parameter  for  this  study  is  the  sea  surface  spectral  reflection 

coefficient  p for  polarization  a.  Theoretical  expressions  for  p for  an 
a « 

absorbing  medium  such  as  sea  water  are  given  by  Friedman  (Reference  C4)  along  with 
a tabulation  of  the  index  of  refraction  n and  the  extinction  coefficient  k of 
water  for  wavelengths  from  1.5  to  15y. 


mLfSrMr 
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(C13) 


b + c cos2  to  -2  [d(n2-k2)  + 2 ekn]  cos  to 
b + c cos2to  +2  [d(n2-k21  + 2 ekn]  cos  to 


2 

P,  b + cos  to  -2d  cos  to 

h = 2 

b -t  cos  to  + 2d  cos  to 


where 


a = 
b2  - 


c = 
d2  - 


e = 


2 2 2 

n + k - sin  to 

2 A ..2  2 
a + 4k  n 

, 2 ,2.2  ..22 

(n  -k)  +4kn 


1/2  (b  + a) 
1/2  (b  - a) 


(C14) 


(C15) 


The  spectral  reflectances  depend  upon  measurements  of  n and  k.  A recent 
tabulation  of  n (X)  and  k (X)  from  0.200  to  200  pm  has  been  performed  by  Hale  and 
Querry  (Reference  C5)  for  pure  water.  Research  by  Hobson  and  Williams 
(Reference  C6)  indicates  that  one  can  also  use  the  pure  water  data  for  ocean  water 
in  the  spectral  range  2-8  pm.  One  should  also  note  that  this  tabulation  is  only 
for  "clean"  water  that  has  no  foreign  surface  films  on  it.  The  presence  of  surface 
contaminants  will  greatly  alter  the  reflectances. 


Several  interesting  features  of  the  k (X)  and  n (X)  data  are  plotted  in 
Figure  C— 2 . The  influence  of  the  two  major  absorption  bands  at  2.94  pm  (3400  cm 
and  6.09  pm  (1640  cm  ) are  clearly  evident. 
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APPENDIX  D 

AFGL  INDUSTRIAL  SITE  SPECTRA 

This  appendix  contains  plots  of  27  AFGL  spectra  provided  to  SAI  by 

**1 

Honeywell.  The  spectra  are  plotted  over  three  regions,  from  1900-3060  cm  , 
3060-4220  cm  \ and  4220-5380  cm  *.  Also  an  atmospheric  correction  has  been  made 
to  the  spectra  over  the  1900-3060  cm  * region,  and  plots  of  these  adjusted  values 
are  included  with  each  data  set.  Thus  each  data  set  consists  of  4 plates.  The 
missing  data  in  the  adjusted  spectra  between  2290  and  2380  cm  results  from  the 
fact  that  LOWTRAN  shows  zero  transmission  in  this  region.  The  ranges  used  to  make 
the  atmospheric  corrections  are  shown  in  the  preface  of  the  appendix; 
meteorological  data  were  taken  from  radiosondes  given  in  Appendix  B. 


I 
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Table  D-l  IDENTIFICATION  OF  SPECTRA  FOR  AFGL  MISSIONS  816  and  815 


M1S8I0N/RUN 


DERIVED 
CANOE  (ft) 


Gary , IN 
2)  Jun.  78 


Wysndott.,  HI 
23  Juna  78 


Ashland,  KY 
22  June  78 


Fairfield,  At 
22  June  78 


Ashland,  KY 
22  June  78 


816/1 

11,500 

1344 

12.6 

Bright  orange  flame 
or  hot  ore  near  blaat 
furnace 

166-169 

816/1 

11,500 

1345 

00.0 

Background-narrow  bridge 
and  river 

170-173 

“*  ■ 

816/2 

11,500 

1348 

45. 1 

Hot.  a lag 

174-177 

816/2 

11,500 

1355 

17.3 

Water  background 

178-181 

816/3 

12,700 

1401 

55.5 

Hot  coke  oven 

182-185 

816/3 

15,300 

1405 

49.6 

Background'park.'.ng  lot 

186-189 

816/3 

15,300 

1406 

33.8 

Stack  flare  (gas  burnoff) 

190-193 

816/4 

14,300 

1415 

07.9 

Stecks  of  Lerge  cepecity 
rolling  .til  1 

194-197 

816/4 

14,300 

1415 

42.2 

Background 

198-201 

816/6 

15,477 

1502 

25.0 

Blast  furnace 

202-20* 

ft  x 6/6 

8,182 

1503 

28,3 

Background 

206-209 

8!6/7 

9,155 

1507 

42,7 

Coke  oven 

210-213 

816,  7 

8,913 

1510 

34.5 

Background-dark  barren 
area 

214-217 

816/8 

11,318 

1527 

34.3 

Stacks  with  railroad 
cars  behind 

218-221 

816/8 

.0,411 

1524 

14.  f 

Background 

222-225 

816/8 

10,421 

1528 

05.4 

Railroad  cars  with 
hot  metal 

226-229 

816/ 10 

10,424 

1540 

00.0 

BASF  st.sck,  nearly 
invisible  flame 

230-233 

816/10 

7,972 

1543 

14.8 

Water  background 

234-237 

815/2 

10,725 

1133 

26.5 

Flame  from  stack 

238-241 

815/2 

12,406 

1136 

31.6 

Background 

232-2*5 

815/2 

9,550 

1143 

13.2 

Refinery  complex 

246-249 

815/4 

9,608 

1257 

26,  2 

Tip  of  stack  (no  visible 
t 1 ame ) 

250-253 

815/4 

'.0,513 

1258 

51.6 

Background 

254-257 

815/4 

10,056 

1303 

28.9 

Grey  building  (hot?) 

256-261 

815/4 

8,510 

1307 

05.9 

Orange  hot  metal  (or 

262-2.65 

1136  11.6 

1 143  69,2 


flame?)  (moving?)  (in 
railroad  car?) 


Background  (barge  on  river! 
Refinery  complex 


266-26'* 

270-27.) 


♦From  wide-fiold  interferometer  (104)|  all  other  spectra  are  from  narrow-field  interferometer  (103). 
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essAifihi na  Vw.-fr  ij. as u.>^U'ds*.«(rim,.).* 'iokAieJ^SMa.  .isn  i! 


GARY,  INDIANA  (816/1)  - 13:44:12.6  - HOT  ORE 


TARGET  RADIANCE 


GARY,  INDIANA  (816/1)  - 13:44:12.6  - HOT  ORE 


TARGET  RADIANCE 


GARY,  INDIANA  (816/1)  - 13:44:12.6  - HOT  ORE 


t I 

K 

i * 
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TARGET  RADIANCE 


EflGH 


GARY,  INDIANA  (316/1)  - 13:45:0.0  - BACKGROUND 


BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/1)  - 13:45:0.0  - BACKGROUND 


1 71 


BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/1)  - 13:45:0.0  - BACKGROUND 


BACKGROUND  RADIANCE 


(816/1) 

GARY. INDIANA  - 13:45:0.0  - BACKGROUND  WITH  ATMOSPHERIC  TRANSMITTANCE 

f CORRECTION) 


1900,00  2045.00  2190.00  2335.00  2480.00  2625.00  2770.00  2915.00 

UflVENUMBER  (Cn-1) 


***  +M* 

(i-uo*ys*suo/n)33Nyiayy 


TARGET  RADIANCE 


BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/2)  - 13:55:17.3  - BACKGROUND 


BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/2)  - 13:55:17.3  - BACKGROUND 
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BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/3)  - 14:01:15.3  - HOT  COKE  OVEN 


GARY,  INDIANA  (816/3)  - 14:01:55.5  - HOT  COKE  OVEN 


TARGET  RADIANCE 


(816/3) 


GARY,  INDIANA  (816/3)  - 14:05:49.6  - BACKGROUND 


1AA 


BACKGROUND  RADIANCE 


GRAY,  INDIANfi  - 14:05:49.6  - BACKGROUND  (WITH  ATMOSPHERIC  TRANSMITTANCE 

CORRECTION) 


CORRECTED  BACKGROUND  RADIANCE 


14:06:33.8  - STACK  FLARE  (GAS  BURN-OFF) 


u-U3sHS*zuo/n)  BONHiayy 


197 


WKWW.f  ^ 


BACKGROUND  RADIANCE 


GARY,  INDIANA  (816/4)  - 14:15:42.2  - BACKGROUND 


BACKGROUND  RADIANCE 


i-uo*bs*suo/m  BONUiaya 


201 


.1. 


WYANDOTTE,  MICHIGAN  (816/6)  - 15:02:25.0  - BLAST  FURNACE 


5090.00 


(816/6) 

WYANDOTTE. MICHIGAN  - 15:02*25.0  - BLAST  FURNACE  (WITH  ATMOSPHERIC  TRANSMITTANCE 

ORRECTION) 


CORRECTED  TARGET  RADIANCE 


MICHIGAN  - (816/6)  - 15:03:28.3  - BACKGROUND 


816/6 


(816/7) 

WYANDOTTE, MICHIGAN  - 15*. 07s 42. 7 - COKE  OVEN  (WITH  ATMOSPHERIC  TRANSMITTANCE 

CORRECTION) 


WYANDOTTE,  MICHIGAN  - (816/7)  - 15:10:34.5  - BACKGROUND 


BRCKGROUND  RADIANCE 


U-U3*JS*2U0/f1)  33NbIQUM 


•217 


A/  j 


lAli* 


*v  : ^ ~W_ 
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STACKS 


WYANDOTTE,  MICHIGAN  - (816/8)  - 15:27:34.3  - STACKS  WITH  RR  CARS 


WYANDOTTE.  MICHIGAN  - (816/8)  - 15:29:14.8  - BACKGROUND 


BACKGROUND  RADIANCE 


u-u3*as*2U3/n)33Nuiaya 


WYANDOTTE,  MICHIGAN  - (816/8)  - 15;?9:14.8  - BACKGROUND 


BACKGROUND  RADIANCE 


CORRECTED  TARGET  RADIANCE 


WYANDOTTE,  MICHIGAN  - (816/10)  - 15:43:14.3  - BACKGROUND 


2'3‘j 


BACKGROUND  RRDIRNCE 


WYANDOTTE,  MICHIGAN  - (816/10)  - 15:43:14.8  - BACKGROUND 


BACKGROUND  RADIANCE 


■mm 


ASHLAND,  KENTUCKY  - (815/2)  - 11:33:26,5  - FLAME  FROM  STACK 


TARGET  RADIANCE 


ASHLAND,  KENTUCKY  - (815/2)  - 11:33:26.5  - FLAME  FROM  STACK 


TARGET  RADIANCE 


ASHLAND  KENTUCKY  - (815/2)  - 11:33:26.5  - FLAME  FROM  STACK 


TRRGET  RRDIRNCE 


(815/2 

ASHLANO, KENTUCKY  - 11 


241 


. ■ v»  fc-  1 v&aa'SAt 


CORRECTED  TARGET  R* 


ASHLAND,  KENTUCKY  - (815/2)  - 11,36:31. 6 - BACKGROUND 


BACKGROUND  RADIANCE 


ASHLAND,  KENTUCKY  - (815/2)  - 11:36:31.6  - BACKGROUND 


— £*,  . .Ltk. 


BACKGROUND  RADIANCE 


ASHLAND  KENTUCKY  - (815/2)  - 11:36:31. 6 - BACKGROUND 


BACKGROUND  RADIANCE 


ASHLAND. KENTUCKY It  36*31.6  - BACKGROUND  (WITH  ATMOSPHERIC  TRANSMITTANCE 

CORRECTION) 


TRRGET  RADIANCE 


mmb 


I 


o 


(l~U0*HS*SU3/n)33NdiabH 


FAIRFIELD  ALABAMA  - (815/41  - 12:57:26.2  - TIP  OF  STACK 


TARGET  RADIANCE 


CORRECTED  TARGET  RADIANCE 


FA IRFiELD,  ALABAMA  - (815/4)  - 12:58:51,6  - BACKGROUND 


BRCKGROUND  RRDIRNCE 


FAIRFIELD, ALABAMA  - 12*58:51.6  - BACKGROUND  (WITH  ATMOSPHERIC  TRANSMITTANCE 

CORRECTION) 


257 


CORRECTED  BACKGROUND  RADIANCE 


FAIRFIELD,  ALABAMA  (815/4)  - 13:07:5.9  - ORANGE  HOT  METAL 


TARGET  RADIANCE 


FAIRFIELD  ALABAMA  - (815/4)  - 13:07:5.9  - ORANGE  HOT  METAL 


TARGET  RADIANCE 


CORRECTED  TARGET  RADIANCE 


(815/2  *) 

ASHLAND, KENTUCKY  - 11*36*  11.8  - BACKGROUND  (WITH  ATMOSPHERIC  TRANSMITTANCE 

CORRECTION) 


CORRECTED  BACKGROUND  RADIANCE 


ASHLAND,  KENTUCKY  - (815/2  *)  - 11:43:49.2  - REFINERY  COMPLEX 


5090, 


flSHLfiNH. KENTUCKY  - Ilt43i49.2  - REFINERY  COMPLEX  (WITH  ATMOSPHERIC  TRANS 

PITTANCE  CORRECTION) 


CORRECTED  TARGET  RADIANCE 


APPENDIX  E 

GRAY BODY  TARGET  TEMPERATURE  DETERMINATION 


For  a radiating,  absorbing  (nonscattering)  atmosphere,  the  spectral 
radiance  Lv(r)  at  wavenumber  V and  optical  depth  t is  given  by 


L v (x)  * L v(0)  e 


B u ( t 1 ) eT  dr'  (E— 1 ) 


where  B^(t)  is  the  Planck  radiation  function  at  optical  depth  t and  L^(0)  is  the 
radiance  at  r * 0.  The  Planck  function  is  given  by 

3 


B.,  = 


av 


V bV  , 
e - 1 


(E— 2) 


where 


and 


—12  2 —4 

a = 1,191062  x 10  Watts/cm  -sr-cm 


, 1.438786 

0 — cmn 


where  u is  the  temperature. 

If  the  radiation  function  B^Ct)  is  independent  of  optical  thickness  then 
Equation  (E— 1)  becomes 


Lv(0  = Lv(0)e~T  + Bv(l  - e T) 


(E— 3) 


or 


Lv<V  * Lv(0)Tv  * 8v(1  - Iv) 


where  T^  is  the  spectral  transmittance. 
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(E— 4) 


□sseas 


If  the  surface  also  emits  radiation  thermally  then  integrating 
Equation  (E— 4)  over  a nairow  spectral  band  gives 


evBvsTv 


\>a(1  ~ V 


(E-5) 


where  e,.  is  the  emissivity  of  the  surface  and  Bv.  and  B , are  the  Planck  fune- 
v ' vs  Va 

tions  for  the  surface  and  atmosphere  respectively. 


Knowing  the  value  of  experimentally  one  mu3t  subtract  the  path 

radiance  and  divide  by  the  transmittance  in  Equation  (E— 5)  in  order  to  find  the 

surfcce  radiance  . 

v vs 

Algorithm  for  Temperature  Determination:  The  spectral  radiance  for  the 
surface  is  given  by 


L V a £vBv 


( E — 6 ) 


where  we  do  not  know  either  the  emissivity  E^or  the  temperature  u.  Nevertheless, 
if  we  know  the  radiance  Lv  and  assume  a grey  surface,  i.e.,  * E we  can  find  the 
temperature.  The  derivative  of  Equation  ( E — 6 ) is 


L\ 


3a 


2 3 , bv 

V - av  be 


bV 

e - 1 


( E — 7 ) 


Dividing  Equation  ( E — 7 ) by  (E— 6)  gives 


VL' 

LV 


3 


xe 


x , 
e — l 


(E— 8) 


where  x * bv.  Since  the  left  hand  side  of  Equation  (E— 8)  is  known,  the  right  hand 
side  can  be  solved  for  x.  The  temperature  is  then  given  by 


_i , 


u(K) 


] .438786 v( cm  *) 


(E— 9) 


where  x is  the  solution  of  Equation  (E-8). 
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I i 


g-^fufggff 


A 

i 


~Z  ,5 


Once  the  effective  temperature  has  been  determined  using  the  above 
algorithm,  the  emissivity  can  then  be  determined  by  inverting  Equation  (E-6), 


£ 


BV 


(E— 10) 


APPENDIX  F 

SMOKE  CLOUD  EXTINCTION  MODEL 


32 

Casperson  has  developed  an  analytical  approximation  to  the  exact 
Lorentz-Mie  treatment  of  attenuation  of  light  beams  propagating  through  ensembles 
of  spherical  particles  with  complex  index  of  refraction.  The  general  form  for  his 
exponential  extinction  coefficient  is 

Y=  (rra2  N /n2)  [2  (n  + 2)(n  + 1)  + Y - Y sin  Y,  -Y.  sin  Y-] 

max  o 1 l j A 5 


where  N is  the  particle  number  density,  a is  the  radius  at  which  the  particle 
o v max 

distribution  function  peaks,  and  n is  the  index  for  the  power-law-exponential 
particle  size  distribution  function.  A typical  value  of  n is  4.6. 


The  can  be  expressed  in  terms  of  a normalized  wavelength 
absorption  coefficient  4»  yielding 

.,2 


(nk)  * and  an 


Yj  (X)  - 


4n2  X" 


(1  + x'V) 


(1  + X'  2 <l>2)  J 


-1 


Y - ^nXn  t liAl 

Ci  + rV>1/2 


2^2 
n X 


(n  + 2)/2 


L n2  X'2  (1  + <|>/ n ) 2 + U 


Y3  (X) 


* (n  + 2)  tan 


-1 


1 


L nX'  (1  + $/n)  J 


-tan  ^(X*  <J>) 


\ (X) 


y5  (X) 


> 2 
~>n  X 

(1  + x" 2 4>2) 


(n  + 1)  tan 


2 x'2 


[n2  X'  2 (1  + <fr/n)2  +1 

1 A 1 -2  tan-1  (x'/t) 

|_n  X"  (1  + <t>/n) 


(n  + 1 ) / 2 


The  absorption  coefficient  (f>  is  related  to  the  exponential  amplitude  absorption 
coefficients  , outside  the  sphere  and  inside  the  sphere  by 


111 


2a  (a„  - 
max  2 


al} 


The  coefficient  K is  related  to  the  index  of  refraction  of  the  particle 
the  index  of  refraction  of  the  surrounding  medium  n^,  and  the  vacuum  wave- 
length X by 


nK  = 4 it  (n„  - n. ) a /X 
2 1 max 


Dalzell  & Sarofim  have  developed  explicit  expressions  for  n and  k that 
are  appropriate  for  carbon  smokes.  Their  dispersion  model  is  based  upon  fitting 
classical  electron  theory  to  room  temperature  measurement  of  n and  a. 


1 - 


Fce2/meo 


2 2 
gc  + w 


- E 

j 


F.  (W? 
_J J_ 


W2)  e2/m  £ 


(to2  - to2)  + (02  g2 


2na  * 


F g e /m  e 

c 6C  ( 


(g  I +oi2) 


♦ E 

j 


F.  (Ogj  e /m  £q 

o o o 
(Wj  - wr+u)  gj 


In  these  equations  F is  the  number  of  effective  electrons  per  unit 

volume,  q is  the  electron  damping  constant,  to  is  the  frequency  of  the  radiation, 

th  2 

j is  the  natural  frequency  of  the  j electron,  and  e /m  eQ  is  the  square  of  the 
electronic  charge  divided  by  the  mass  of  the  electron  and  the  electric  permit- 
tivity. Three  electrons,  two  bound  (subscript  j)  and  one  free  (subscript  c)  were 
employed  to  fit  the  experimental  data.  The  constants  used  in  the  final  dispersion 
model  for  propane  soot  are  presented  in  Table  F.l. 
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Table  F.l 


CONSTANTS  IN  THE  DISPERSION  MODEL 

Electron 

electrons  /m 

q , sec  ^ 

-1 

tij,  sec 

C 

4. 06xl027 

6 . OOxlO15 

* * « 

1 

2.6'xlO27 

6.00xl015 

1.25xl015 

2 

2 . 86xl028 

7.25xl015 

7. 25xl015 
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